docs/examples/driver_examples/Qcodes example with Lakeshore 325.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# QCoDeS Example with Lakeshore 325\n",
    "\n",
    "Here provided is an example session with model 325 of the Lakeshore temperature controller "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib notebook"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "\n",
    "from qcodes.instrument_drivers.Lakeshore import LakeshoreModel325"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Connected to: LSCI 325 (serial:LSA2251, firmware:1.8/1.1) in 1.30s\n"
     ]
    }
   ],
   "source": [
    "lake = LakeshoreModel325(\"lake\", \"GPIB0::12::INSTR\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Sensor commands"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'OK'"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Check that the sensor is in the correct status\n",
    "lake.sensor_A.status()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "23.102"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# What temperature is it reading?\n",
    "lake.sensor_A.temperature()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'K'"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "lake.sensor_A.temperature.unit"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "# We can access the sensor objects through the sensor list as well\n",
    "assert lake.sensor_A is lake.sensor[0]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Heater commands "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'A'"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# In a closed loop configuration, heater 1 reads from...\n",
    "lake.heater_1.input_channel()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'Kelvin'"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "lake.heater_1.unit()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "P =  400.0\n",
      "I =  40.0\n",
      "D =  10.0\n"
     ]
    }
   ],
   "source": [
    "# Get the PID values\n",
    "print(\"P = \", lake.heater_1.P())\n",
    "print(\"I = \", lake.heater_1.I())\n",
    "print(\"D = \", lake.heater_1.D())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'Off'"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Is the heater on?\n",
    "lake.heater_1.output_range()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Loading and updating sensor calibration values"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "curve = lake.sensor_A.curve"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "curve_data = curve.get_data()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "dict_keys(['Temperature (K)', 'log Ohm'])"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "curve_data.keys()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": "/* Put everything inside the global mpl namespace */\nwindow.mpl = {};\n\n\nmpl.get_websocket_type = function() {\n    if (typeof(WebSocket) !== 'undefined') {\n        return WebSocket;\n    } else if (typeof(MozWebSocket) !== 'undefined') {\n        return MozWebSocket;\n    } else {\n        alert('Your browser does not have WebSocket support.' +\n              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n              'Firefox 4 and 5 are also supported but you ' +\n              'have to enable WebSockets in about:config.');\n    };\n}\n\nmpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n    this.id = figure_id;\n\n    this.ws = websocket;\n\n    this.supports_binary = (this.ws.binaryType != undefined);\n\n    if (!this.supports_binary) {\n        var warnings = document.getElementById(\"mpl-warnings\");\n        if (warnings) {\n            warnings.style.display = 'block';\n            warnings.textContent = (\n                \"This browser does not support binary websocket messages. \" +\n                    \"Performance may be slow.\");\n        }\n    }\n\n    this.imageObj = new Image();\n\n    this.context = undefined;\n    this.message = undefined;\n    this.canvas = undefined;\n    this.rubberband_canvas = undefined;\n    this.rubberband_context = undefined;\n    this.format_dropdown = undefined;\n\n    this.image_mode = 'full';\n\n    this.root = $('<div/>');\n    this._root_extra_style(this.root)\n    this.root.attr('style', 'display: inline-block');\n\n    $(parent_element).append(this.root);\n\n    this._init_header(this);\n    this._init_canvas(this);\n    this._init_toolbar(this);\n\n    var fig = this;\n\n    this.waiting = false;\n\n    this.ws.onopen =  function () {\n            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n            fig.send_message(\"send_image_mode\", {});\n            if (mpl.ratio != 1) {\n                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n            }\n            fig.send_message(\"refresh\", {});\n        }\n\n    this.imageObj.onload = function() {\n            if (fig.image_mode == 'full') {\n                // Full images could contain transparency (where diff images\n                // almost always do), so we need to clear the canvas so that\n                // there is no ghosting.\n                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n            }\n            fig.context.drawImage(fig.imageObj, 0, 0);\n        };\n\n    this.imageObj.onunload = function() {\n        this.ws.close();\n    }\n\n    this.ws.onmessage = this._make_on_message_function(this);\n\n    this.ondownload = ondownload;\n}\n\nmpl.figure.prototype._init_header = function() {\n    var titlebar = $(\n        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n        'ui-helper-clearfix\"/>');\n    var titletext = $(\n        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n        'text-align: center; padding: 3px;\"/>');\n    titlebar.append(titletext)\n    this.root.append(titlebar);\n    this.header = titletext[0];\n}\n\n\n\nmpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n\n}\n\n\nmpl.figure.prototype._root_extra_style = function(canvas_div) {\n\n}\n\nmpl.figure.prototype._init_canvas = function() {\n    var fig = this;\n\n    var canvas_div = $('<div/>');\n\n    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n\n    function canvas_keyboard_event(event) {\n        return fig.key_event(event, event['data']);\n    }\n\n    canvas_div.keydown('key_press', canvas_keyboard_event);\n    canvas_div.keyup('key_release', canvas_keyboard_event);\n    this.canvas_div = canvas_div\n    this._canvas_extra_style(canvas_div)\n    this.root.append(canvas_div);\n\n    var canvas = $('<canvas/>');\n    canvas.addClass('mpl-canvas');\n    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n\n    this.canvas = canvas[0];\n    this.context = canvas[0].getContext(\"2d\");\n\n    var backingStore = this.context.backingStorePixelRatio ||\n\tthis.context.webkitBackingStorePixelRatio ||\n\tthis.context.mozBackingStorePixelRatio ||\n\tthis.context.msBackingStorePixelRatio ||\n\tthis.context.oBackingStorePixelRatio ||\n\tthis.context.backingStorePixelRatio || 1;\n\n    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n\n    var rubberband = $('<canvas/>');\n    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n\n    var pass_mouse_events = true;\n\n    canvas_div.resizable({\n        start: function(event, ui) {\n            pass_mouse_events = false;\n        },\n        resize: function(event, ui) {\n            fig.request_resize(ui.size.width, ui.size.height);\n        },\n        stop: function(event, ui) {\n            pass_mouse_events = true;\n            fig.request_resize(ui.size.width, ui.size.height);\n        },\n    });\n\n    function mouse_event_fn(event) {\n        if (pass_mouse_events)\n            return fig.mouse_event(event, event['data']);\n    }\n\n    rubberband.mousedown('button_press', mouse_event_fn);\n    rubberband.mouseup('button_release', mouse_event_fn);\n    // Throttle sequential mouse events to 1 every 20ms.\n    rubberband.mousemove('motion_notify', mouse_event_fn);\n\n    rubberband.mouseenter('figure_enter', mouse_event_fn);\n    rubberband.mouseleave('figure_leave', mouse_event_fn);\n\n    canvas_div.on(\"wheel\", function (event) {\n        event = event.originalEvent;\n        event['data'] = 'scroll'\n        if (event.deltaY < 0) {\n            event.step = 1;\n        } else {\n            event.step = -1;\n        }\n        mouse_event_fn(event);\n    });\n\n    canvas_div.append(canvas);\n    canvas_div.append(rubberband);\n\n    this.rubberband = rubberband;\n    this.rubberband_canvas = rubberband[0];\n    this.rubberband_context = rubberband[0].getContext(\"2d\");\n    this.rubberband_context.strokeStyle = \"#000000\";\n\n    this._resize_canvas = function(width, height) {\n        // Keep the size of the canvas, canvas container, and rubber band\n        // canvas in synch.\n        canvas_div.css('width', width)\n        canvas_div.css('height', height)\n\n        canvas.attr('width', width * mpl.ratio);\n        canvas.attr('height', height * mpl.ratio);\n        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n\n        rubberband.attr('width', width);\n        rubberband.attr('height', height);\n    }\n\n    // Set the figure to an initial 600x600px, this will subsequently be updated\n    // upon first draw.\n    this._resize_canvas(600, 600);\n\n    // Disable right mouse context menu.\n    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n        return false;\n    });\n\n    function set_focus () {\n        canvas.focus();\n        canvas_div.focus();\n    }\n\n    window.setTimeout(set_focus, 100);\n}\n\nmpl.figure.prototype._init_toolbar = function() {\n    var fig = this;\n\n    var nav_element = $('<div/>')\n    nav_element.attr('style', 'width: 100%');\n    this.root.append(nav_element);\n\n    // Define a callback function for later on.\n    function toolbar_event(event) {\n        return fig.toolbar_button_onclick(event['data']);\n    }\n    function toolbar_mouse_event(event) {\n        return fig.toolbar_button_onmouseover(event['data']);\n    }\n\n    for(var toolbar_ind in mpl.toolbar_items) {\n        var name = mpl.toolbar_items[toolbar_ind][0];\n        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n        var image = mpl.toolbar_items[toolbar_ind][2];\n        var method_name = mpl.toolbar_items[toolbar_ind][3];\n\n        if (!name) {\n            // put a spacer in here.\n            continue;\n        }\n        var button = $('<button/>');\n        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n                        'ui-button-icon-only');\n        button.attr('role', 'button');\n        button.attr('aria-disabled', 'false');\n        button.click(method_name, toolbar_event);\n        button.mouseover(tooltip, toolbar_mouse_event);\n\n        var icon_img = $('<span/>');\n        icon_img.addClass('ui-button-icon-primary ui-icon');\n        icon_img.addClass(image);\n        icon_img.addClass('ui-corner-all');\n\n        var tooltip_span = $('<span/>');\n        tooltip_span.addClass('ui-button-text');\n        tooltip_span.html(tooltip);\n\n        button.append(icon_img);\n        button.append(tooltip_span);\n\n        nav_element.append(button);\n    }\n\n    var fmt_picker_span = $('<span/>');\n\n    var fmt_picker = $('<select/>');\n    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n    fmt_picker_span.append(fmt_picker);\n    nav_element.append(fmt_picker_span);\n    this.format_dropdown = fmt_picker[0];\n\n    for (var ind in mpl.extensions) {\n        var fmt = mpl.extensions[ind];\n        var option = $(\n            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n        fmt_picker.append(option)\n    }\n\n    // Add hover states to the ui-buttons\n    $( \".ui-button\" ).hover(\n        function() { $(this).addClass(\"ui-state-hover\");},\n        function() { $(this).removeClass(\"ui-state-hover\");}\n    );\n\n    var status_bar = $('<span class=\"mpl-message\"/>');\n    nav_element.append(status_bar);\n    this.message = status_bar[0];\n}\n\nmpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n    // which will in turn request a refresh of the image.\n    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n}\n\nmpl.figure.prototype.send_message = function(type, properties) {\n    properties['type'] = type;\n    properties['figure_id'] = this.id;\n    this.ws.send(JSON.stringify(properties));\n}\n\nmpl.figure.prototype.send_draw_message = function() {\n    if (!this.waiting) {\n        this.waiting = true;\n        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n    }\n}\n\n\nmpl.figure.prototype.handle_save = function(fig, msg) {\n    var format_dropdown = fig.format_dropdown;\n    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n    fig.ondownload(fig, format);\n}\n\n\nmpl.figure.prototype.handle_resize = function(fig, msg) {\n    var size = msg['size'];\n    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n        fig._resize_canvas(size[0], size[1]);\n        fig.send_message(\"refresh\", {});\n    };\n}\n\nmpl.figure.prototype.handle_rubberband = function(fig, msg) {\n    var x0 = msg['x0'] / mpl.ratio;\n    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n    var x1 = msg['x1'] / mpl.ratio;\n    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n    x0 = Math.floor(x0) + 0.5;\n    y0 = Math.floor(y0) + 0.5;\n    x1 = Math.floor(x1) + 0.5;\n    y1 = Math.floor(y1) + 0.5;\n    var min_x = Math.min(x0, x1);\n    var min_y = Math.min(y0, y1);\n    var width = Math.abs(x1 - x0);\n    var height = Math.abs(y1 - y0);\n\n    fig.rubberband_context.clearRect(\n        0, 0, fig.canvas.width, fig.canvas.height);\n\n    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n}\n\nmpl.figure.prototype.handle_figure_label = function(fig, msg) {\n    // Updates the figure title.\n    fig.header.textContent = msg['label'];\n}\n\nmpl.figure.prototype.handle_cursor = function(fig, msg) {\n    var cursor = msg['cursor'];\n    switch(cursor)\n    {\n    case 0:\n        cursor = 'pointer';\n        break;\n    case 1:\n        cursor = 'default';\n        break;\n    case 2:\n        cursor = 'crosshair';\n        break;\n    case 3:\n        cursor = 'move';\n        break;\n    }\n    fig.rubberband_canvas.style.cursor = cursor;\n}\n\nmpl.figure.prototype.handle_message = function(fig, msg) {\n    fig.message.textContent = msg['message'];\n}\n\nmpl.figure.prototype.handle_draw = function(fig, msg) {\n    // Request the server to send over a new figure.\n    fig.send_draw_message();\n}\n\nmpl.figure.prototype.handle_image_mode = function(fig, msg) {\n    fig.image_mode = msg['mode'];\n}\n\nmpl.figure.prototype.updated_canvas_event = function() {\n    // Called whenever the canvas gets updated.\n    this.send_message(\"ack\", {});\n}\n\n// A function to construct a web socket function for onmessage handling.\n// Called in the figure constructor.\nmpl.figure.prototype._make_on_message_function = function(fig) {\n    return function socket_on_message(evt) {\n        if (evt.data instanceof Blob) {\n            /* FIXME: We get \"Resource interpreted as Image but\n             * transferred with MIME type text/plain:\" errors on\n             * Chrome.  But how to set the MIME type?  It doesn't seem\n             * to be part of the websocket stream */\n            evt.data.type = \"image/png\";\n\n            /* Free the memory for the previous frames */\n            if (fig.imageObj.src) {\n                (window.URL || window.webkitURL).revokeObjectURL(\n                    fig.imageObj.src);\n            }\n\n            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n                evt.data);\n            fig.updated_canvas_event();\n            fig.waiting = false;\n            return;\n        }\n        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n            fig.imageObj.src = evt.data;\n            fig.updated_canvas_event();\n            fig.waiting = false;\n            return;\n        }\n\n        var msg = JSON.parse(evt.data);\n        var msg_type = msg['type'];\n\n        // Call the  \"handle_{type}\" callback, which takes\n        // the figure and JSON message as its only arguments.\n        try {\n            var callback = fig[\"handle_\" + msg_type];\n        } catch (e) {\n            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n            return;\n        }\n\n        if (callback) {\n            try {\n                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n                callback(fig, msg);\n            } catch (e) {\n                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n            }\n        }\n    };\n}\n\n// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\nmpl.findpos = function(e) {\n    //this section is from http://www.quirksmode.org/js/events_properties.html\n    var targ;\n    if (!e)\n        e = window.event;\n    if (e.target)\n        targ = e.target;\n    else if (e.srcElement)\n        targ = e.srcElement;\n    if (targ.nodeType == 3) // defeat Safari bug\n        targ = targ.parentNode;\n\n    // jQuery normalizes the pageX and pageY\n    // pageX,Y are the mouse positions relative to the document\n    // offset() returns the position of the element relative to the document\n    var x = e.pageX - $(targ).offset().left;\n    var y = e.pageY - $(targ).offset().top;\n\n    return {\"x\": x, \"y\": y};\n};\n\n/*\n * return a copy of an object with only non-object keys\n * we need this to avoid circular references\n * http://stackoverflow.com/a/24161582/3208463\n */\nfunction simpleKeys (original) {\n  return Object.keys(original).reduce(function (obj, key) {\n    if (typeof original[key] !== 'object')\n        obj[key] = original[key]\n    return obj;\n  }, {});\n}\n\nmpl.figure.prototype.mouse_event = function(event, name) {\n    var canvas_pos = mpl.findpos(event)\n\n    if (name === 'button_press')\n    {\n        this.canvas.focus();\n        this.canvas_div.focus();\n    }\n\n    var x = canvas_pos.x * mpl.ratio;\n    var y = canvas_pos.y * mpl.ratio;\n\n    this.send_message(name, {x: x, y: y, button: event.button,\n                             step: event.step,\n                             guiEvent: simpleKeys(event)});\n\n    /* This prevents the web browser from automatically changing to\n     * the text insertion cursor when the button is pressed.  We want\n     * to control all of the cursor setting manually through the\n     * 'cursor' event from matplotlib */\n    event.preventDefault();\n    return false;\n}\n\nmpl.figure.prototype._key_event_extra = function(event, name) {\n    // Handle any extra behaviour associated with a key event\n}\n\nmpl.figure.prototype.key_event = function(event, name) {\n\n    // Prevent repeat events\n    if (name == 'key_press')\n    {\n        if (event.which === this._key)\n            return;\n        else\n            this._key = event.which;\n    }\n    if (name == 'key_release')\n        this._key = null;\n\n    var value = '';\n    if (event.ctrlKey && event.which != 17)\n        value += \"ctrl+\";\n    if (event.altKey && event.which != 18)\n        value += \"alt+\";\n    if (event.shiftKey && event.which != 16)\n        value += \"shift+\";\n\n    value += 'k';\n    value += event.which.toString();\n\n    this._key_event_extra(event, name);\n\n    this.send_message(name, {key: value,\n                             guiEvent: simpleKeys(event)});\n    return false;\n}\n\nmpl.figure.prototype.toolbar_button_onclick = function(name) {\n    if (name == 'download') {\n        this.handle_save(this, null);\n    } else {\n        this.send_message(\"toolbar_button\", {name: name});\n    }\n};\n\nmpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n    this.message.textContent = tooltip;\n};\nmpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n\nmpl.extensions = [\"eps\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\"];\n\nmpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n    // Create a \"websocket\"-like object which calls the given IPython comm\n    // object with the appropriate methods. Currently this is a non binary\n    // socket, so there is still some room for performance tuning.\n    var ws = {};\n\n    ws.close = function() {\n        comm.close()\n    };\n    ws.send = function(m) {\n        //console.log('sending', m);\n        comm.send(m);\n    };\n    // Register the callback with on_msg.\n    comm.on_msg(function(msg) {\n        //console.log('receiving', msg['content']['data'], msg);\n        // Pass the mpl event to the overriden (by mpl) onmessage function.\n        ws.onmessage(msg['content']['data'])\n    });\n    return ws;\n}\n\nmpl.mpl_figure_comm = function(comm, msg) {\n    // This is the function which gets called when the mpl process\n    // starts-up an IPython Comm through the \"matplotlib\" channel.\n\n    var id = msg.content.data.id;\n    // Get hold of the div created by the display call when the Comm\n    // socket was opened in Python.\n    var element = $(\"#\" + id);\n    var ws_proxy = comm_websocket_adapter(comm)\n\n    function ondownload(figure, format) {\n        window.open(figure.imageObj.src);\n    }\n\n    var fig = new mpl.figure(id, ws_proxy,\n                           ondownload,\n                           element.get(0));\n\n    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n    // web socket which is closed, not our websocket->open comm proxy.\n    ws_proxy.onopen();\n\n    fig.parent_element = element.get(0);\n    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n    if (!fig.cell_info) {\n        console.error(\"Failed to find cell for figure\", id, fig);\n        return;\n    }\n\n    var output_index = fig.cell_info[2]\n    var cell = fig.cell_info[0];\n\n};\n\nmpl.figure.prototype.handle_close = function(fig, msg) {\n    var width = fig.canvas.width/mpl.ratio\n    fig.root.unbind('remove')\n\n    // Update the output cell to use the data from the current canvas.\n    fig.push_to_output();\n    var dataURL = fig.canvas.toDataURL();\n    // Re-enable the keyboard manager in IPython - without this line, in FF,\n    // the notebook keyboard shortcuts fail.\n    IPython.keyboard_manager.enable()\n    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n    fig.close_ws(fig, msg);\n}\n\nmpl.figure.prototype.close_ws = function(fig, msg){\n    fig.send_message('closing', msg);\n    // fig.ws.close()\n}\n\nmpl.figure.prototype.push_to_output = function(remove_interactive) {\n    // Turn the data on the canvas into data in the output cell.\n    var width = this.canvas.width/mpl.ratio\n    var dataURL = this.canvas.toDataURL();\n    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n}\n\nmpl.figure.prototype.updated_canvas_event = function() {\n    // Tell IPython that the notebook contents must change.\n    IPython.notebook.set_dirty(true);\n    this.send_message(\"ack\", {});\n    var fig = this;\n    // Wait a second, then push the new image to the DOM so\n    // that it is saved nicely (might be nice to debounce this).\n    setTimeout(function () { fig.push_to_output() }, 1000);\n}\n\nmpl.figure.prototype._init_toolbar = function() {\n    var fig = this;\n\n    var nav_element = $('<div/>')\n    nav_element.attr('style', 'width: 100%');\n    this.root.append(nav_element);\n\n    // Define a callback function for later on.\n    function toolbar_event(event) {\n        return fig.toolbar_button_onclick(event['data']);\n    }\n    function toolbar_mouse_event(event) {\n        return fig.toolbar_button_onmouseover(event['data']);\n    }\n\n    for(var toolbar_ind in mpl.toolbar_items){\n        var name = mpl.toolbar_items[toolbar_ind][0];\n        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n        var image = mpl.toolbar_items[toolbar_ind][2];\n        var method_name = mpl.toolbar_items[toolbar_ind][3];\n\n        if (!name) { continue; };\n\n        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n        button.click(method_name, toolbar_event);\n        button.mouseover(tooltip, toolbar_mouse_event);\n        nav_element.append(button);\n    }\n\n    // Add the status bar.\n    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n    nav_element.append(status_bar);\n    this.message = status_bar[0];\n\n    // Add the close button to the window.\n    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n    button.click(function (evt) { fig.handle_close(fig, {}); } );\n    button.mouseover('Stop Interaction', toolbar_mouse_event);\n    buttongrp.append(button);\n    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n    titlebar.prepend(buttongrp);\n}\n\nmpl.figure.prototype._root_extra_style = function(el){\n    var fig = this\n    el.on(\"remove\", function(){\n\tfig.close_ws(fig, {});\n    });\n}\n\nmpl.figure.prototype._canvas_extra_style = function(el){\n    // this is important to make the div 'focusable\n    el.attr('tabindex', 0)\n    // reach out to IPython and tell the keyboard manager to turn it's self\n    // off when our div gets focus\n\n    // location in version 3\n    if (IPython.notebook.keyboard_manager) {\n        IPython.notebook.keyboard_manager.register_events(el);\n    }\n    else {\n        // location in version 2\n        IPython.keyboard_manager.register_events(el);\n    }\n\n}\n\nmpl.figure.prototype._key_event_extra = function(event, name) {\n    var manager = IPython.notebook.keyboard_manager;\n    if (!manager)\n        manager = IPython.keyboard_manager;\n\n    // Check for shift+enter\n    if (event.shiftKey && event.which == 13) {\n        this.canvas_div.blur();\n        // select the cell after this one\n        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n        IPython.notebook.select(index + 1);\n    }\n}\n\nmpl.figure.prototype.handle_save = function(fig, msg) {\n    fig.ondownload(fig, null);\n}\n\n\nmpl.find_output_cell = function(html_output) {\n    // Return the cell and output element which can be found *uniquely* in the notebook.\n    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n    // IPython event is triggered only after the cells have been serialised, which for\n    // our purposes (turning an active figure into a static one), is too late.\n    var cells = IPython.notebook.get_cells();\n    var ncells = cells.length;\n    for (var i=0; i<ncells; i++) {\n        var cell = cells[i];\n        if (cell.cell_type === 'code'){\n            for (var j=0; j<cell.output_area.outputs.length; j++) {\n                var data = cell.output_area.outputs[j];\n                if (data.data) {\n                    // IPython >= 3 moved mimebundle to data attribute of output\n                    data = data.data;\n                }\n                if (data['text/html'] == html_output) {\n                    return [cell, data, j];\n                }\n            }\n        }\n    }\n}\n\n// Register the function which deals with the matplotlib target/channel.\n// The kernel may be null if the page has been refreshed.\nif (IPython.notebook.kernel != null) {\n    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n}\n",
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAACebSURBVHhe7d0LlKRleSfwGklABaYjQQJyCcga8YIwK0YhgJcEFXSVS05AcnaDGjUmGt0rTMztJMJoNkcxXmIkHg3egueE4CoiGFkDMoio3AyChh2R4ZYQcZqgQoKzz1NV30x1dd26q3q66nt/v5y/9X1VPcN0pbr73+/7fm81AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAqJc17VuWJ5+/J0QeaJ4BALNi98hdka3Ns8IogOPZN7K5dQgAzJj9Ine2DsuiAI5nbWTLHXfc0Vi7Ng8BgGk3Pz/f2H///fNwLk/zoDQK4HiaBTAogAAwI7IAzs1l9yu3AD6qfQsAQCEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIBT7O4tP2xsvO2+5i0AwKTUtQCuj2yNnNs86++UyM2Rh9q3J0WmwgXXfrfxC2+7vHH6edc0b/McAGAS6lgAnxV5beTG5ll/R0YuiHwkclj79pORZ0dWVY74rb/wpsaPs8KGvP2dC79hJBAAmIi6FcDdIh+LvCZyf94xwJsjn49siNzSvv1CJO9fVZvue3Bb+as8snVr4zv3/aB9BgCwfHUrgO+NXBz5u+bZYDkCeFnrcJtLI0e1DnvaJbK2I7tHJu6gPXdtPGpN+6RtpzVrGgfu+dj2GQDA8tWpAJ4WeWYk1/+NYu/Iva3DbfI87+8n/+4tHdkcmbh95h7T2HDyoc3Sl/L2nJOf3rwfAGBcXeNMM2v/yFcjL4zckHeEL0auj/Sb0n048muRTzTPWn418sHIo5tni+UIYKaSI4Cbt2zZ0li7NgcEJyvX/OW0b478KX8AMBnz8/ONubm5PMz/mc+D0tRlBDBH/vaKfC3y7+08N/Lb7eOdIt3uiXSP9uXf0T0q2CmvFs4XSpUHIismS9+RB/+08gcATFRdCmBevHFo5PCO5IhgXhCSx49Eul0dOa51uE2OIG5sHQIA1FNdCmCOxH2jKw9G/qV9nM6P5JW+lXdFsvCdGTmkfftLkWF7BwIAzLS6XQU8yAGRfVqHTTnSlxeOvDKSewaeETk1ck0EAKC26nIRyGrJKz+2rORFILknYG4LYx0gAEyGi0AUwHGtWAHMt36r3g0k9wTMbWFOfVYOYgIA41AAy5oCnhneCg4AWEkK4BTyVnAAwEpSAKdQr7eCy3NvBQcATIICOIXygo9c89d+J7imrVsbjSu+9c/tMwCA5VMAp9SxP/f4aH3tk5CH1gECAJOgAE6pXAfYtQzQOkAAYCIUwCnVax3gTmvWWAcIAIxNAZxS1TrAqgTm7TknP92G0ADA2BTAKZcXf3TeAgCMSwGcUtVm0FXvcxEIADApCuCUshk0ALBSFMAp5SIQAGClKIBTKi/2OGndvu2zlhPXPcFFIADA2BTAKZVr/f72ujvbZy0XXXeXNYAAwNgUwCllDSAAsFIUwCllDSAAsFIUwCllDSAAsFIUwCllDSAAsFIUwCllDSAAsFIUwCmVawC7lgA21sQd1gACAONSAGdJ14ggAMByKIBTKqeAu/tenpsCBgDGpQBOKdvAAAArRQGcUraBAQBWigI4pWwDAwCsFAVwStkGBgBYKQrglOq1BjDduPn77SMAgOVRAKdUrvU78/hD2mfb/cnnbjUNDACMRQGcYofuO9c+2s40MAAwLgVwink3EABgJdSlAL4+cmNkvp2rI8dH+jkjkpdYdOfRkemW/0oAgDHUpQBujpwVOaKdyyOfijwt0k8WxX268qPI1PBuIADASqhLAfx05LORb7Xzlsi/Rp4T6Se71D1dmSqmgAGAlVDHNYA7RU6L7BrJqeB+dovcHsnRw89E1kWG2SWytiO7R3YsU8AAwJjqVAAPjeSo30OR90dOitwc6eWWSK4DfFnkFZGc+r0q8qTIIOsjWzqS5XHFmAIGAFZCnQrgrZHDIznt++eRv4o8NdLLlyMfjdwQuTLyK5GcOn5jZJANkdybpcp+kRWz6845mLnYY3eu48AtALCj1KlJPBz5x8hXIzlSl+XuTZFR/DhybWTYCGCOLlZXGmceiKyYBx9+pH200A8ezn8uAMDy1HkoKa+fyDV7o8iPzdHDu5tnU8JFIADASqhLATwnckzkwEiuBTw78rzIxyLp/EhO31b+IPKiyBMjWfw+2L7NtYPTzUUgAMCY6lIAfybykUiuA/xC5NmRF0c+H0kHRHKfv8pPRT4Q+Wbkssi+kWMjX4lMDReBAAAroXuGkaXJrWC2hMbatXk4WXdv+WHjqA2XLyiBOQW88awXNPaZe0z7HgBgKebn5xtzc83328//yTX9xXE56awxBQwAjEkBnGL9poA/dNWm1gkAwDIogFMsrwJ+VI9J+r+8YlNzehgAYDkUwCmW6/xeffRB7bPtchdAF4IAAMulAE65lz6j8+Ll7bwbCACwXFrElPNuIADApCmAU877AQMAk6ZFTDkjgADApCmAU84IIAAwaVrElOs3AnjxTXe3jwAAlkYBnHL2AgQAJk0BnHL2AgQAJk0BnAH2AgQAJkmDmAGuBAYAJkkBnAGuBAYAJkmDmAFGAAGASVIAZ4ARQABgkjSIGWAvQABgkhTAGWAvQABgkhTAGWAvQABgkhTAGfGqHgVwzZpG48A9H9s+AwAYjQI4QxbNAm9t3wIALIECOCM23ffgor6X5x+6alPrBABgRArgjOh3IcgHXAgCACyRAjgj8kKQ037+gPbZQl+//f72EQDAcArgDDnyiXu0jxbaai0gALAECuAMOWCP3lf87r/HY9pHAADDKYAzxDuCAACToADOEO8IAgBMggI4Q7wjCAAwCXUpgK+P3BiZb+fqyPGRQU6J3Bx5qH17UmTqvfQZ+7SPFnrszro8ADCaurSGzZGzIke0c3nkU5GnRXo5MnJB5CORw9q3n4w8OzLVrAMEAMZVlwL46chnI99q5y2Rf408J9LLmyOfj2yI3NK+/UIk759qNoQGAMZVx3nDnSKnRXaN5FRwLzkCeFnrcJtLI0e1DqeXDaEBgHHVqQAeGslRv1zT9/5IrunLtX297B25t3W4TZ7n/YPsElnbkd0jO1y/DaGvuu2+9hEAQH91KoC3Rg6P5LTvn0f+KvLUSD/d75+RE6vD3lNjfWRLR3Lt4Q53xIF7NP+x3T5+zR2mgQGAoepUAB+O/GPkq5EsajdE3hTp5Z5I92jfXpHuUcFuuVZwriP7RXa4nAZ+xbN7TwO/+wvfbh8BAPRWxzWAlRwkyynbXnJt4HGtw21eGNnYOuwrp5errWYyD0RWRb9p4E98xSggADBYXQrgOZFjIgdGci3g2ZHnRT4WSedHcvSu8q5IFr4zI4e0b38pcm5kJuQ0cC85h21TaABgkLoUwJ+J5F5+uQ4wt3PJ/fxeHMmtXlLOl3buoJwjfXml8CsjuYH0GZFTI9dEZkJOA7/h+Qe3zxb60rf/uX0EALBYr2sJGF1eCbwlNNauzcMda+Nt9zVOP693Z716/QuaJREAWGh+fr4xN5dL+Zvr+XNJV3HqvAaw9nJT6H4N3p6AAEA/CuAMa14N/PP7t88W+vzNeaEzAMBiCuCMe+MvPql9tNBF19/d+Iu/v619BgCwnQI443IU8DXHHNQ+W2jDJbfYEgYAWEQBrIFXHd27AKa3RwkEAOikANZAjgKe3mct4EXX32UqGABYQAGsiX5rAZOpYACgkwJYEzkKuP6EfFOT3mwLAwBUFMAaed2xBzdOPLzzDU+2sy0MAFBRAGvmzOOf0j5ayLYwAEBFAawZ28IAAMMogDVkWxgAYBAFsIaGbQvzvy9VAgGgZApgTQ3aFua9//c2JRAACqYA1tSwbWGyBLooBADKpADWWG4L81vPP7h9tpiLQgCgTApgzf3PFx0ysAT+j09erwQCQGEUwAJkCey3QfRVt32vceSGyxt/cYXpYAAohQJYiH4bRFc2fPYWF4YAQCEUwEIMuygkuToYAMqgABZk2EUhydXBAFB/CmBhcj3g+uMHjwS+7ZJbGjfccX/7DACoGwWwQK977sGNq9e/oHHi4U9o37PQ1sjL37vRhSEAUFMKYKFyTeC5p60bvE/gZ28xHQwANaQAFq6aEl7TPu9mOhgA6kcBpDklfNFvHdU+Wying09878bG2RffbMNoAKgJBZCmw/Z/XN9tYrIEnnflpsZRGy5XBAGgBhRAtsltYgZNB3cWwQuu/W7rTgBg5iiALFBNB/crgSmL4PoLbzISCAAzSgFkkZwOftsphzYeNaAF/jha4Fs/YzoYAGbRoIEehlsb2RIaa9fmYb1kufvQVZsa512xqTnq10u+gH79mIMarzr6oObWMgAw7ebn5xtzc3N5mP8znwelqUsBXB85OZJXMeSQ1MbImZFbI/2cEflQ63CBbDE/ah0OVesCWMkimKN9F990T/uexRRBAGaFAlifKeDnRt4beU7kuMhPRC6L7BoZJP+fvk9XRi1/xchC97svferQdYEuEAGA2VCXAvjiyIcj/xC5IfLKyAGRZ0YGyd6Sw1qdoYcsgcPWBabqAhGbRwPA9KrrRSDNcd3wvfZtP7tFbo9sjnwmsi5CH6c+64DGVWe9oPHaYw8aOBqYF4jYPBoApldd1gB2ys/pU5HHRY7JO/rI6eL/ELkpkgv43hQ5IXJY5NuRXnZpp7J7ZHPd1wD2MsoFIin/n3HWCYc0Dt13rnHQnrtaHwjAqrMGsJ4FMNcCviRydCRH9kaVo6Ffj1wR+e28o4c/jPxB63C7EgtgZdQimHL6+NVHu1AEgNWlANavAL47cmLk2MimvGOJzovsFzm+ebaYEcA+cs1fTvsOK4EpX3SuGAZgtSiA9VkDmJ3iPZHcCuYFkeWUv/w7Do/c3Tzr7aFIvlCqPBAhjLJ5dCVLovcWBoDVU5cRwPdFTo+8PNK599+WSNUuzo/cGck9A1NO5X45kuv9cvgup33/c+QXIl+JjCL/XO33AVyKakr4L6/Y1Phx+75hjAgCsCMZAaxPAew385jbweT2MOmLke9EcgPo9M5IjhjuHcmieF0k1/hdHRmVAthHFsHv3PeDxo2bv994++duaV4ZPEy+GHMUMa82BoCVogDWpwCuFgVwBEu9UORdpx3eOOLAPYwGArAiFEAFcFwK4BK4YhiAaaAAKoDjUgCXYSlFMF+g9hEEYJIUQAVwXArgGJZ6wYhRQQAmQQFUAMelAE5AFsGv335/4w0fv27oiGAyKgjAOBRABXBcCuAEXXDtdxvrL7xppCuGK0YFAVgqBVABHJcCOGFLnRauGBUEYFQKoAI4LgVwhWQRXOo+gpUcFTzzeGUQgN4UQAVwXArgDrDcUcHUOUWcNt33oFIIUDgFUAEclwK4A40zKpjyxZ5/xAghQNkUQAVwXArgKhlnVLCTi0gAyqMAKoDjUgBXWeeo4J987tbGI1uXOCzYll8ILiIBKIMCqACOSwGcIuNOEVdMEQPUmwKoAI5LAZxS3VPEa/KVHoVwqZ1QGQSoHwVQARyXAjjlqlHBA/d8bPN8nBFC6wUB6kEBVADHpQDOqHEuIskvmredcmjj2J97vG1lAGaQAqgAjksBnHHLvYikc0q51zRx/r3KIcB0UgAVwHEpgDWy3DLYKcvgSev2bfztdXc2p5itIQSYPgqgAjguBbCmOsvgOFcUd1IGAaaDAqgAjksBLED3esEscjk4OE4nzL9jw8nWEQKsBgVQARyXAliQalQwryi+4lv/3PidC7+x7I2n07B1hACsDAVQARyXAliwXmsGd4pWd+K6JzQuuu6uZZdDo4MAK0sBVADHpQDS1Dk6mGWtVzlcikGjg0kxBFg+BVABHJcCyFC9yuBy1xFWX7BVMTRSCLB0CqACOC4FkCVZ6XWEWQhPfdYBzf+OUgjQmwKoAI5LAWQskxwdTLkG8X+9+Mnbtq5RCgEWUwAVwHEpgEzMpEYH84u6808NKoVJMQRKowAqgONSAFkxvUYHu6d8u0cK875em1a3/9g2WQq/dNbzm0Vz/YU3GS0EiqIAKoDjUgDZITpHB1OvkcLmSN/xT268/ZKF71zSrxS+9/R1jTd+4roFj5lCBkqgACqA41IAWXWd5TCL2QXXfndoKcz7zz3tsCiA17fv2S6/KXT2xUGlMCmGwKxRABXAcSmATKVhpfCck5/e3D7mF952+YJiaAoZKIECqACOSwFkZnSXwjTKaOGkppCTYghMAwWwPgVwfeTkyCGRH0Y2Rs6M3BoZ5JTIH0cOjtwWeUvkbyOjUgCZecNGCycxhTxotDAphsCOpADWpwB+LvLXkWsjPxE5O3Jo5KmRByO9HBm5MvJ7kSx9J0X+KHJ05JrIKBRAamlYKVzqFHK/0ULFEFgNCmB9p4AfH/mnyHMjV+QdPVwQydZ2fPOsJYvk/ZFXNM+GUwApRncpTOOOFiqGwGpQAOP7avu2bpr/Xw3fa9/2kiOAl7UOt7k0clTrsKddItn0quwegSJk0Try4J9eULiykGVZ+8RrntO8fd2xBzeLWpa4lLc5WnjEgXs0S1ynfOzHURw7y1/KMvn12+/fVv5S3mbRzOKXpTNHHk8/75rmbZ5X8vGNt93XvAWgvzoWwPwx847IlyLfyDv62Dtyb+twmzzP+/vJtYZbOrI5AkXrLobdpTDP87HVLIZJOQTYro4F8D2RZ0RGmcbt+hGzaP16tw2RHF2ssl8E6DLKaOGOKobJqCHAQnUrgO+OvCzy/Miw0bl7It2jfXtFukcFOz0UybUCVR6IACPa0cUw1yxmsTOdDLBQ17fSmZWfR5a/vJL3eZFvR4bJi0ByDd8JzbOWSyLfj7gIBKZAlq9hF570uyI5H8tCmReLZMHr5gIUKJeLQOpTAN8XOT3y8kjn3n+5Tq/69f38yJ2RXMeX8mKPvEI49/77VCT/7FsjtoGBKTdqMczSlh/bqxy6MhnKpQDWpwB2fKte4JWRD7cOG1+MfCdyRvOs5ZcjWfqeGKk2gr4wMioFEKZIr2KYljJqqBhC/SmA9SmAq0UBhBkx6qjhjiqGaVA5VBxh5SiACuC4FECYcatRDPO/k/+NfuVw0GNJOYTxKIAK4LgUQKiplSqGebVz/p29/kyWw9TvsWHFMSmHMJwCGN8/2rcAdMjytBJb1mT5y4LWWfBSlsosnIMey3JXlb+Ut1lI8/6U5TDLo21tgGEUQIAlGLcY5v05OtevHA56bLnlcFAxTIPK4aDHgNnV9W2GJTIFDAyUxal7KjllCeueTq6mcvs9ln9Xv+nhLIdL3e8w/z3537IWkdKYAlYAx6UAAsvWrxympRbHfuXQWkRYTAE0BQywarIYdU8nV/o9liUrS1jndHPKj5uFtYgpP67XtHK/+4HJUwABZsxSymG/Ypj37+i1iKlfOVxuaQSWRwEEqJFe5XCpo4Z5/6DHJl0Ob7jj/hUZUUyDHoOSdX0Js0TWAAIzL8tRr/WGqd9jWcKWshYxi2eWw14XqvzuS57SeOvF32yfbTdsnWL+e/Lf0W8t4qDHUv5b899knWJ5rAE0AghQvCw/K70WMe/vN3L4rAMfN/Hp5kGPJaOKlE4BBGBZJlUOD9v/cROfbl5ucUyDyuGw4pgGFcRBj8GO1PVlxRKZAgZYoiw/vaaV+92fJWup082p32P9pqKHTTmnfo9V/978ty53Sjrl52RaeuWZAjYCCMAOlsWm18hhv/uXOqKY9w96bCVGFVOWt+VOSadho4v5sYNGD4c9Dp26vgRYIiOAAKssC0+vkcPU77FJjyrm353lq9/I4tb4v36PZekd9N/Nv3vY6OGwx/PvHzSyOOzxujECaAQQgBmXhaXXyGHq99ikRxXToJHFQY+lQaOLWc4GjR4Oe3zYyOKwx/PvMbJYPwogAEXK4raUcpgGPTaoIA56LC13WjqtZnlM+bH9CuKgx1hdXS83lsgUMAALZNlZ6pR0ynK11Gnp/DsGPZ7lcNDU86Bp6/w3Dvrvpvw395t6HvRYJf/t/aaeBz02LlPARgABYKKyrCx1Sjr1G13Mjx00ejjo8WFTz+OMPGZB6zd6OOixyqDRxUGPMRkKIABMiSxtS52WTqtRHgcVxJUuj4xPAQSAGZClrN/oYer3+EqVx0EFcSXLI5OhAAJAzfUrh5V+jw8qj/mx/QrioMfSOOWRyeh6ilkiF4EAULScms3RuSxo3QVy0GO5rq/XRS/DHpsEF4EogONSAAFgmZZbHselACqA41IAAWDGKIDWAAIAFEcBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDC1KkAHhv5dOSuSL6JzImRQZ4XyY/rziERAIDaqlMB3DVyQ+QNzbPRPTmyT0e+HQEAqK06FcBLIr8bubB5Nrp/itzTkUciAAC1ZQ1go3Fd5O7IFyLPzzsG2CWSb/lRZfcIAMBMKbkAZul7beSUyMmRWyNZAnMtYT/rI1s6sjkCADBT6vpewHkxx0mRi5pno8uLSPLPvqx5tliOAGYqOQK42XsBA8Ds8F7ApoC7fTnypNZhTw9F8oVS5YEIAMBMUQAXWhfJqWEAgNqqUwHcLXJ4O+mgSB4f0DxrNDZEzm8dNr05knsF5ojf0yL5eK4HfE8EAKC26lQAj4jkFb2Z9I5IHv9R86y1x19VBtPOkT+N3Bi5MnJ05CWRpW4jAwAwU+p6EciOkld+bHERCADMDheBWAMIAFAcBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEAClOnAnhs5NORuyJbIydGhnlu5GuRH0X+X+Q3IgAAtVanArhr5IbIG5pnwx0U+Wzkysi6yDmRP4ucEgEAqK017du6yRHAkyIXNc96e3vkZZGnNM9a3h85LHJk82y4tZEtobF2bR4CANNufn6+MTc3l4f5P/N5UJqS1wBmybusdbjNpZEjIj/ZPAMAqKGSC+DekXtbh9vk+U9E9myeLbZLJIf6quweAQCYKaVfBZxTxZ2qKfHu+yvrI1s6sjkCADBTSi6A90RyFLDTXpF/j/xL82yxDZFcL1BlvwgAwEwpuQBeHTmudbjNCyNfjfxb82yxhyK5WLTKAxEAgJlSpwK4W+TwdlJu85LHBzTPWqN357cOm/KK35+NvCOSVwK/KvLqyJ9GAABqq04FMK/eva6dlMUuj/+oedZo7BOpymDaFDkh8rzI9ZHfi/x25G8iAAC1Vdd9AHcU+wACwIyxD6CrgAEAiqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMLUrQD+ZmRT5EeRr0WOifRzRmRrjzw6AgBQW3UqgKdGzo2cHVkXuTJySeSASD/zkX26kuURAKC26lQA/1vkg5G/jHwz8ubIHZHXR/rJEb97ugIAUGt1KYA7R54Zuax5tl2eH9U67Gm3yO2RzZHPRHLkEACg1upSAPeM7BS5t3m2XZ7v3Tpc5JZIrgN8WeQVkZz6vSrypEg/u0TWdmT3CADATKnbRSA5pdtpTaT7vsqXIx+N3BDJ9YK/EvlW5I2RftZHtnQkRw4BAGZKXQrgfZFHIt2jfXtFukcF+/lx5NrIoBHADZG5juwXAQCYKXUpgA9HctuX45pn2+X5xtbhUDlaeHjk7uZZbw9F8srhKg9EAABmSp2mgN8R+fXIqyJPibwzklvAvD+Szo/kCF7lDyIvijwxksUvryDO2+rjAQBqqU4F8IJIbv3y+5HrI8dGTojkVb4py2Du81f5qcgHIrllTF4tvG8k/8xXIgAAtZXTnixfXgm8JTTWrs1DAGDazc/PN+bmcil/cz1/LukqTt2uAgYAYAgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABSmbgXwNyObIj+KfC1yTGSQUyI3Rx5q354UAQCotToVwFMj50bOjqyLXBm5JHJApJcjIxdEPhI5rH37ycizI6vuwLMu3hYAgEla076tg2siX4+8vnnW8s3IRZH1zbOFsvytjRzfPGv5XOT+yCuaZ8Pln98SGmvX5uFk9Cp933nbS9pHAMA45ufnG3Nzc3mY/zOfB6WpywjgzpFnRi5rnm2X50e1DhfJEcDuj7800u/j0y6RbHpVdo9MVL8RPyOBAMCk1KUA7hnZKXJv82y7PN+7dbhI3r+Uj085krilI5sjAAAzpW4XgWxt31Zyirv7vk5L/fgNkRwurrJfBABgptSlAN4XeSTSPXq3V6R7lK9yT2QpH5/yauFcK1DlgchE9VvrZw0gADApdSmAD0dy25fjmmfb5fnG1uEiV0e6P/6FkX4fv8N0lz3lDwCYpDpdBZzbwORWLr8RyXL32shrIk+L3B45P3JnpLoiOC/2uCLylsinIi+PvDVydCSvKB7FilwFDACsHFcB12sNYG7r8ubI70eujxwbOSGS5S/lfoD7tA6bcqTvtMgrIzdGzohkiRy1/AEAzKQ6jQCuBiOAADBjjADW7ypgAACGUAABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIXxTiDjab4TyB133OGdQABgRuQ7gey///55WOw7gSiA49k3srl1CADMmP0id7YOy6IAjiefvydEHmieTdbukSyX+eJcib9/1ng+FvJ8LOT5WMxzspDnYyHPR+s5uCuytXkGUyLnlPNFaW65xfOxkOdjIc/HYp6ThTwfC3k+CuciEACAwiiAAACF2al9y3R6JPLF9i2ej26ej4U8H4t5ThbyfCzk+QAAAAAAAAAAAAAAAABYqt+MbIr8KPK1yDGREvxhJDcm7cw9kUq+80p+TO7c/sNIXr32tEhdHBv5dKTamf7ESKdRPv/HRT4S2dJOHv9UZFYNe04+HOl8vWS+HOm0S+TdkfsiD0b+TyTf/WAWrY9cG8l3bvinyEWRJ0c6jfL5HhDJ5zUfz4/7s8jOkVkzyvORXyfdr5G/jnSqy9fN6yM3RvK9bTNXR46PVEp6bcDMOTXycOTXI0+JnBv510h+UdZdlptvRPbuyOMjlTMj+U3t5MjTI/lNPItBvp1PHeQ36rdG8vPLH1LdZWeUz/+SyE2RI9vJ4/xmPquGPSdZAPNz7nzN7BHp9OeRfMurX4qsi1weuT4yi9tgfS5yRiSL/2GRz0Ruj+waqQz7fPM2Xxd5fz6eH5fvhZrFYNaM8nxkAfxApPM1MhfpVJevm/8UOSHyc+2cHcmfJ9UviiW9NmDmXBPJL9JO34xsaB3WWhbA/GbUS45+3R3JElTJ32a/H3ld86xeusvOKJ9//sKQf+7ZzbOW50Tyvu5RkVnUrwDmqE8/+YM+fwDmL1aVJ0Ry37MXNc9mW/6ClM9LjpSmUT7fLNV5nvdXTovkjMOsvy1Y9/ORsgDmL9L91P3r5nuRV0dKf23QxTuBTJccZn9m5LLm2XZ5flTrsPaeFMlRrZwCzxGuJ0bSQZH8zb3zuXko8veREp6bUT7/HLnI6av8JaKS06F5X52fo+dFcvrvW5HzIntFKvn19JORzuctX1850lyH56Qaycof8mmUzzdfJ3me91cujeQvFPnnZ1n381H51UhOZ/5D5E8jnaPmdf26ydG8LG85GppTwaW/NuiiAE6XPSP5RXtv82y7PM8f/nWX34D/SyR/G31NJD/njZGfbh+nUp+bUT7/vM0i1C3vq+tzlFN3+cP9BZH/HnlWJKev8gdWys87Rz3ub55tV4fXTY4KvyPypUj+0E6jfL552/06yo/PPzfLz0mv5yN9LPKKSP6i8MeRUyIXRir5Odfp6+bQSC4byl8Q3x85KXJzJD+XUl8b9KAATqeceuiU39i676uj/GH+N5Fcg/J3kZdE0q+1b1Opz01l2Off67mo83N0QeTiSP7AzzVbOYWVa5+q104/dXhO3hN5RiTLzTAlvE76PR85KpzfT/I1krMKvxzJtW3/MVKp0/Nxa+TwSE5j53Kiv4o8NdJPCa8NelAAp0tOUeT6i+7ftHJKq/u3shLkVWhZBnNauLoauNTnZpTPPz/mZ1qHC+S6qFJeP7lOMi8CyNdMyuckl1bkVZ6dZv11k4vyXxZ5fiQX9VdG+XzzY7pfR/nxOT04q89Jv+ejl69H/i3S+Rqp09dNjtb9Y+SrkbxK+obImyKlvjboQwGcLvmFm9u+HNc82y7Pcyq0NDmNlwu084d6rgnMb06dz01+M3tupITnZpTPP9f55Bqon2+eteTC9ryvlNdPLhfYP5KvmZRfT/nDvvN52yeSV1HP4nOSIzE50pVXRee0d74uOo3y+ebrJM/z/soLIzllmH9+lgx7PnrJK2Kz0FSvkbp/3eRzlN9LS3ttwMyptoF5VSTLzzsjuZ7jZyN1l4uzs9DkBQ/5DTin9HLbk+pzzytg86rXXNOS36Q+HsnFynXZBma3SE7dZHK65b+2j6stgEb5/HMaPX/jz+mfTO4JNsvbwAx6TvKxfM3kwvUDI7nGK3+Q5QhQ53OS02B3RH4xkltbfCEyq9vAvC+Sr4H8OsmRmiqPiVSGfb7VVh85LZqP58flx8/iVh/Dno+DI78fOSKSr5HcIiV3VchRwM7//9fl6+acSO4bm59rrgXMbWByVqkqfSW9NmAm5UbQ34lUv3V1bmlQZ9W+dlmAc++pXA/YuXYlf5PNrWLyN/fcliCvgM0iVBdZYLLkdCe3OkmjfP65B95HI9VGsHk8yxtBD3pO8od8XqGYi/XzNZNTv3l/jgB2enQkf4D9S+QHkfzB3v0xs6LXc5HJvfAqo3y+WaBzz7x8PD8uP766cGaW9HouMtXzkZ93fp3k55jfT3Nq9F2R7r0i6/J188FI9bMjvy6yyHWO+JX02gAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACAGmg0/j/Edp7X3539/wAAAABJRU5ErkJggg==\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots()\n",
    "ax.plot(curve_data[\"Temperature (K)\"], curve_data['log Ohm'], '.')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'CX-1050-SD-HT-1'"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "curve.curve_name()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "curve_x = lake.curve[23]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "curve_x_data = curve_x.get_data()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "dict_keys(['Temperature (K)', 'log Ohm'])"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "curve_x_data.keys()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": "/* Put everything inside the global mpl namespace */\nwindow.mpl = {};\n\n\nmpl.get_websocket_type = function() {\n    if (typeof(WebSocket) !== 'undefined') {\n        return WebSocket;\n    } else if (typeof(MozWebSocket) !== 'undefined') {\n        return MozWebSocket;\n    } else {\n        alert('Your browser does not have WebSocket support.' +\n              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n              'Firefox 4 and 5 are also supported but you ' +\n              'have to enable WebSockets in about:config.');\n    };\n}\n\nmpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n    this.id = figure_id;\n\n    this.ws = websocket;\n\n    this.supports_binary = (this.ws.binaryType != undefined);\n\n    if (!this.supports_binary) {\n        var warnings = document.getElementById(\"mpl-warnings\");\n        if (warnings) {\n            warnings.style.display = 'block';\n            warnings.textContent = (\n                \"This browser does not support binary websocket messages. \" +\n                    \"Performance may be slow.\");\n        }\n    }\n\n    this.imageObj = new Image();\n\n    this.context = undefined;\n    this.message = undefined;\n    this.canvas = undefined;\n    this.rubberband_canvas = undefined;\n    this.rubberband_context = undefined;\n    this.format_dropdown = undefined;\n\n    this.image_mode = 'full';\n\n    this.root = $('<div/>');\n    this._root_extra_style(this.root)\n    this.root.attr('style', 'display: inline-block');\n\n    $(parent_element).append(this.root);\n\n    this._init_header(this);\n    this._init_canvas(this);\n    this._init_toolbar(this);\n\n    var fig = this;\n\n    this.waiting = false;\n\n    this.ws.onopen =  function () {\n            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n            fig.send_message(\"send_image_mode\", {});\n            if (mpl.ratio != 1) {\n                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n            }\n            fig.send_message(\"refresh\", {});\n        }\n\n    this.imageObj.onload = function() {\n            if (fig.image_mode == 'full') {\n                // Full images could contain transparency (where diff images\n                // almost always do), so we need to clear the canvas so that\n                // there is no ghosting.\n                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n            }\n            fig.context.drawImage(fig.imageObj, 0, 0);\n        };\n\n    this.imageObj.onunload = function() {\n        this.ws.close();\n    }\n\n    this.ws.onmessage = this._make_on_message_function(this);\n\n    this.ondownload = ondownload;\n}\n\nmpl.figure.prototype._init_header = function() {\n    var titlebar = $(\n        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n        'ui-helper-clearfix\"/>');\n    var titletext = $(\n        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n        'text-align: center; padding: 3px;\"/>');\n    titlebar.append(titletext)\n    this.root.append(titlebar);\n    this.header = titletext[0];\n}\n\n\n\nmpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n\n}\n\n\nmpl.figure.prototype._root_extra_style = function(canvas_div) {\n\n}\n\nmpl.figure.prototype._init_canvas = function() {\n    var fig = this;\n\n    var canvas_div = $('<div/>');\n\n    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n\n    function canvas_keyboard_event(event) {\n        return fig.key_event(event, event['data']);\n    }\n\n    canvas_div.keydown('key_press', canvas_keyboard_event);\n    canvas_div.keyup('key_release', canvas_keyboard_event);\n    this.canvas_div = canvas_div\n    this._canvas_extra_style(canvas_div)\n    this.root.append(canvas_div);\n\n    var canvas = $('<canvas/>');\n    canvas.addClass('mpl-canvas');\n    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n\n    this.canvas = canvas[0];\n    this.context = canvas[0].getContext(\"2d\");\n\n    var backingStore = this.context.backingStorePixelRatio ||\n\tthis.context.webkitBackingStorePixelRatio ||\n\tthis.context.mozBackingStorePixelRatio ||\n\tthis.context.msBackingStorePixelRatio ||\n\tthis.context.oBackingStorePixelRatio ||\n\tthis.context.backingStorePixelRatio || 1;\n\n    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n\n    var rubberband = $('<canvas/>');\n    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n\n    var pass_mouse_events = true;\n\n    canvas_div.resizable({\n        start: function(event, ui) {\n            pass_mouse_events = false;\n        },\n        resize: function(event, ui) {\n            fig.request_resize(ui.size.width, ui.size.height);\n        },\n        stop: function(event, ui) {\n            pass_mouse_events = true;\n            fig.request_resize(ui.size.width, ui.size.height);\n        },\n    });\n\n    function mouse_event_fn(event) {\n        if (pass_mouse_events)\n            return fig.mouse_event(event, event['data']);\n    }\n\n    rubberband.mousedown('button_press', mouse_event_fn);\n    rubberband.mouseup('button_release', mouse_event_fn);\n    // Throttle sequential mouse events to 1 every 20ms.\n    rubberband.mousemove('motion_notify', mouse_event_fn);\n\n    rubberband.mouseenter('figure_enter', mouse_event_fn);\n    rubberband.mouseleave('figure_leave', mouse_event_fn);\n\n    canvas_div.on(\"wheel\", function (event) {\n        event = event.originalEvent;\n        event['data'] = 'scroll'\n        if (event.deltaY < 0) {\n            event.step = 1;\n        } else {\n            event.step = -1;\n        }\n        mouse_event_fn(event);\n    });\n\n    canvas_div.append(canvas);\n    canvas_div.append(rubberband);\n\n    this.rubberband = rubberband;\n    this.rubberband_canvas = rubberband[0];\n    this.rubberband_context = rubberband[0].getContext(\"2d\");\n    this.rubberband_context.strokeStyle = \"#000000\";\n\n    this._resize_canvas = function(width, height) {\n        // Keep the size of the canvas, canvas container, and rubber band\n        // canvas in synch.\n        canvas_div.css('width', width)\n        canvas_div.css('height', height)\n\n        canvas.attr('width', width * mpl.ratio);\n        canvas.attr('height', height * mpl.ratio);\n        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n\n        rubberband.attr('width', width);\n        rubberband.attr('height', height);\n    }\n\n    // Set the figure to an initial 600x600px, this will subsequently be updated\n    // upon first draw.\n    this._resize_canvas(600, 600);\n\n    // Disable right mouse context menu.\n    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n        return false;\n    });\n\n    function set_focus () {\n        canvas.focus();\n        canvas_div.focus();\n    }\n\n    window.setTimeout(set_focus, 100);\n}\n\nmpl.figure.prototype._init_toolbar = function() {\n    var fig = this;\n\n    var nav_element = $('<div/>')\n    nav_element.attr('style', 'width: 100%');\n    this.root.append(nav_element);\n\n    // Define a callback function for later on.\n    function toolbar_event(event) {\n        return fig.toolbar_button_onclick(event['data']);\n    }\n    function toolbar_mouse_event(event) {\n        return fig.toolbar_button_onmouseover(event['data']);\n    }\n\n    for(var toolbar_ind in mpl.toolbar_items) {\n        var name = mpl.toolbar_items[toolbar_ind][0];\n        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n        var image = mpl.toolbar_items[toolbar_ind][2];\n        var method_name = mpl.toolbar_items[toolbar_ind][3];\n\n        if (!name) {\n            // put a spacer in here.\n            continue;\n        }\n        var button = $('<button/>');\n        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n                        'ui-button-icon-only');\n        button.attr('role', 'button');\n        button.attr('aria-disabled', 'false');\n        button.click(method_name, toolbar_event);\n        button.mouseover(tooltip, toolbar_mouse_event);\n\n        var icon_img = $('<span/>');\n        icon_img.addClass('ui-button-icon-primary ui-icon');\n        icon_img.addClass(image);\n        icon_img.addClass('ui-corner-all');\n\n        var tooltip_span = $('<span/>');\n        tooltip_span.addClass('ui-button-text');\n        tooltip_span.html(tooltip);\n\n        button.append(icon_img);\n        button.append(tooltip_span);\n\n        nav_element.append(button);\n    }\n\n    var fmt_picker_span = $('<span/>');\n\n    var fmt_picker = $('<select/>');\n    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n    fmt_picker_span.append(fmt_picker);\n    nav_element.append(fmt_picker_span);\n    this.format_dropdown = fmt_picker[0];\n\n    for (var ind in mpl.extensions) {\n        var fmt = mpl.extensions[ind];\n        var option = $(\n            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n        fmt_picker.append(option)\n    }\n\n    // Add hover states to the ui-buttons\n    $( \".ui-button\" ).hover(\n        function() { $(this).addClass(\"ui-state-hover\");},\n        function() { $(this).removeClass(\"ui-state-hover\");}\n    );\n\n    var status_bar = $('<span class=\"mpl-message\"/>');\n    nav_element.append(status_bar);\n    this.message = status_bar[0];\n}\n\nmpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n    // which will in turn request a refresh of the image.\n    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n}\n\nmpl.figure.prototype.send_message = function(type, properties) {\n    properties['type'] = type;\n    properties['figure_id'] = this.id;\n    this.ws.send(JSON.stringify(properties));\n}\n\nmpl.figure.prototype.send_draw_message = function() {\n    if (!this.waiting) {\n        this.waiting = true;\n        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n    }\n}\n\n\nmpl.figure.prototype.handle_save = function(fig, msg) {\n    var format_dropdown = fig.format_dropdown;\n    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n    fig.ondownload(fig, format);\n}\n\n\nmpl.figure.prototype.handle_resize = function(fig, msg) {\n    var size = msg['size'];\n    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n        fig._resize_canvas(size[0], size[1]);\n        fig.send_message(\"refresh\", {});\n    };\n}\n\nmpl.figure.prototype.handle_rubberband = function(fig, msg) {\n    var x0 = msg['x0'] / mpl.ratio;\n    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n    var x1 = msg['x1'] / mpl.ratio;\n    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n    x0 = Math.floor(x0) + 0.5;\n    y0 = Math.floor(y0) + 0.5;\n    x1 = Math.floor(x1) + 0.5;\n    y1 = Math.floor(y1) + 0.5;\n    var min_x = Math.min(x0, x1);\n    var min_y = Math.min(y0, y1);\n    var width = Math.abs(x1 - x0);\n    var height = Math.abs(y1 - y0);\n\n    fig.rubberband_context.clearRect(\n        0, 0, fig.canvas.width, fig.canvas.height);\n\n    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n}\n\nmpl.figure.prototype.handle_figure_label = function(fig, msg) {\n    // Updates the figure title.\n    fig.header.textContent = msg['label'];\n}\n\nmpl.figure.prototype.handle_cursor = function(fig, msg) {\n    var cursor = msg['cursor'];\n    switch(cursor)\n    {\n    case 0:\n        cursor = 'pointer';\n        break;\n    case 1:\n        cursor = 'default';\n        break;\n    case 2:\n        cursor = 'crosshair';\n        break;\n    case 3:\n        cursor = 'move';\n        break;\n    }\n    fig.rubberband_canvas.style.cursor = cursor;\n}\n\nmpl.figure.prototype.handle_message = function(fig, msg) {\n    fig.message.textContent = msg['message'];\n}\n\nmpl.figure.prototype.handle_draw = function(fig, msg) {\n    // Request the server to send over a new figure.\n    fig.send_draw_message();\n}\n\nmpl.figure.prototype.handle_image_mode = function(fig, msg) {\n    fig.image_mode = msg['mode'];\n}\n\nmpl.figure.prototype.updated_canvas_event = function() {\n    // Called whenever the canvas gets updated.\n    this.send_message(\"ack\", {});\n}\n\n// A function to construct a web socket function for onmessage handling.\n// Called in the figure constructor.\nmpl.figure.prototype._make_on_message_function = function(fig) {\n    return function socket_on_message(evt) {\n        if (evt.data instanceof Blob) {\n            /* FIXME: We get \"Resource interpreted as Image but\n             * transferred with MIME type text/plain:\" errors on\n             * Chrome.  But how to set the MIME type?  It doesn't seem\n             * to be part of the websocket stream */\n            evt.data.type = \"image/png\";\n\n            /* Free the memory for the previous frames */\n            if (fig.imageObj.src) {\n                (window.URL || window.webkitURL).revokeObjectURL(\n                    fig.imageObj.src);\n            }\n\n            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n                evt.data);\n            fig.updated_canvas_event();\n            fig.waiting = false;\n            return;\n        }\n        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n            fig.imageObj.src = evt.data;\n            fig.updated_canvas_event();\n            fig.waiting = false;\n            return;\n        }\n\n        var msg = JSON.parse(evt.data);\n        var msg_type = msg['type'];\n\n        // Call the  \"handle_{type}\" callback, which takes\n        // the figure and JSON message as its only arguments.\n        try {\n            var callback = fig[\"handle_\" + msg_type];\n        } catch (e) {\n            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n            return;\n        }\n\n        if (callback) {\n            try {\n                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n                callback(fig, msg);\n            } catch (e) {\n                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n            }\n        }\n    };\n}\n\n// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\nmpl.findpos = function(e) {\n    //this section is from http://www.quirksmode.org/js/events_properties.html\n    var targ;\n    if (!e)\n        e = window.event;\n    if (e.target)\n        targ = e.target;\n    else if (e.srcElement)\n        targ = e.srcElement;\n    if (targ.nodeType == 3) // defeat Safari bug\n        targ = targ.parentNode;\n\n    // jQuery normalizes the pageX and pageY\n    // pageX,Y are the mouse positions relative to the document\n    // offset() returns the position of the element relative to the document\n    var x = e.pageX - $(targ).offset().left;\n    var y = e.pageY - $(targ).offset().top;\n\n    return {\"x\": x, \"y\": y};\n};\n\n/*\n * return a copy of an object with only non-object keys\n * we need this to avoid circular references\n * http://stackoverflow.com/a/24161582/3208463\n */\nfunction simpleKeys (original) {\n  return Object.keys(original).reduce(function (obj, key) {\n    if (typeof original[key] !== 'object')\n        obj[key] = original[key]\n    return obj;\n  }, {});\n}\n\nmpl.figure.prototype.mouse_event = function(event, name) {\n    var canvas_pos = mpl.findpos(event)\n\n    if (name === 'button_press')\n    {\n        this.canvas.focus();\n        this.canvas_div.focus();\n    }\n\n    var x = canvas_pos.x * mpl.ratio;\n    var y = canvas_pos.y * mpl.ratio;\n\n    this.send_message(name, {x: x, y: y, button: event.button,\n                             step: event.step,\n                             guiEvent: simpleKeys(event)});\n\n    /* This prevents the web browser from automatically changing to\n     * the text insertion cursor when the button is pressed.  We want\n     * to control all of the cursor setting manually through the\n     * 'cursor' event from matplotlib */\n    event.preventDefault();\n    return false;\n}\n\nmpl.figure.prototype._key_event_extra = function(event, name) {\n    // Handle any extra behaviour associated with a key event\n}\n\nmpl.figure.prototype.key_event = function(event, name) {\n\n    // Prevent repeat events\n    if (name == 'key_press')\n    {\n        if (event.which === this._key)\n            return;\n        else\n            this._key = event.which;\n    }\n    if (name == 'key_release')\n        this._key = null;\n\n    var value = '';\n    if (event.ctrlKey && event.which != 17)\n        value += \"ctrl+\";\n    if (event.altKey && event.which != 18)\n        value += \"alt+\";\n    if (event.shiftKey && event.which != 16)\n        value += \"shift+\";\n\n    value += 'k';\n    value += event.which.toString();\n\n    this._key_event_extra(event, name);\n\n    this.send_message(name, {key: value,\n                             guiEvent: simpleKeys(event)});\n    return false;\n}\n\nmpl.figure.prototype.toolbar_button_onclick = function(name) {\n    if (name == 'download') {\n        this.handle_save(this, null);\n    } else {\n        this.send_message(\"toolbar_button\", {name: name});\n    }\n};\n\nmpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n    this.message.textContent = tooltip;\n};\nmpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n\nmpl.extensions = [\"eps\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\"];\n\nmpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n    // Create a \"websocket\"-like object which calls the given IPython comm\n    // object with the appropriate methods. Currently this is a non binary\n    // socket, so there is still some room for performance tuning.\n    var ws = {};\n\n    ws.close = function() {\n        comm.close()\n    };\n    ws.send = function(m) {\n        //console.log('sending', m);\n        comm.send(m);\n    };\n    // Register the callback with on_msg.\n    comm.on_msg(function(msg) {\n        //console.log('receiving', msg['content']['data'], msg);\n        // Pass the mpl event to the overriden (by mpl) onmessage function.\n        ws.onmessage(msg['content']['data'])\n    });\n    return ws;\n}\n\nmpl.mpl_figure_comm = function(comm, msg) {\n    // This is the function which gets called when the mpl process\n    // starts-up an IPython Comm through the \"matplotlib\" channel.\n\n    var id = msg.content.data.id;\n    // Get hold of the div created by the display call when the Comm\n    // socket was opened in Python.\n    var element = $(\"#\" + id);\n    var ws_proxy = comm_websocket_adapter(comm)\n\n    function ondownload(figure, format) {\n        window.open(figure.imageObj.src);\n    }\n\n    var fig = new mpl.figure(id, ws_proxy,\n                           ondownload,\n                           element.get(0));\n\n    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n    // web socket which is closed, not our websocket->open comm proxy.\n    ws_proxy.onopen();\n\n    fig.parent_element = element.get(0);\n    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n    if (!fig.cell_info) {\n        console.error(\"Failed to find cell for figure\", id, fig);\n        return;\n    }\n\n    var output_index = fig.cell_info[2]\n    var cell = fig.cell_info[0];\n\n};\n\nmpl.figure.prototype.handle_close = function(fig, msg) {\n    var width = fig.canvas.width/mpl.ratio\n    fig.root.unbind('remove')\n\n    // Update the output cell to use the data from the current canvas.\n    fig.push_to_output();\n    var dataURL = fig.canvas.toDataURL();\n    // Re-enable the keyboard manager in IPython - without this line, in FF,\n    // the notebook keyboard shortcuts fail.\n    IPython.keyboard_manager.enable()\n    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n    fig.close_ws(fig, msg);\n}\n\nmpl.figure.prototype.close_ws = function(fig, msg){\n    fig.send_message('closing', msg);\n    // fig.ws.close()\n}\n\nmpl.figure.prototype.push_to_output = function(remove_interactive) {\n    // Turn the data on the canvas into data in the output cell.\n    var width = this.canvas.width/mpl.ratio\n    var dataURL = this.canvas.toDataURL();\n    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n}\n\nmpl.figure.prototype.updated_canvas_event = function() {\n    // Tell IPython that the notebook contents must change.\n    IPython.notebook.set_dirty(true);\n    this.send_message(\"ack\", {});\n    var fig = this;\n    // Wait a second, then push the new image to the DOM so\n    // that it is saved nicely (might be nice to debounce this).\n    setTimeout(function () { fig.push_to_output() }, 1000);\n}\n\nmpl.figure.prototype._init_toolbar = function() {\n    var fig = this;\n\n    var nav_element = $('<div/>')\n    nav_element.attr('style', 'width: 100%');\n    this.root.append(nav_element);\n\n    // Define a callback function for later on.\n    function toolbar_event(event) {\n        return fig.toolbar_button_onclick(event['data']);\n    }\n    function toolbar_mouse_event(event) {\n        return fig.toolbar_button_onmouseover(event['data']);\n    }\n\n    for(var toolbar_ind in mpl.toolbar_items){\n        var name = mpl.toolbar_items[toolbar_ind][0];\n        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n        var image = mpl.toolbar_items[toolbar_ind][2];\n        var method_name = mpl.toolbar_items[toolbar_ind][3];\n\n        if (!name) { continue; };\n\n        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n        button.click(method_name, toolbar_event);\n        button.mouseover(tooltip, toolbar_mouse_event);\n        nav_element.append(button);\n    }\n\n    // Add the status bar.\n    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n    nav_element.append(status_bar);\n    this.message = status_bar[0];\n\n    // Add the close button to the window.\n    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n    button.click(function (evt) { fig.handle_close(fig, {}); } );\n    button.mouseover('Stop Interaction', toolbar_mouse_event);\n    buttongrp.append(button);\n    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n    titlebar.prepend(buttongrp);\n}\n\nmpl.figure.prototype._root_extra_style = function(el){\n    var fig = this\n    el.on(\"remove\", function(){\n\tfig.close_ws(fig, {});\n    });\n}\n\nmpl.figure.prototype._canvas_extra_style = function(el){\n    // this is important to make the div 'focusable\n    el.attr('tabindex', 0)\n    // reach out to IPython and tell the keyboard manager to turn it's self\n    // off when our div gets focus\n\n    // location in version 3\n    if (IPython.notebook.keyboard_manager) {\n        IPython.notebook.keyboard_manager.register_events(el);\n    }\n    else {\n        // location in version 2\n        IPython.keyboard_manager.register_events(el);\n    }\n\n}\n\nmpl.figure.prototype._key_event_extra = function(event, name) {\n    var manager = IPython.notebook.keyboard_manager;\n    if (!manager)\n        manager = IPython.keyboard_manager;\n\n    // Check for shift+enter\n    if (event.shiftKey && event.which == 13) {\n        this.canvas_div.blur();\n        // select the cell after this one\n        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n        IPython.notebook.select(index + 1);\n    }\n}\n\nmpl.figure.prototype.handle_save = function(fig, msg) {\n    fig.ondownload(fig, null);\n}\n\n\nmpl.find_output_cell = function(html_output) {\n    // Return the cell and output element which can be found *uniquely* in the notebook.\n    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n    // IPython event is triggered only after the cells have been serialised, which for\n    // our purposes (turning an active figure into a static one), is too late.\n    var cells = IPython.notebook.get_cells();\n    var ncells = cells.length;\n    for (var i=0; i<ncells; i++) {\n        var cell = cells[i];\n        if (cell.cell_type === 'code'){\n            for (var j=0; j<cell.output_area.outputs.length; j++) {\n                var data = cell.output_area.outputs[j];\n                if (data.data) {\n                    // IPython >= 3 moved mimebundle to data attribute of output\n                    data = data.data;\n                }\n                if (data['text/html'] == html_output) {\n                    return [cell, data, j];\n                }\n            }\n        }\n    }\n}\n\n// Register the function which deals with the matplotlib target/channel.\n// The kernel may be null if the page has been refreshed.\nif (IPython.notebook.kernel != null) {\n    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n}\n",
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAB/WSURBVHhe7d0NlJ11fSfwO4bykpQZXYIlFkIinuILwW5dLCgUoeKqeKykPSeAnlbW5XRxe2q7W1eDa9euSui+oBbUttajS6vC7mlKa6mWti5CjPXErivkaNGNyRIwSmfBDAYEG2Z/v2fuk1yGuTN35s5M7n3+n88539znuXMD+sww853/y3NbAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAANMtI+5GFyev3rMjD1RkAMCyOj3w7MlmdFUYB7M+PR+6bOgQAhszJkfunDsuiAPZnNLJ/7969rdHRPAQABt3ExETrlFNOycOxPM2D0iiA/akKYFAAAWBIZAEcG8vuV24BfFr7EQCAQiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCNKUAXhW5K5J38858MfKqyGx+PvK1yGPtx0siAACN15QCeF/k7ZF/1s7nIn8aeUFkJudEbo78YeSF7cf/HvnpyEDYt//R1vZd49UjAMBiavJ7AT8YeWvko9XZk2X5yzfv7Rwl/Gzkochl1VlvluS9gG/ecW9r89a7W09MRkOPz9CWjRtam85a2/4oANAP7wXczDWAKyKXRlZFcip4JjkCeNvU4SF/GXnJ1GFXx0Sy6dU5PrKocsSvLn8pH6/eutNIIACwaJpUADdEvh/JNX2/G8k1fbm2byYnRb47dXhInufzs9kc2d+RnHpeVLvHDxwqf7WDk5OtPeOPtM8AAPrTpAJ4T+QnI2dHPhz5b5HnR7qZVrOq6fDpz023JZLDxXVOjiyq9atXVdO+nVaMjLTWrV7ZPgMA6E+TCuDjkf8T+XIkR+q+GnlLZCbfiUwf7XtmZPqo4HQ5uljvNM48HFlUa8aOq9b8ZelL+XjNxjOq5wEAFkOTN4H8TWRv5I3V2ZPlJpBcv/fq6mzKZyLfixzxTSAp1/zltG+O/Cl/ALB4bAJpzgjgNZHzIusiuRbwvZGXRT4RSTdGcvq29oHIKyJvizy3/fjyyPsjAyFL3zmnnaD8AQCLrikF8McieS+/XAeYI395P79XRv4qkvIeKmumDivbI7lT+IpI3kA6Rwk3Rb4UAQBotCZPAS+HJZsCBgCWhingZm0CAQCgBwogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwTSmAmyM7Ig9HHojcEjk9Mpdfi9wTeTSyN/K+yLERAIDGakoBPD/ywcjZkYsiR0Vui6yKdPP6yLWR34o8L/KmyKbIlggAQGONtB+b5sRIjgRmMbwjn5jBDZEsfj9bnU35r5EXR86rzuY2GtkfWqOjeQgADLqJiYnW2NhYHuYfE3lQmqauAaw+q+HB9uNMtkVeFMnCl54deXXk1upsZsdEsunVOT4CADBUmlgAc1TzukgWvJ35RBc3Rd4Zydf9MLIr8j8jOS3cTa413N+R+yIAAEOliQUwp3bPjFxWnXX3ssg7Im+O/FRkY+Q1kSyF3eT6wBxdrHNyBABgqDRtDeD1kddFfiayO5+YxZ2Rv428tTqb8obI70d+NPJEPjEHawABYMhYA9icEcAssjnyl6N4F0bmKn9pZWR6yTsYyX9WUzfHAAA0pgDmLWBy9O7ySN4L8KR2jovUbox03uLl05GrIpdG1kfy9jHvjvxZJIsgAEAjNaUAZpHLYdzbI/s6kvf1q62NrJk6rLwnkrd9ycevRT4a+cvIL0cAABrLVGd/rAEEgCFjDWAzdwEDADALBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogANu3/5HW9t3jVePAACLYaT9yMKMRvaH1uhoHi6um3fc29q89e7WE5PR1OMztWXjhtams9a2PwoALMTExERrbGwsD/OPiTwojRHAAZUjfnX5S/l49dadRgIBgL4pgANq9/iBQ+WvdnBysrVn/JH2GQDAwiiAA2r96lXVtG+nFSMjrXWrV7bPAAAWRgEcUGvGjqvW/GXpS/l4zcYzqucBAPphE0h/lnQTSMo1fzntmyN/yh8A9M8mECOAAy9L3zmnnaD8AQCLRgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwjSlAG6O7Ig8HHkgckvk9Mhcnh75YGRf5AeRr0deHQEAaKymFMDzI1nkzo5cFDkqcltkVaSboyN/FVkX+YVIFsYrI/dHAAAaa6T92DQnRnIkMIvhHfnEDP5V5K2R50Z+mE8swGhkf2iNjuYhADDoJiYmWmNjY3mYf0zkQWmaugaw+qyGB9uPM3lt5IuRHDn8bmRn5OrIikg3x0Sy6dU5PgIAMFSaWABzVPO6yLZIlrpunh3Jqd8sfLnu7z2Rfxt5R6SbXGu4vyP3RQAAhkoTp4BzRO/iyLmR2QraNyLHRtZHDuYT4d9Eclp4TXX2VDkCmKnlCOB9poABYHiYAm7eCOD1kZzavSAy1+hc7vzNEliXv5S7gE+K5AaRmTwWyS+UOrnrGABgqDSlAOZI5g2RjZELI7sjc/lC5DmRzmvwE5Esho9XZwAADdSUApjTvm+IXB7JUbkcxcscF6ndGNkydVj5cOSEyAciWfxy2jg3geQ/CwCgsZpSAK+K5Dz+7ZEcwauzKVJbG+lc27c38orIWZG7Ir8TyTJ4bQQAoLGaeh/A5ZI7P9wHEACGiE0gzdsEAgDAHBRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRTAIbFv/6Ot7bvGq0cAgH6MtB9ZmNHI/tAaHc3DpXHzjntbm7fe3XpiMhp7fMa2bNzQ2nTW2vZHAYD5mJiYaI2NjeVh/jGRB6UxAjjgcsSvLn8pH6/eutNIIACwYArggNs9fuBQ+asdnJxs7Rl/pH0GADA/CuCAW796VTXt22nFyEhr3eqV7TMAgPlRAAfcmrHjqjV/WfpSPl6z8YzqeQCAhbAJpD/Lsgkk5Zq/nPbNkT/lDwAWziYQI4BDI0vfOaedoPwBAH1TAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVpSgHcHNkReTjyQOSWyOmRXl0amYzk3wMAaLSmFMDzIx+MnB25KHJU5LbIqshcTo38l8id1dmA27f/0db2XePVIwDAQoy0H5vmxEiOBGYxvCOf6GJF5PORj0XOizw98rpIr0Yj+0NrdDQPl9bNO+5tbd56d+uJyWju8ZnbsnFDa9NZa9sfBQB6MTEx0RobG8vD/GMiD0rT1DWA1Wc1PNh+7OY3I/8Q+Wh1NrdjItn06hwfWRY54leXv5SPV2/daSQQAJi3JhbAHNW8LrItsjOf6OKlkTdFrqzOepNrDfd35L7Istg9fuBQ+asdnJxs7Rl/pH0GANCbJhbAGyJnRi6rzmaWI3d/FMnyN55P9GhLJEcX65wcWRbrV6+qpn07rRgZaa1bvbJ9BgDQm6YVwOsjr41cEJltdO60yLrIpyP/2M4vRvLv5nF+fCaPRXKtQJ3cdbws1owdV635y9KX8vGajWdUzwMAzEdTNoHk/48sf5dEXhb5ZmQ2x0aeM3V4yHsiOTL4lsg3Io9H5rKsm0BSrvnLad8c+VP+AGD+bAJpTgH8UOTyyM9F7skn2nKdXr1L4sbI/ZFcxzeTj0cGehcwANA/BbA5U8BXRfKTeHtkX0c2RWp5v5Q1U4cAAOVqygjgkWIEEACGjBHA5t4HEACALhRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGARxi+/Y/2tq+a7x6BADo1Uj7kYUZjewPrdHRPFw+N++4t7V5692tJyajxcdnccvGDa1NZ61tfxQA6GZiYqI1NjaWh/nHRB6UxgjgEMoRv7r8pXy8eutOI4EAQE8UwCG0e/zAofJXOzg52doz/kj7DACgOwVwCK1fvaqa9u20YmSktW71yvYZAEB3CuAQWjN2XLXmL0tfysdrNp5RPQ8AMBebQPpzxDaBpFzzl9O+OfKn/AFAb2wCMQI41LL0nXPaCcofADAvCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKYAPkO4Js3zVePQIAzMVbwfXniL4VXLp5x72tzVvvbj0xGW0+Ppv5HsGbzlrb/igAMJ23gjMCONRyxK8ufykfr96600ggADArBXCI7R4/cKj81Q5OTrb2jD/SPgMAeCoFcIitX72qmvbttGJkpLVu9cr2GQDAUymAQ2zN2HHVmr8sfSkfr9l4RvU8AEA3NoH054hvAkm55i+nfXPkT/kDgNnZBGIEsBGy9J1z2gnKHwDQEwUQAKAwCiAAQGEUQACAwiiAAACFUQAbxHsCAwC9cBuY/gzEbWCS9wQGgN64DYwRwEbwnsAAwHwogA3gPYEBgPlQABvAewIDAPOhADaA9wQGAObDJpD+DMwmkOQ9gQFgbjaBGAFsFO8JDAD0QgEEACiMAggAUBgFsIG8IwgAMBubQPozUJtAkncEAYDZ2QRiBLBRvCMIANALBbBBvCMIANALBbBBvCMIANCLphTAzZEdkYcjD0RuiZwemc2VkTsjD7Xz15EXR4aWdwQBAHrRlE0gn43cFMkSeFTkvZENkedHDkRm8onIFyLbIz+I/LvIxsgLIvdHejFwm0CSdwQBgO5sAmnuLuATIzkSeH7kjnyiBysiORL4K5Eb84keDGQB7JRlMNcG5vSwMggACmBq6hrA6rMaHmw/9iIXyv1IZD5/Z6DlLWFeeu3nWpd/5EvVY54DADSxAOao5nWRbZGd+USPro3k1G+uBezmmEgO9dU5PjKQ3BIGAOimiQXwhsiZkcuqs97k+r98fa4BzPWA3eRmk/0duS8ykNwSBgDopmkF8PrIayMXRHotZ78RuTryishd+cQstkRyernOyZGB5JYwAEA3TSmAWXVy5C9H8C6M7I704q2Rd0ZeGflyPjGHxyK5WLRO3nZmILklDADQTVN2AX8ocnnk5yL35BNtOU1bL3rLnb25xi+ncVNO+747kn8vbwdT+347vRiKXcA57bvy6Ke1Djx+0G5gAIpnF3BzCuC01W6HXBH5+NRh6/bInsgbq7Op41OnDp/ktyLvmjqc08AXwJS7f+sNITktnCODm85a2/4oAJRFAWzufQCXy1CMAOYtYDo3hOR08La3X2AkEIAiKYDN3AVMB7uBAYDpFMCGsxsYAJhOAWy46buB8xP+pnPXVccAQJmsAezPUGwCSbkW8GPb9rT+YNu3bAYBoGjWABoBLEpd/lI+ems4ACiTAlgIm0EAgJoCWAibQQCAmgJYCJtBAICaTSD9GZpNIDWbQQAonU0gRgCLZDMIAJRNASxMt80gt961TwkEgEIogIWZaTNIes+tX6/eM/jmHfe2nwEAmkoBLMz0zSCdTAcDQBkUwALlho9tb7+g9e8vfl77mcPcGxAAmk8BLFSOBF585pqnTAfnF8TKo31ZAECT+UlfsJmmg5+IXPKh7dYCAkCDKYCFy+ngrW8+p9W5JNBaQABoNgWQ1oHHD7YmZ7g1zN/teah9BgA0iQJI11vD/OpNXzEVDAANpAByaC3g9C+GnAre/Md3t/78rm+bDgaABplh3Id5GLr3Ap5NFr1f+eRX2mdP5j2DAWgK7wVsBJAOLzr1GTNOBad6NPCre60LBIBhpwByyGzvEpLyFjGvc4sYABh6poD706gp4Fqu98sdwLkJJEf+psvfGv7kX7+k9cJTnjH1BAAMEVPARgCZQY4EvuaFz5pxY0gyEggAw00BpKvc8JEjfTPNCOd9A60JBIDhpAAyq5zmvdZIIAA0igLInHoZCXSvQAAYHjaB9KeRm0C6yZG+LHs58jeT/GK68rz1rSvOXV+tIwSAQWQTiALYr6IKYMo1fzntO/29gzspggAMMgXQFDDzVK8J7HavwJTd8Pfv3N166bWfsz4QAAaQEcD+FDcCWJvrXoG1/A3DPQMBGCRGAI0AskBz3SuwVu8U/r3P72pt3zVuowgADAAjgP0pdgSwU5a6j23b0/rInd+qpn9nY30gAEeaEUAFsF8KYAdFEIBhoAAqgP1SAGfQy07hWmcRTLvHD7TWr16lFAKwZBRABbBfCmAXufv36q07Wwd7aYGh/kLMVxsdBGApKYAKYL8UwFnklPCe8Udad93/vdZv/8Xfd72BdDeKIABLQQFUAPulAPZoPusDp3tafJW+7ZXPbW04eay16ugVrQOPHzRNDMCCKYAKYL8UwHnqpwh2qkcHLz5zjUIIwLwogApgvxTABaqL4B9s+1Z1I+n8QuynEHaOEiqDAMxGAVQA+6UA9qleJ7hu9crqfClGB00bA9BJAVQA+6UALoHpo4OLxVpCAJICqAD2SwFcQvXoYO4i/k+fuafnW8rMx/TRwiyEyf0IAZpLAVQA+6UALpO6DK48+mmtW+/6zqKPDtbq/yDyH60cAjSTAqgA9ksBPEI6C+Ejjz+x4HsNzsds5dCUMsDwUAAVwH4pgANkqdYOzsdM6wynPxpJBDiyFEAFsF8K4ACaPjrYOUq4VGsJ56PXkcSkJAIsPgVQAeyXAjhk5lpLmP9BHNl6OGW+JXGmUUalEWBmCqAC2C8FcMhNHy3svB/hoJbDWmdJnG6+U9EKJFASBVAB7JcC2GDzKYfDaLYCmR/LUccrzl1fnc9UEnstkvXHlEpgUCiACmC/FMBC9bPOMP+jG5buOFtJrM1VJFN+LI/nO5W90Ncom8BsFMDD359ZGAWQp+hWDuczkpj/YQ5LSezXbAWyNt/X5HFn2VyKkjnba5RPGGwK4OHvmSyMAkhfupXFXkti/nGENzUPhc5y2M1ivWa29ZfLXUS7vWaujymwNJ0CePj7GQujALLk5iqJ/UxF5x8K5PLppUAux2vm+thirP/s5TX9/v3leE0+pxA3jwJ4+PsAC6MAMrB6nYqe/prO2+NU3yDij9lK4myvyY/N8lcZYNXnNcz2+ev3Ncvx76gt9DX53HIV4l5esxz/jl5es9z/Oxa7hCuAh7/eWRgFkEaqy+Nso4xzFcnOj/U7ld3ra2b74Q79qL4Gw2xfY8vxmhL/d+Syii0bN7Q2nbW2/Uz/FMDD137YbY5sjDw38mhke+RtkXsis/n5yLsjp0V2Rd4R+ZNIrxRA6NFCprIX8ppu93DMPxajZM71GmDxrRgZaW17+wWLNhKoALa/nzXAZyM3RXZEjoq8N7Ih8vzIgchMzoncGXlnJEvfJZH/GDk38qVILxRAGEDTy+Zil8zZXjPb+svFKpn9vqaXvw+D5lNXnt0657QT2mf9UQDb3wca6MTIA5HzI3fkEzO4OZKt7VXV2ZQskg9FLqvO5qYAAk+xXKOdC31Nt48t5vrP2myv6ffv15bzNRwZRgAXX/X13kDPiXwzkqOAO/OJGdwbeV87tV+P/Frk1Opsbgog0CiLtf6zl9csx7+jn9csdyGuzfaafv9+rd/X9Pv3a728JsvfNRvPsAZwkVXXvmHy/9OfRp4ROS+f6OLxyBsjn6zOplwe+VjkmOrsqfL5zo8dH7lPAQRopuUsxL28ptT/HXYBL74mFsAPRi6O5Fq++/KJLrIA/lLkU9XZlNdHPho5tjp7qndF/sPU4WEKIAAMDwWw1Xpa+7Epro+8NnJBZLbyl74TOWnq8JBnRr47dTijLZH8YqlzcgQAYKg0pQDmSOYNkbwVzIWR3ZG5fDFy0dThIa+I5C1kunkskr8p1Hk4AgAwVJpSAHPa9w2RXMOXpSxH9jKdiwZujOQIXu0DkSx8eb/AvH9gPr488v4IAEBjNaUAXhXJKdnbI/s6silSy+1Da6YOKznSd2nkishdkdwQkq/v9R6AAABDqYmbQJaT28AAwJCxCaR5m0AAAJiDAggAUBgFEACgMAogAEBhFEAAgMIogAAAhXEbmP5Ut4HZu3ev28AAwJDI28CccsopeVjsbWAUwP78eGSu9xwGAAZTvqf//VOHZVEA+5PX71mRpXhP4OMjWS7zi9N7Di8d13l5uM7Lw3VeHq7z8ljq65z//G9HJqszGBA5p5xflOaWl5brvDxc5+XhOi8P13l5uM5LyCYQAIDCKIAAAIVZ0X5kMB2M3N5+ZOm4zsvDdV4ervPycJ2Xh+sMAAAAAAAAAAAAAAAADJY3R3ZHfhD5u8h5ERZuc2RHJO8k/0DklsjpkU7HRK6PjEcORP4sknefZ+HyuudNXN9fnU1xnRdHvg3lH0X+X+SRyP+OvChSy3cpelck3+Xg0UjuonxBhN4dFXlPJL8X5zX8VuQ3I523T3OdF+ZnIp+O1O/C8bpIp16u6zMifxjZ304ePz0CQ2tT5PHIv4w8L5I/PL8fWRthYT4beWMkv4G8MPLnkf8bWRWpfTiSbzn08sg/jXwukj9U3SppYc6K5A/Or0Y6C6Dr3L/8wbcn8rHIiyPrIj8bOS1Se1sk3+B+Y+SMyE2R/GGab31Fb94RyV9ULo7kNf6FSP4S+ZZIzXVemFdFslzndZupAPZyXT8TuTtyTjt5nKUShtaXIvlDstPXI1umDlkEJ0bym07+FprGIlm6s3zX8j2e875T/7w6Yz5+NPKNSJa8/M29LoCu8+K4NnLn1OGMcvRkXyR/iNZy5PV7kV+uzuhF/qL40anDQ/44kiNNyXVeHNMLYC/XNQdH8u/9dHU25exIPjd9docuvBPIYDk6ktM4t1Vnh+X5S6YOWQRZRNKD7ce85j8S6bzu+dvmzojrPn8fjNwa+evq7DDXeXG8NvLlyP+I5JKGr0SujNTWR06KdF7nxyKfj7jOvdsWyZHVn6jOpmYPzo38RXXmOi+VXq5rjvjltG8OmNT+NpLPufY9UgAHy+pIToV9tzo7LM/zPwj6l79dXhfJb+5ZPFJe2xyZeqg6O8x1n79LI1n0cv3fdK7z4nh25KrINyM5cvq7kd+J/GIk1dfS95H+/HbkU5G/j/wwkkU7R7PzueQ6L41erms+5i8/0+Vzrn2PFMDBlMPYnbK0TH+Ohbkhcmbksupsdq77/JwS+UDk9ZHcwNQr13l+8vv2/4pcHclS8nuRj0SyFHbyfaQ/uVThDZHLIz8V+aXIb7QfO7nOS2Ou6zrTNXbt50EBHCy54DjXQ03/DeaZkem/DTF/ufs0p88uiORGhNp3Ijn9novrO7nu85Mjf3nNcuf6P7ZzfuRX28d5LV3n/uX6qK9NHR6S64TrjWL59Zx8H+nPf47kesvcgJAbDHLt3/si9ei267w0ermu+Zofmzp8klzf7dr3SAEcLDk9lj88L6rODsvz7VOHLED+Vpgjf7mj7MJI7k7tlNc8p3g6r/uaSO4+c9179zeRDZGf7EiuVftEx7Hr3L8vRKYvdM91armzPeXXd/6A7LzOWbyzjLvOvVsZeWLq8JD8Bb3+uek6L41erusXI7mWO3fB13JDSD7n2jO06tvA/ItI7nTK3zjzNjCnRliYD0VyB1l+A8nfKuscF6nlzuu9kVz0nbcnyTLj9iT969wFnFzn/uUtdrJI5xTwcyI5RZn3VMyp91ruoMyv+UsiWbA/GXF7kvn5eCRnCurbwOS1/IdIrg2suc4Lk3cKqH9JzCnbX28f16PYvVzXvA1M3mYqd/9m7oq4DQxDL28Enff5yp1POTpV366EhclvMDMl7w1YOzaSU8T1jXXzG0muaaM/0wug67w4XhPJaclca5nTv527gFOOeueNdHO6OF+TOyjzBym9y7KRX7s5spo3I94VyXvX5WhUzXVemJdFZvqenKU79XJd/0kkb4ae9wvM5LEbQQMAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADQMK3W/weZS0BCEiWYiwAAAABJRU5ErkJggg==\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "temp = np.linspace(0, 100, 200)\n",
    "new_data = {\"Temperature (K)\": temp, \"log Ohm\": 1/(temp+1)+2}\n",
    "\n",
    "fig, ax = plt.subplots()\n",
    "ax.plot(new_data[\"Temperature (K)\"], new_data[\"log Ohm\"], '.')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "curve_x.format(\"log Ohm/K\")\n",
    "curve_x.set_data(new_data)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'log Ohm/K'"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "curve_x.format()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "curve_x_data = curve_x.get_data()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": "/* Put everything inside the global mpl namespace */\nwindow.mpl = {};\n\n\nmpl.get_websocket_type = function() {\n    if (typeof(WebSocket) !== 'undefined') {\n        return WebSocket;\n    } else if (typeof(MozWebSocket) !== 'undefined') {\n        return MozWebSocket;\n    } else {\n        alert('Your browser does not have WebSocket support.' +\n              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n              'Firefox 4 and 5 are also supported but you ' +\n              'have to enable WebSockets in about:config.');\n    };\n}\n\nmpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n    this.id = figure_id;\n\n    this.ws = websocket;\n\n    this.supports_binary = (this.ws.binaryType != undefined);\n\n    if (!this.supports_binary) {\n        var warnings = document.getElementById(\"mpl-warnings\");\n        if (warnings) {\n            warnings.style.display = 'block';\n            warnings.textContent = (\n                \"This browser does not support binary websocket messages. \" +\n                    \"Performance may be slow.\");\n        }\n    }\n\n    this.imageObj = new Image();\n\n    this.context = undefined;\n    this.message = undefined;\n    this.canvas = undefined;\n    this.rubberband_canvas = undefined;\n    this.rubberband_context = undefined;\n    this.format_dropdown = undefined;\n\n    this.image_mode = 'full';\n\n    this.root = $('<div/>');\n    this._root_extra_style(this.root)\n    this.root.attr('style', 'display: inline-block');\n\n    $(parent_element).append(this.root);\n\n    this._init_header(this);\n    this._init_canvas(this);\n    this._init_toolbar(this);\n\n    var fig = this;\n\n    this.waiting = false;\n\n    this.ws.onopen =  function () {\n            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n            fig.send_message(\"send_image_mode\", {});\n            if (mpl.ratio != 1) {\n                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n            }\n            fig.send_message(\"refresh\", {});\n        }\n\n    this.imageObj.onload = function() {\n            if (fig.image_mode == 'full') {\n                // Full images could contain transparency (where diff images\n                // almost always do), so we need to clear the canvas so that\n                // there is no ghosting.\n                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n            }\n            fig.context.drawImage(fig.imageObj, 0, 0);\n        };\n\n    this.imageObj.onunload = function() {\n        this.ws.close();\n    }\n\n    this.ws.onmessage = this._make_on_message_function(this);\n\n    this.ondownload = ondownload;\n}\n\nmpl.figure.prototype._init_header = function() {\n    var titlebar = $(\n        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n        'ui-helper-clearfix\"/>');\n    var titletext = $(\n        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n        'text-align: center; padding: 3px;\"/>');\n    titlebar.append(titletext)\n    this.root.append(titlebar);\n    this.header = titletext[0];\n}\n\n\n\nmpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n\n}\n\n\nmpl.figure.prototype._root_extra_style = function(canvas_div) {\n\n}\n\nmpl.figure.prototype._init_canvas = function() {\n    var fig = this;\n\n    var canvas_div = $('<div/>');\n\n    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n\n    function canvas_keyboard_event(event) {\n        return fig.key_event(event, event['data']);\n    }\n\n    canvas_div.keydown('key_press', canvas_keyboard_event);\n    canvas_div.keyup('key_release', canvas_keyboard_event);\n    this.canvas_div = canvas_div\n    this._canvas_extra_style(canvas_div)\n    this.root.append(canvas_div);\n\n    var canvas = $('<canvas/>');\n    canvas.addClass('mpl-canvas');\n    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n\n    this.canvas = canvas[0];\n    this.context = canvas[0].getContext(\"2d\");\n\n    var backingStore = this.context.backingStorePixelRatio ||\n\tthis.context.webkitBackingStorePixelRatio ||\n\tthis.context.mozBackingStorePixelRatio ||\n\tthis.context.msBackingStorePixelRatio ||\n\tthis.context.oBackingStorePixelRatio ||\n\tthis.context.backingStorePixelRatio || 1;\n\n    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n\n    var rubberband = $('<canvas/>');\n    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n\n    var pass_mouse_events = true;\n\n    canvas_div.resizable({\n        start: function(event, ui) {\n            pass_mouse_events = false;\n        },\n        resize: function(event, ui) {\n            fig.request_resize(ui.size.width, ui.size.height);\n        },\n        stop: function(event, ui) {\n            pass_mouse_events = true;\n            fig.request_resize(ui.size.width, ui.size.height);\n        },\n    });\n\n    function mouse_event_fn(event) {\n        if (pass_mouse_events)\n            return fig.mouse_event(event, event['data']);\n    }\n\n    rubberband.mousedown('button_press', mouse_event_fn);\n    rubberband.mouseup('button_release', mouse_event_fn);\n    // Throttle sequential mouse events to 1 every 20ms.\n    rubberband.mousemove('motion_notify', mouse_event_fn);\n\n    rubberband.mouseenter('figure_enter', mouse_event_fn);\n    rubberband.mouseleave('figure_leave', mouse_event_fn);\n\n    canvas_div.on(\"wheel\", function (event) {\n        event = event.originalEvent;\n        event['data'] = 'scroll'\n        if (event.deltaY < 0) {\n            event.step = 1;\n        } else {\n            event.step = -1;\n        }\n        mouse_event_fn(event);\n    });\n\n    canvas_div.append(canvas);\n    canvas_div.append(rubberband);\n\n    this.rubberband = rubberband;\n    this.rubberband_canvas = rubberband[0];\n    this.rubberband_context = rubberband[0].getContext(\"2d\");\n    this.rubberband_context.strokeStyle = \"#000000\";\n\n    this._resize_canvas = function(width, height) {\n        // Keep the size of the canvas, canvas container, and rubber band\n        // canvas in synch.\n        canvas_div.css('width', width)\n        canvas_div.css('height', height)\n\n        canvas.attr('width', width * mpl.ratio);\n        canvas.attr('height', height * mpl.ratio);\n        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n\n        rubberband.attr('width', width);\n        rubberband.attr('height', height);\n    }\n\n    // Set the figure to an initial 600x600px, this will subsequently be updated\n    // upon first draw.\n    this._resize_canvas(600, 600);\n\n    // Disable right mouse context menu.\n    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n        return false;\n    });\n\n    function set_focus () {\n        canvas.focus();\n        canvas_div.focus();\n    }\n\n    window.setTimeout(set_focus, 100);\n}\n\nmpl.figure.prototype._init_toolbar = function() {\n    var fig = this;\n\n    var nav_element = $('<div/>')\n    nav_element.attr('style', 'width: 100%');\n    this.root.append(nav_element);\n\n    // Define a callback function for later on.\n    function toolbar_event(event) {\n        return fig.toolbar_button_onclick(event['data']);\n    }\n    function toolbar_mouse_event(event) {\n        return fig.toolbar_button_onmouseover(event['data']);\n    }\n\n    for(var toolbar_ind in mpl.toolbar_items) {\n        var name = mpl.toolbar_items[toolbar_ind][0];\n        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n        var image = mpl.toolbar_items[toolbar_ind][2];\n        var method_name = mpl.toolbar_items[toolbar_ind][3];\n\n        if (!name) {\n            // put a spacer in here.\n            continue;\n        }\n        var button = $('<button/>');\n        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n                        'ui-button-icon-only');\n        button.attr('role', 'button');\n        button.attr('aria-disabled', 'false');\n        button.click(method_name, toolbar_event);\n        button.mouseover(tooltip, toolbar_mouse_event);\n\n        var icon_img = $('<span/>');\n        icon_img.addClass('ui-button-icon-primary ui-icon');\n        icon_img.addClass(image);\n        icon_img.addClass('ui-corner-all');\n\n        var tooltip_span = $('<span/>');\n        tooltip_span.addClass('ui-button-text');\n        tooltip_span.html(tooltip);\n\n        button.append(icon_img);\n        button.append(tooltip_span);\n\n        nav_element.append(button);\n    }\n\n    var fmt_picker_span = $('<span/>');\n\n    var fmt_picker = $('<select/>');\n    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n    fmt_picker_span.append(fmt_picker);\n    nav_element.append(fmt_picker_span);\n    this.format_dropdown = fmt_picker[0];\n\n    for (var ind in mpl.extensions) {\n        var fmt = mpl.extensions[ind];\n        var option = $(\n            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n        fmt_picker.append(option)\n    }\n\n    // Add hover states to the ui-buttons\n    $( \".ui-button\" ).hover(\n        function() { $(this).addClass(\"ui-state-hover\");},\n        function() { $(this).removeClass(\"ui-state-hover\");}\n    );\n\n    var status_bar = $('<span class=\"mpl-message\"/>');\n    nav_element.append(status_bar);\n    this.message = status_bar[0];\n}\n\nmpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n    // which will in turn request a refresh of the image.\n    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n}\n\nmpl.figure.prototype.send_message = function(type, properties) {\n    properties['type'] = type;\n    properties['figure_id'] = this.id;\n    this.ws.send(JSON.stringify(properties));\n}\n\nmpl.figure.prototype.send_draw_message = function() {\n    if (!this.waiting) {\n        this.waiting = true;\n        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n    }\n}\n\n\nmpl.figure.prototype.handle_save = function(fig, msg) {\n    var format_dropdown = fig.format_dropdown;\n    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n    fig.ondownload(fig, format);\n}\n\n\nmpl.figure.prototype.handle_resize = function(fig, msg) {\n    var size = msg['size'];\n    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n        fig._resize_canvas(size[0], size[1]);\n        fig.send_message(\"refresh\", {});\n    };\n}\n\nmpl.figure.prototype.handle_rubberband = function(fig, msg) {\n    var x0 = msg['x0'] / mpl.ratio;\n    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n    var x1 = msg['x1'] / mpl.ratio;\n    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n    x0 = Math.floor(x0) + 0.5;\n    y0 = Math.floor(y0) + 0.5;\n    x1 = Math.floor(x1) + 0.5;\n    y1 = Math.floor(y1) + 0.5;\n    var min_x = Math.min(x0, x1);\n    var min_y = Math.min(y0, y1);\n    var width = Math.abs(x1 - x0);\n    var height = Math.abs(y1 - y0);\n\n    fig.rubberband_context.clearRect(\n        0, 0, fig.canvas.width, fig.canvas.height);\n\n    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n}\n\nmpl.figure.prototype.handle_figure_label = function(fig, msg) {\n    // Updates the figure title.\n    fig.header.textContent = msg['label'];\n}\n\nmpl.figure.prototype.handle_cursor = function(fig, msg) {\n    var cursor = msg['cursor'];\n    switch(cursor)\n    {\n    case 0:\n        cursor = 'pointer';\n        break;\n    case 1:\n        cursor = 'default';\n        break;\n    case 2:\n        cursor = 'crosshair';\n        break;\n    case 3:\n        cursor = 'move';\n        break;\n    }\n    fig.rubberband_canvas.style.cursor = cursor;\n}\n\nmpl.figure.prototype.handle_message = function(fig, msg) {\n    fig.message.textContent = msg['message'];\n}\n\nmpl.figure.prototype.handle_draw = function(fig, msg) {\n    // Request the server to send over a new figure.\n    fig.send_draw_message();\n}\n\nmpl.figure.prototype.handle_image_mode = function(fig, msg) {\n    fig.image_mode = msg['mode'];\n}\n\nmpl.figure.prototype.updated_canvas_event = function() {\n    // Called whenever the canvas gets updated.\n    this.send_message(\"ack\", {});\n}\n\n// A function to construct a web socket function for onmessage handling.\n// Called in the figure constructor.\nmpl.figure.prototype._make_on_message_function = function(fig) {\n    return function socket_on_message(evt) {\n        if (evt.data instanceof Blob) {\n            /* FIXME: We get \"Resource interpreted as Image but\n             * transferred with MIME type text/plain:\" errors on\n             * Chrome.  But how to set the MIME type?  It doesn't seem\n             * to be part of the websocket stream */\n            evt.data.type = \"image/png\";\n\n            /* Free the memory for the previous frames */\n            if (fig.imageObj.src) {\n                (window.URL || window.webkitURL).revokeObjectURL(\n                    fig.imageObj.src);\n            }\n\n            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n                evt.data);\n            fig.updated_canvas_event();\n            fig.waiting = false;\n            return;\n        }\n        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n            fig.imageObj.src = evt.data;\n            fig.updated_canvas_event();\n            fig.waiting = false;\n            return;\n        }\n\n        var msg = JSON.parse(evt.data);\n        var msg_type = msg['type'];\n\n        // Call the  \"handle_{type}\" callback, which takes\n        // the figure and JSON message as its only arguments.\n        try {\n            var callback = fig[\"handle_\" + msg_type];\n        } catch (e) {\n            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n            return;\n        }\n\n        if (callback) {\n            try {\n                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n                callback(fig, msg);\n            } catch (e) {\n                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n            }\n        }\n    };\n}\n\n// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\nmpl.findpos = function(e) {\n    //this section is from http://www.quirksmode.org/js/events_properties.html\n    var targ;\n    if (!e)\n        e = window.event;\n    if (e.target)\n        targ = e.target;\n    else if (e.srcElement)\n        targ = e.srcElement;\n    if (targ.nodeType == 3) // defeat Safari bug\n        targ = targ.parentNode;\n\n    // jQuery normalizes the pageX and pageY\n    // pageX,Y are the mouse positions relative to the document\n    // offset() returns the position of the element relative to the document\n    var x = e.pageX - $(targ).offset().left;\n    var y = e.pageY - $(targ).offset().top;\n\n    return {\"x\": x, \"y\": y};\n};\n\n/*\n * return a copy of an object with only non-object keys\n * we need this to avoid circular references\n * http://stackoverflow.com/a/24161582/3208463\n */\nfunction simpleKeys (original) {\n  return Object.keys(original).reduce(function (obj, key) {\n    if (typeof original[key] !== 'object')\n        obj[key] = original[key]\n    return obj;\n  }, {});\n}\n\nmpl.figure.prototype.mouse_event = function(event, name) {\n    var canvas_pos = mpl.findpos(event)\n\n    if (name === 'button_press')\n    {\n        this.canvas.focus();\n        this.canvas_div.focus();\n    }\n\n    var x = canvas_pos.x * mpl.ratio;\n    var y = canvas_pos.y * mpl.ratio;\n\n    this.send_message(name, {x: x, y: y, button: event.button,\n                             step: event.step,\n                             guiEvent: simpleKeys(event)});\n\n    /* This prevents the web browser from automatically changing to\n     * the text insertion cursor when the button is pressed.  We want\n     * to control all of the cursor setting manually through the\n     * 'cursor' event from matplotlib */\n    event.preventDefault();\n    return false;\n}\n\nmpl.figure.prototype._key_event_extra = function(event, name) {\n    // Handle any extra behaviour associated with a key event\n}\n\nmpl.figure.prototype.key_event = function(event, name) {\n\n    // Prevent repeat events\n    if (name == 'key_press')\n    {\n        if (event.which === this._key)\n            return;\n        else\n            this._key = event.which;\n    }\n    if (name == 'key_release')\n        this._key = null;\n\n    var value = '';\n    if (event.ctrlKey && event.which != 17)\n        value += \"ctrl+\";\n    if (event.altKey && event.which != 18)\n        value += \"alt+\";\n    if (event.shiftKey && event.which != 16)\n        value += \"shift+\";\n\n    value += 'k';\n    value += event.which.toString();\n\n    this._key_event_extra(event, name);\n\n    this.send_message(name, {key: value,\n                             guiEvent: simpleKeys(event)});\n    return false;\n}\n\nmpl.figure.prototype.toolbar_button_onclick = function(name) {\n    if (name == 'download') {\n        this.handle_save(this, null);\n    } else {\n        this.send_message(\"toolbar_button\", {name: name});\n    }\n};\n\nmpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n    this.message.textContent = tooltip;\n};\nmpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n\nmpl.extensions = [\"eps\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\"];\n\nmpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n    // Create a \"websocket\"-like object which calls the given IPython comm\n    // object with the appropriate methods. Currently this is a non binary\n    // socket, so there is still some room for performance tuning.\n    var ws = {};\n\n    ws.close = function() {\n        comm.close()\n    };\n    ws.send = function(m) {\n        //console.log('sending', m);\n        comm.send(m);\n    };\n    // Register the callback with on_msg.\n    comm.on_msg(function(msg) {\n        //console.log('receiving', msg['content']['data'], msg);\n        // Pass the mpl event to the overriden (by mpl) onmessage function.\n        ws.onmessage(msg['content']['data'])\n    });\n    return ws;\n}\n\nmpl.mpl_figure_comm = function(comm, msg) {\n    // This is the function which gets called when the mpl process\n    // starts-up an IPython Comm through the \"matplotlib\" channel.\n\n    var id = msg.content.data.id;\n    // Get hold of the div created by the display call when the Comm\n    // socket was opened in Python.\n    var element = $(\"#\" + id);\n    var ws_proxy = comm_websocket_adapter(comm)\n\n    function ondownload(figure, format) {\n        window.open(figure.imageObj.src);\n    }\n\n    var fig = new mpl.figure(id, ws_proxy,\n                           ondownload,\n                           element.get(0));\n\n    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n    // web socket which is closed, not our websocket->open comm proxy.\n    ws_proxy.onopen();\n\n    fig.parent_element = element.get(0);\n    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n    if (!fig.cell_info) {\n        console.error(\"Failed to find cell for figure\", id, fig);\n        return;\n    }\n\n    var output_index = fig.cell_info[2]\n    var cell = fig.cell_info[0];\n\n};\n\nmpl.figure.prototype.handle_close = function(fig, msg) {\n    var width = fig.canvas.width/mpl.ratio\n    fig.root.unbind('remove')\n\n    // Update the output cell to use the data from the current canvas.\n    fig.push_to_output();\n    var dataURL = fig.canvas.toDataURL();\n    // Re-enable the keyboard manager in IPython - without this line, in FF,\n    // the notebook keyboard shortcuts fail.\n    IPython.keyboard_manager.enable()\n    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n    fig.close_ws(fig, msg);\n}\n\nmpl.figure.prototype.close_ws = function(fig, msg){\n    fig.send_message('closing', msg);\n    // fig.ws.close()\n}\n\nmpl.figure.prototype.push_to_output = function(remove_interactive) {\n    // Turn the data on the canvas into data in the output cell.\n    var width = this.canvas.width/mpl.ratio\n    var dataURL = this.canvas.toDataURL();\n    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n}\n\nmpl.figure.prototype.updated_canvas_event = function() {\n    // Tell IPython that the notebook contents must change.\n    IPython.notebook.set_dirty(true);\n    this.send_message(\"ack\", {});\n    var fig = this;\n    // Wait a second, then push the new image to the DOM so\n    // that it is saved nicely (might be nice to debounce this).\n    setTimeout(function () { fig.push_to_output() }, 1000);\n}\n\nmpl.figure.prototype._init_toolbar = function() {\n    var fig = this;\n\n    var nav_element = $('<div/>')\n    nav_element.attr('style', 'width: 100%');\n    this.root.append(nav_element);\n\n    // Define a callback function for later on.\n    function toolbar_event(event) {\n        return fig.toolbar_button_onclick(event['data']);\n    }\n    function toolbar_mouse_event(event) {\n        return fig.toolbar_button_onmouseover(event['data']);\n    }\n\n    for(var toolbar_ind in mpl.toolbar_items){\n        var name = mpl.toolbar_items[toolbar_ind][0];\n        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n        var image = mpl.toolbar_items[toolbar_ind][2];\n        var method_name = mpl.toolbar_items[toolbar_ind][3];\n\n        if (!name) { continue; };\n\n        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n        button.click(method_name, toolbar_event);\n        button.mouseover(tooltip, toolbar_mouse_event);\n        nav_element.append(button);\n    }\n\n    // Add the status bar.\n    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n    nav_element.append(status_bar);\n    this.message = status_bar[0];\n\n    // Add the close button to the window.\n    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n    button.click(function (evt) { fig.handle_close(fig, {}); } );\n    button.mouseover('Stop Interaction', toolbar_mouse_event);\n    buttongrp.append(button);\n    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n    titlebar.prepend(buttongrp);\n}\n\nmpl.figure.prototype._root_extra_style = function(el){\n    var fig = this\n    el.on(\"remove\", function(){\n\tfig.close_ws(fig, {});\n    });\n}\n\nmpl.figure.prototype._canvas_extra_style = function(el){\n    // this is important to make the div 'focusable\n    el.attr('tabindex', 0)\n    // reach out to IPython and tell the keyboard manager to turn it's self\n    // off when our div gets focus\n\n    // location in version 3\n    if (IPython.notebook.keyboard_manager) {\n        IPython.notebook.keyboard_manager.register_events(el);\n    }\n    else {\n        // location in version 2\n        IPython.keyboard_manager.register_events(el);\n    }\n\n}\n\nmpl.figure.prototype._key_event_extra = function(event, name) {\n    var manager = IPython.notebook.keyboard_manager;\n    if (!manager)\n        manager = IPython.keyboard_manager;\n\n    // Check for shift+enter\n    if (event.shiftKey && event.which == 13) {\n        this.canvas_div.blur();\n        // select the cell after this one\n        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n        IPython.notebook.select(index + 1);\n    }\n}\n\nmpl.figure.prototype.handle_save = function(fig, msg) {\n    fig.ondownload(fig, null);\n}\n\n\nmpl.find_output_cell = function(html_output) {\n    // Return the cell and output element which can be found *uniquely* in the notebook.\n    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n    // IPython event is triggered only after the cells have been serialised, which for\n    // our purposes (turning an active figure into a static one), is too late.\n    var cells = IPython.notebook.get_cells();\n    var ncells = cells.length;\n    for (var i=0; i<ncells; i++) {\n        var cell = cells[i];\n        if (cell.cell_type === 'code'){\n            for (var j=0; j<cell.output_area.outputs.length; j++) {\n                var data = cell.output_area.outputs[j];\n                if (data.data) {\n                    // IPython >= 3 moved mimebundle to data attribute of output\n                    data = data.data;\n                }\n                if (data['text/html'] == html_output) {\n                    return [cell, data, j];\n                }\n            }\n        }\n    }\n}\n\n// Register the function which deals with the matplotlib target/channel.\n// The kernel may be null if the page has been refreshed.\nif (IPython.notebook.kernel != null) {\n    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n}\n",
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAB/oSURBVHhe7d0NlKV1fR/wu13Cy67MSFkImwDuiif4AibGrmFVghCxUXJQNjlnYfUkpgknxeZoX2LNYkzTiiw9bYnKi2k9HFMSFfqywVgbJcYSwDWe1Vphi9WELmEhK2QK2cEFwcD093vmPst1uHfmztyZu/c+/8/nnO/e57lzLyzPDDPf+f+f5/+0AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAaJZV7UeWJo/fD0Ueq/YAgHFxbOSvIjPVXmEUwMH8cOSB2U0AYMycHHlwdrMsCuBgJiIH9u3b15qYyE0AYNRNT0+3TjnllNyczN3cKI0COJiqAAYFEADGRBbAycnsfuUWwL/TfgQAoBAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgME0pgJdF7orkat6ZL0XeGJnPz0buiTzZfrwoAgDQeE0pgA9Efj3y99r5QuRTkZdFutkcuTnye5EfbT/+p8hPREbC/gNPtHbdO1U9AgAspybfC/iRyLsjN1R73y/LX968t3OU8LORRyOXVHv9WZF7Ad+8+/7W9p13t56ZiYYen6EdW85sbd10avujAMAg3Au4mecAro5cHFkbyangbnIE8NbZzUM+F3n17GZPR0Wy6dU5NrKscsSvLn8pHy/fucdIIACwbJpUAM+MfCeS5/T9TiTP6ctz+7o5KfLQ7OYhuZ/Pz2d75EBHcup5We2dOnio/NWenplp3Tf1eHsPAGAwTSqA34z8WOSsyEci/zHy0kgvc2pWNR0+97m5dkRyuLjOyZFltXHd2mrat9PqVataG9atae8BAAymSQXwqchfRL4SyZG6r0feFenm25G5o30nRuaOCs6Vo4v1lcaZxyLLav3kMdU5f1n6Uj5eueWM6nkAgOXQ5ItA/iSyL/L2au/75UUgef7em6q9WX8U+ZvIYb8IJOU5fzntmyN/yh8ALB8XgTRnBPDKyNmRDZE8F/ADkddFPh5JN0Zy+rb2ocgbIu+JvLj9+PrIByMjIUvf5tOOV/4AgGXXlAL4g5Fcyy/PA8yRv1zP76cjfxxJuYbK+tnNyq5IXin8i5FcQDpHCbdGvhwBAGi0Jk8BD8OKTQEDACvDFHCzLgIBAKAPCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDBNKYDbI7sjj0UejtwSOT2ykH8c+Wbkici+yG9Hjo4AADRWUwrgOZHrImdFzo8cEbk1sjbSy1sjV0X+ZeQlkV+KbI3siAAANNaq9mPTnBDJkcAshrfnE11cG8ni91PV3qx/F3lV5Oxqb2ETkQOhNTGRmwDAqJuenm5NTk7mZv4xnRulaeo5gNVnNTzSfuzmzsgrI1n40gsjb4p8ptrr7qhINr06x0YAAMZKEwtgjmpeHcmCtyef6OGmyPsi+brvRe6N/I9ITgv3kucaHujIAxEAgLHSxAKYU7svj1xS7fX2ush7I++I/HhkS+RnIlkKe8nzA3N0sc7JEQCAsdK0cwCvibwl8pORvfnEPO6I/Fnk3dXerLdF/kPkeZFn8okFOAcQAMaMcwCbMwKYRTZH/nIU77zIQuUvrYnMLXlPR/Kf1dSLYwAAGlMAcwmYHL3bFsm1AE9q55hI7cZI5xIvn45cFrk4sjGSy8e8P/KHkSyCAACN1JQCmEUuh3Fvi+zvSK7rVzs1sn52s3JFJJd9ycd7IjdEPhf5lQgAQGOZ6hyMcwABYMw4B7CZVwEDADAPBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMArgiNt/4InWrnunqkcAgOWwqv3I0kxEDoTWxERuLq+bd9/f2r7z7tYzM9HU4zO1Y8uZra2bTm1/FABYiunp6dbk5GRu5h/TuVEaI4AjKkf86vKX8vHynXuMBAIAA1MAR9TeqYOHyl/t6ZmZ1n1Tj7f3AACWRgEcURvXra2mfTutXrWqtWHdmvYeAMDSKIAjav3kMdU5f1n6Uj5eueWM6nkAgEG4CGQwK3oRSMpz/nLaN0f+lD8AGJyLQIwAjrwsfZtPO175AwCWjQIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABSmKQVwe2R35LHIw5FbIqdHFvL8yHWR/ZHvRr4ReVMEAKCxmlIAz4lkkTsrcn7kiMitkbWRXo6M/HFkQ+TnIlkYL408GAEAaKxV7cemOSGSI4FZDG/PJ7r4h5F3R14c+V4+sQQTkQOhNTGRmwDAqJuenm5NTk7mZv4xnRulaeo5gNVnNTzSfuzmwsiXIjly+FBkT+TyyOpIL0dFsunVOTYCADBWmlgAc1Tz6sidkSx1vbwwklO/WfjyvL8rIv8s8t5IL3mu4YGOPBABABgrTZwCzhG9CyKvjcxX0L4VOTqyMfJ0PhH+aSSnhddXe8+VI4CZWo4APmAKGADGhyng5o0AXhPJqd1zIwuNzuWVv1kC6/KX8irgkyJ5gUg3T0byC6VOXnUMADBWmlIAcyTz2siWyHmRvZGFfDHyokjnMfiRSBbDp6o9AIAGakoBzGnft0W2RXJULkfxMsdEajdGdsxuVj4SOT7yoUgWv5w2zotA8p8FANBYTSmAl0VyHv+2SI7g1dkaqZ0a6Ty3b1/kDZFNkbsiH45kGbwqAgDQWE1dB3BY8soP6wACwBhxEUjzLgIBAGABCiAAQGEUQACAwiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAVwTOw/8ERr171T1SMAwCBWtR9ZmonIgdCamMjNlXHz7vtb23fe3XpmJhp7fMZ2bDmztXXTqe2PAgCLMT093ZqcnMzN/GM6N0pjBHDE5YhfXf5SPl6+c4+RQABgyRTAEbd36uCh8ld7emamdd/U4+09AIDFUQBH3MZ1a6tp306rV61qbVi3pr0HALA4CuCIWz95THXOX5a+lI9Xbjmjeh4AYClcBDKYoVwEkvKcv5z2zZE/5Q8Als5FIEYAx0aWvs2nHa/8AQADUwABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGAQQAKIwCCABQGAUQAKAwCiAAQGEUQACAwiiAAACFUQABAArTlAK4PbI78ljk4cgtkdMj/bo4MhPJ9wEANFpTCuA5kesiZ0XOjxwRuTWyNrKQF0T+beSOam/E7T/wRGvXvVPVIwDAUqxqPzbNCZEcCcxieHs+0cPqyJ9GPhY5O/L8yFsi/ZqIHAitiYncXFk3776/tX3n3a1nZqK5x2dux5YzW1s3ndr+KADQj+np6dbk5GRu5h/TuVGapp4DWH1WwyPtx15+M/LXkRuqvYUdFcmmV+fYyFDkiF9d/lI+Xr5zj5FAAGDRmlgAc1Tz6sidkT35RA+vifxS5NJqrz95ruGBjjwQGYq9UwcPlb/a0zMzrfumHm/vAQD0p4kF8NrIyyOXVHvd5cjd70ey/E3lE33aEcnRxTonR4Zi47q11bRvp9WrVrU2rFvT3gMA6E/TCuA1kQsj50bmG507LbIh8unI37bz85F8b27nx7t5MpLnCtTJq46HYv3kMdU5f1n6Uj5eueWM6nkAgMVoykUg+d+R5e+iyOsifx6Zz9GRF81uHnJFJEcG3xX5VuSpyEKGehFIynP+cto3R/6UPwBYPBeBNKcAXh/ZFnlz5Jv5RFuep1dfJXFj5MFInsfXze9GRvoqYABgcApgc6aAL4vkJ/G2yP6ObI3Ucr2U9bObAADlasoI4OFiBBAAxowRwOauAwgAQA8KIABAYRRAAIDCKIAAAIVRAAEACqMAAgAURgEEACiMAggAUBgFEACgMAogAEBhFEAAgMIogAAAhVEAAQAKowACABRGARxj+w880dp171T1CADQr1XtR5ZmInIgtCYmcnN4bt59f2v7zrtbz8xEi4/P4o4tZ7a2bjq1/VEAoJfp6enW5ORkbuYf07lRGiOAYyhH/Oryl/Lx8p17jAQCAH1RAMfQ3qmDh8pf7emZmdZ9U4+39wAAelMAx9DGdWurad9Oq1etam1Yt6a9BwDQmwI4htZPHlOd85elL+XjlVvOqJ4HAFiIi0AGc9guAkl5zl9O++bIn/IHAP1xEYgRwLGWpW/zaccrfwDAoiiAAACFUQABAAqjAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYRRAAIDCKIAAAIVRABsgbwm3696p6hEAYCHuBTyYw3ov4HTz7vtb23fe3XpmJtp8fDZ3bDmztXXTqe2PAgBzuRewEcCxliN+dflL+Xj5zj1GAgGAeSmAY2zv1MFD5a/29MxM676px9t7AADPpQCOsY3r1lbTvp1Wr1rV2rBuTXsPAOC5FMAxtn7ymOqcvyx9KR+v3HJG9TwAQC8uAhnMYb8IJOU5fzntmyN/yh8AzM9FIEYAGyFL3+bTjlf+AIC+KIAAAIVRAAEACqMANog7ggAA/XARyGBG4iKQ5I4gANAfF4EYAWwEdwQBABZDAWwAdwQBABZDAWwAdwQBABZDAWwAdwQBABbDRSCDGZmLQJI7ggDAwlwEYgSwUdwRBADohwIIAFAYBRAAoDAKYAO5IwgAMB8XgQxmpC4CSe4IAgDzcxGIEcBGcUcQAKAfCmCDuCMIANAPBbBB3BEEAOhHUwrg9sjuyGORhyO3RE6PzOfSyB2RR9v5fORVkbHljiAAQD+achHIZyM3RbIEHhH5QOTMyEsjByPdfDzyxciuyHcj/zyyJfKyyIORfozcRSDJHUEAoDcXgTT3KuATIjkSeE7k9nyiD6sjORL4q5Eb84k+jGQB7JRlMM8NzOlhZRAAFMDU1HMAq89qeKT92I88Ue4HIot5z0jLJWFec9UXWts++uXqMfcBAJpYAHNU8+rInZE9+USfrork1G+eC9jLUZEc6qtzbGQkWRIGAOiliQXw2sjLI5dUe/3J8//y9XkOYJ4P2EtebHKgIw9ERpIlYQCAXppWAK+JXBg5N9JvOfu1yOWRN0TuyifmsSOS08t1To6MJEvCAAC9NKUAZtXJkb8cwTsvsjfSj3dH3hf56chX8okFPBnJk0Xr5LIzI8mSMABAL025Cvj6yLbImyPfzCfacpq2Puktr+zNc/xyGjfltO/7I/m+XA6m9p12+jEWVwFbEgYAnuUq4OaMAF4WyU/ibZH9HdkaqZ0aWT+7WXlH5MjIf4l0vienhBsjS9/m046vtnfdO+UiEACgsesADsvIjwCmXP6lviI4zwvMqeGtm7IPA0B5jAA28ypgOlgOBgCYSwFsuF7LwXz1vrzpCQBQIgWw4botB5PeedPX3BkEAAqlADZcvRzM3E+0qWAAKJcCWIC84OPD217R3nuWO4MAQJkUwEK88gXHuTMIAFBRAAvhziAAQM06gIMZi3UAO7kzCAClsw6gEcDi1HcGyccsg+4OAgDlMQI4mLEbAay5OwgApTICaASwSDni5+4gAFAuBbBA7g4CAGVTAAvk7iAAUDYFsED1kjBzP/mmggGgDApgodwdBADKpQAWrNvdQfILYs2RviwAoMn8pC9YPRVc3x0kPRO56PpdzgUEgAZTAAuXU8E737G51dEBnQsIAA2nANI6+NTTrZkey8JkCXS3EABoli6LgbAIY3snkE5Z7l5z1ReeszZg/cWRT7tbCABN4U4gRgAJvZaFyeJXd0LTwgDQHAoglV7LwnRytxAAaAYFkEO6LQszl7uFAMD4UwA5ZO6yMN26oKlgABh/C4z3sIBGXAQyV5a7vBvIhnVrWl/9y0dbv/qJr7U/8qxPXnpWa/Npx7f3AGB8uAjECCBd5Ehglrt8nO9uIVkULREDAOPHCOBgGjkCOFee85fTvnkRSK3+wslnLBEDwDgxAmgEkD50u1tIFr+6DjovEADGiwJIX7rdLaRTjg5+/p6HDk0Jmx4GgNFlCngwRUwBpyxy3e4W0o3pYQBGmSlgI4D0qZ8lYmpZ/OqeaHoYAEaPEcDBFDMCWMsit9ASMd1YNgaAUWEE0Aggi7TQEjHd5BfZ4099rzon8Ov7HnVuIAAcZkYAB1PcCOBcPZeIiT/mu2jEuYEAHC5GABXAQRVfAFM9LZyLQz/+1DPV9PDD099tveX6XfOWwDyf8M5fP7caTQSAYVEAFcBBKYA95DTvto9+ub3X229c8JLWpg3HVcvMrD1y9fc9bly3VjkEYNkpgArgoBTAHhazbEwvpokBWAkKoItAWCFzl41ZCkvIAMDKMAI4GCOAC5h7fmA+7r7v0dYVn/lG+xUL6zZNnNPDae/UQVPHACyKEUAFcFAK4BIsx/Rw/YXb7R9h6hiA+SiApoA5DJZjejiLX6/+aOoYAOZnBHAwRgAH0G16eKnTxN24+wgA3RgBVAAHpQCukEGniXNo+w/+0atbJ04c3fU8QecSApRLAVQAB6UArqBudxmpVV+48cd8C03XX9zzdcj5XuNcQoBmUgAVwEEpgCus1zRx3m0kdX4s7zf8yzd+dd5SuFjuVgLQPAqgAjgoBXCE9Hv3kcW69pJXtF654biBppJ7fazzNUomwHAogArgoBTAEbIcy8t0M980cW2+1/TzftPNAMOjAFoGhgaZu7xM/rnQSjP9vCaL20Kdcr7X9PN+S9cAMEwL/OhjAUYAR1B93mC38wT7OZfwLx5+rPW+T/3v6vlhs3QNwMozAqgADkoBbKCVmkpeSA7Hf3jbK1qnHHdMz/MEU+e5hP2eX9jrNYO+f77X5HPOawRGkQKoAA5KAWyoziVoqv9J4o8Fl5zp8Zq+3x/m65zdXrPU99UGfX+t22uc1wiMKgXw2e/bLI0C2GCDTiX387F83PfIE6133vS1oY84DoNldIBRpAAqgINSABnYSi1fMyre/+aXtU478XmHbSp6Ma9ZiY8N8zXD/nso9owrBVABHJQCyMAO1zmHw9Ztmrg238dqw3jNSnysNozXDPPvYYqfcaYAWgYGDru5y9f0kh+d+5Juz80132sGfX+tn9dkmehVKOb7WG0Yr1mJj9WG8ZpB31/r5zWWLoLxtsC3bBZgBJBlU59z2O08wV7nEq70+YlLec3hXEaH4fuNC17S2rThuCVNIc99zXJOT/d6zTD+Hcv1mpX8d4zy363ba5b7dAMjgArgoBRAmKOUKW16q3+wzPcl0O01S33fXPO9Zhj/jtoo/11H+e8210qcbqAANmcKeHtkd+SxyMORWyKnRxbys5F7Ik+2Hy+KAAPoNaWde6MwFT3fa1biY7VhvGbQ99f6ec188of5fD/QU7fXLPV9c833mkHfXxvGawZ9f63ba5b6vtqg76/18xqnG6yMAf4XHymfjdwUyRJ4ROQDkTMjL40cjHSzOXJH5H2RP4hk+ftXkddG+r0k0wgg9DB3SnvYU9FLec1KfGyYrxnW32P3fY+2rvjMN6rXwLAs552SjAA2pwDOdUIkRwLPidyeT3RxcyRb2xurvVlZJB+NXFLtLUwBBIpjmp9hW+41RRXA5l4FXH1WwyPtx25yBPDW2c1DPhd59ewmAN0sdOV6P1PI3V6z1PfNNd9rBn1/bRivGfT9tW6vWer7aoO+v9bPa/Lr7MotZyz7hSClW+Cwj6X8b/pU5LjI2flED09F3h75RLU3a1vkY5Gjqr3nyuc7P3Zs5AEjgECJel25vtRp5mFMYQ/j37Fcr1nJf8co/926vcZVwMuviQXwusgFkTyX74F8oocsgL8Q+WS1N+utkRsiR1d7z/VbkX8xu/ksBRAAxocC2Lwp4GsiF0bOjcxX/tK3IyfNbh5yYuSh2c2udkTyi6XOyREAgLHSlAKYI5nXRrZEzovsjSzkS5HzZzcPeUNk1+xmV7lcTP6mUCeXnQEAGCtNKYA57fu2SJ7Dl6UsR/YynScN3BjJEbzahyJZ+N4TeXH78fWRD0YAABqrKQXwskhOyd4W2d+RrZFaLiG+fnazkiN9F0d+MXJXJC8Iydf3uwYgAMBYauJFIMNkHUAAGDMuAmnuOoAAAPSgAAIAFEYBBAAojAIIAFAYBRAAoDAKIABAYSwDM5hqGZh9+/ZZBgYAxkQuA3PKKafkZrHLwCiAg/nhyEL3HAYARlPe0//B2c2yKICDyeP3Q5GVuCfwsZEsl/nF6Z7DK8dxHh7Hejgc5+FxrIdnJY51/jP/KjJT7cGIyDnl/KI0t7yyHOfhcayHw3EeHsd6eBzrZeYiEACAwiiAAACFWd1+ZDQ9Hbmt/cjKcZyHx7EeDsd5eBzr4XGsAQAAAAAAAAAAAAAAgMPvHZG9ke9Gvho5O8LSbY/sjuTq8Q9HbomcHul0VOSayFTkYOQPI7niPIPJY5+Lt36w2pvlWC+fvB3l70f+X+TxyP+KvDJSy7sV/VYk73bwRCSvoHxZhP4dEbkikt+T8xj+38hvRjqXUXOcl+YnI5+O1HfjeEukUz/H9bjI70UOtJPbz4/A2NkaeSryy5GXRPIH53cip0ZYms9G3h7Jbxw/Gvlvkb+MrI3UPhLJ2wy9PvKKyBci+cPUUklLtymSPzS/HuksgI718sgffPdFPhZ5VWRD5Kcip0Vq74nkje63RM6I3BTJH6Z5Cyz6895I/rJyQSSP8c9F8pfJd0VqjvPSvDGS5TqPW7cC2M9x/aPI3ZHN7eR2lkoYO1+O5A/ITt+I7JjdZBmcEMlvNvnbZ5qMZOnO8l3LezznWlN/v9pjsZ4X+VYkS17+1l4XQMd6+VwVuWN2s6scPdkfyR+itRx9/ZvIr1R79CN/YbxhdvOQ/xrJkabkOC+PuQWwn+OagyT5vp+o9madFcnn5s7yMIc7gYyWIyM5fXNrtfes3H/17CbLIEtIeqT9mMf8ByKdxz1/y9wTcdyX5rrIZyKfr/ae5VgvnwsjX4n850ie2vC1yKWR2sbISZHOY/1k5E8jjnX/7ozkyOqPVHuzswivjfz3as9xXin9HNcc8ctp3xw4qf1ZJJ9z7BegAI6WdZGcBnuo2ntW7uf/CAwuf6u8OpLf1LN0pDy2OSr1aLX3LMd9aS6OZNHL8//mcqyXzwsjl0X+PJKjp78T+XDk5yOpPp6+nwzmX0c+Gfk/ke9FsmjniHY+lxznldHPcc3H/OVnrnzOsV+AAjiacvi6U5aWuc+xNNdGXh65pNqbn+O+eKdEPhR5ayQvYuqXY714+f37f0Yuj2Qp+feRj0ayFHby/WQwebrC2yLbIj8e+YXIr7UfOznOK2Oh49rtGDv2fVAAR0ueaJznQs39zeXEyNzfgli8vPI0p83OjeRFCLVvR3L6PU+q7+S4L16O/OVxy6vX/7adcyLvbG/n8XSsl0eeH3XP7OYheb5wfcFYfl0n308G828ieb5lXoCQFxjkuX+/HalHuB3nldHPcc3X/ODs5vfJ87wd+wUogKMlp8byB+f51d6zcn/X7CZLkL8N5shfXkl2XiSvTO2UxzyndjqP+/pIXnXmuC/On0TOjPxYR/I8tY93bDvWy+OLkbknuud5anmFe8qv8/wB2Xmss3xnIXes+7cm8szs5iH5i3r989NxXhn9HNcvRfKc7rwKvpYXhORzjj1jp14G5h9E8gqn/E0zl4F5QYSluT6SV47lN478bbLOMZFaXnm9L5Ine+fSJFlkLE2yPDqvAk6O9fLIZXayTOcU8IsiOUWZ6yrm9Hstr6DMr/2LIlmyPxGxPMni/G4kZwzqZWDyWP51JM8NrDnOS5OrBdS/KOaU7T9pb9ej2P0c11wGJpeayqt/M3dFLAPD2MqFoHN9r7ziKUen6uVKWJr8xtItuTZg7ehIThHXC+rmN5A8n43BzS2AjvXy+ZlITkvm+ZY5/dt5FXDK0e9cSDeni/M1eQVl/iClf1k28us3R1ZzMeJ7I7l2XY5G1RznpXldpNv35izdqZ/j+ncjuRh6rheYyW0LQQMAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADQEK3W/wcmtTzXrmmDzgAAAABJRU5ErkJggg==\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots()\n",
    "ax.plot(curve_x_data[\"Temperature (K)\"], curve_x_data['log Ohm'], '.')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Go to a set point"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [],
   "source": [
    "import time\n",
    "\n",
    "import numpy\n",
    "from IPython.display import display\n",
    "from ipywidgets import widgets\n",
    "from matplotlib import pyplot as plt\n",
    "\n",
    "\n",
    "def live_plot_temperature_reading(channel_to_read, read_period=0.2, n_reads=1000):\n",
    "    \"\"\"\n",
    "    Live plot the temperature reading from a Lakeshore sensor channel\n",
    "\n",
    "    Args:\n",
    "        channel_to_read:\n",
    "            Lakeshore channel object to read the temperature from\n",
    "        read_period:\n",
    "            time in seconds between two reads of the temperature\n",
    "        n_reads:\n",
    "            total number of reads to perform\n",
    "    \"\"\"\n",
    "\n",
    "    # Make a widget for a text display that is contantly being updated\n",
    "    text = widgets.Text()\n",
    "    display(text)\n",
    "\n",
    "    fig, ax = plt.subplots(1)\n",
    "    line, = ax.plot([], [], '*-')\n",
    "    ax.set_xlabel('Time, s')\n",
    "    ax.set_ylabel(f'Temperature, {channel_to_read.temperature.unit}')\n",
    "    fig.show()\n",
    "    plt.ion()\n",
    "\n",
    "    for i in range(n_reads):\n",
    "        time.sleep(read_period)\n",
    "\n",
    "        # Update the text field\n",
    "        text.value = f'T = {channel_to_read.temperature()}'\n",
    "\n",
    "        # Add new point to the data that is being plotted\n",
    "        line.set_ydata(numpy.append(line.get_ydata(), channel_to_read.temperature()))\n",
    "        line.set_xdata(numpy.arange(0, len(line.get_ydata()), 1)*read_period)\n",
    "\n",
    "        ax.relim()  # Recalculate limits\n",
    "        ax.autoscale_view(True, True, True)  # Autoscale\n",
    "        fig.canvas.draw()  # Redraw"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "5bf837db2f5b4c9d82489aaeaf839847",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "Text(value='')"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/javascript": "/* Put everything inside the global mpl namespace */\nwindow.mpl = {};\n\n\nmpl.get_websocket_type = function() {\n    if (typeof(WebSocket) !== 'undefined') {\n        return WebSocket;\n    } else if (typeof(MozWebSocket) !== 'undefined') {\n        return MozWebSocket;\n    } else {\n        alert('Your browser does not have WebSocket support.' +\n              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n              'Firefox 4 and 5 are also supported but you ' +\n              'have to enable WebSockets in about:config.');\n    };\n}\n\nmpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n    this.id = figure_id;\n\n    this.ws = websocket;\n\n    this.supports_binary = (this.ws.binaryType != undefined);\n\n    if (!this.supports_binary) {\n        var warnings = document.getElementById(\"mpl-warnings\");\n        if (warnings) {\n            warnings.style.display = 'block';\n            warnings.textContent = (\n                \"This browser does not support binary websocket messages. \" +\n                    \"Performance may be slow.\");\n        }\n    }\n\n    this.imageObj = new Image();\n\n    this.context = undefined;\n    this.message = undefined;\n    this.canvas = undefined;\n    this.rubberband_canvas = undefined;\n    this.rubberband_context = undefined;\n    this.format_dropdown = undefined;\n\n    this.image_mode = 'full';\n\n    this.root = $('<div/>');\n    this._root_extra_style(this.root)\n    this.root.attr('style', 'display: inline-block');\n\n    $(parent_element).append(this.root);\n\n    this._init_header(this);\n    this._init_canvas(this);\n    this._init_toolbar(this);\n\n    var fig = this;\n\n    this.waiting = false;\n\n    this.ws.onopen =  function () {\n            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n            fig.send_message(\"send_image_mode\", {});\n            if (mpl.ratio != 1) {\n                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n            }\n            fig.send_message(\"refresh\", {});\n        }\n\n    this.imageObj.onload = function() {\n            if (fig.image_mode == 'full') {\n                // Full images could contain transparency (where diff images\n                // almost always do), so we need to clear the canvas so that\n                // there is no ghosting.\n                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n            }\n            fig.context.drawImage(fig.imageObj, 0, 0);\n        };\n\n    this.imageObj.onunload = function() {\n        this.ws.close();\n    }\n\n    this.ws.onmessage = this._make_on_message_function(this);\n\n    this.ondownload = ondownload;\n}\n\nmpl.figure.prototype._init_header = function() {\n    var titlebar = $(\n        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n        'ui-helper-clearfix\"/>');\n    var titletext = $(\n        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n        'text-align: center; padding: 3px;\"/>');\n    titlebar.append(titletext)\n    this.root.append(titlebar);\n    this.header = titletext[0];\n}\n\n\n\nmpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n\n}\n\n\nmpl.figure.prototype._root_extra_style = function(canvas_div) {\n\n}\n\nmpl.figure.prototype._init_canvas = function() {\n    var fig = this;\n\n    var canvas_div = $('<div/>');\n\n    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n\n    function canvas_keyboard_event(event) {\n        return fig.key_event(event, event['data']);\n    }\n\n    canvas_div.keydown('key_press', canvas_keyboard_event);\n    canvas_div.keyup('key_release', canvas_keyboard_event);\n    this.canvas_div = canvas_div\n    this._canvas_extra_style(canvas_div)\n    this.root.append(canvas_div);\n\n    var canvas = $('<canvas/>');\n    canvas.addClass('mpl-canvas');\n    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n\n    this.canvas = canvas[0];\n    this.context = canvas[0].getContext(\"2d\");\n\n    var backingStore = this.context.backingStorePixelRatio ||\n\tthis.context.webkitBackingStorePixelRatio ||\n\tthis.context.mozBackingStorePixelRatio ||\n\tthis.context.msBackingStorePixelRatio ||\n\tthis.context.oBackingStorePixelRatio ||\n\tthis.context.backingStorePixelRatio || 1;\n\n    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n\n    var rubberband = $('<canvas/>');\n    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n\n    var pass_mouse_events = true;\n\n    canvas_div.resizable({\n        start: function(event, ui) {\n            pass_mouse_events = false;\n        },\n        resize: function(event, ui) {\n            fig.request_resize(ui.size.width, ui.size.height);\n        },\n        stop: function(event, ui) {\n            pass_mouse_events = true;\n            fig.request_resize(ui.size.width, ui.size.height);\n        },\n    });\n\n    function mouse_event_fn(event) {\n        if (pass_mouse_events)\n            return fig.mouse_event(event, event['data']);\n    }\n\n    rubberband.mousedown('button_press', mouse_event_fn);\n    rubberband.mouseup('button_release', mouse_event_fn);\n    // Throttle sequential mouse events to 1 every 20ms.\n    rubberband.mousemove('motion_notify', mouse_event_fn);\n\n    rubberband.mouseenter('figure_enter', mouse_event_fn);\n    rubberband.mouseleave('figure_leave', mouse_event_fn);\n\n    canvas_div.on(\"wheel\", function (event) {\n        event = event.originalEvent;\n        event['data'] = 'scroll'\n        if (event.deltaY < 0) {\n            event.step = 1;\n        } else {\n            event.step = -1;\n        }\n        mouse_event_fn(event);\n    });\n\n    canvas_div.append(canvas);\n    canvas_div.append(rubberband);\n\n    this.rubberband = rubberband;\n    this.rubberband_canvas = rubberband[0];\n    this.rubberband_context = rubberband[0].getContext(\"2d\");\n    this.rubberband_context.strokeStyle = \"#000000\";\n\n    this._resize_canvas = function(width, height) {\n        // Keep the size of the canvas, canvas container, and rubber band\n        // canvas in synch.\n        canvas_div.css('width', width)\n        canvas_div.css('height', height)\n\n        canvas.attr('width', width * mpl.ratio);\n        canvas.attr('height', height * mpl.ratio);\n        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n\n        rubberband.attr('width', width);\n        rubberband.attr('height', height);\n    }\n\n    // Set the figure to an initial 600x600px, this will subsequently be updated\n    // upon first draw.\n    this._resize_canvas(600, 600);\n\n    // Disable right mouse context menu.\n    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n        return false;\n    });\n\n    function set_focus () {\n        canvas.focus();\n        canvas_div.focus();\n    }\n\n    window.setTimeout(set_focus, 100);\n}\n\nmpl.figure.prototype._init_toolbar = function() {\n    var fig = this;\n\n    var nav_element = $('<div/>')\n    nav_element.attr('style', 'width: 100%');\n    this.root.append(nav_element);\n\n    // Define a callback function for later on.\n    function toolbar_event(event) {\n        return fig.toolbar_button_onclick(event['data']);\n    }\n    function toolbar_mouse_event(event) {\n        return fig.toolbar_button_onmouseover(event['data']);\n    }\n\n    for(var toolbar_ind in mpl.toolbar_items) {\n        var name = mpl.toolbar_items[toolbar_ind][0];\n        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n        var image = mpl.toolbar_items[toolbar_ind][2];\n        var method_name = mpl.toolbar_items[toolbar_ind][3];\n\n        if (!name) {\n            // put a spacer in here.\n            continue;\n        }\n        var button = $('<button/>');\n        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n                        'ui-button-icon-only');\n        button.attr('role', 'button');\n        button.attr('aria-disabled', 'false');\n        button.click(method_name, toolbar_event);\n        button.mouseover(tooltip, toolbar_mouse_event);\n\n        var icon_img = $('<span/>');\n        icon_img.addClass('ui-button-icon-primary ui-icon');\n        icon_img.addClass(image);\n        icon_img.addClass('ui-corner-all');\n\n        var tooltip_span = $('<span/>');\n        tooltip_span.addClass('ui-button-text');\n        tooltip_span.html(tooltip);\n\n        button.append(icon_img);\n        button.append(tooltip_span);\n\n        nav_element.append(button);\n    }\n\n    var fmt_picker_span = $('<span/>');\n\n    var fmt_picker = $('<select/>');\n    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n    fmt_picker_span.append(fmt_picker);\n    nav_element.append(fmt_picker_span);\n    this.format_dropdown = fmt_picker[0];\n\n    for (var ind in mpl.extensions) {\n        var fmt = mpl.extensions[ind];\n        var option = $(\n            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n        fmt_picker.append(option)\n    }\n\n    // Add hover states to the ui-buttons\n    $( \".ui-button\" ).hover(\n        function() { $(this).addClass(\"ui-state-hover\");},\n        function() { $(this).removeClass(\"ui-state-hover\");}\n    );\n\n    var status_bar = $('<span class=\"mpl-message\"/>');\n    nav_element.append(status_bar);\n    this.message = status_bar[0];\n}\n\nmpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n    // which will in turn request a refresh of the image.\n    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n}\n\nmpl.figure.prototype.send_message = function(type, properties) {\n    properties['type'] = type;\n    properties['figure_id'] = this.id;\n    this.ws.send(JSON.stringify(properties));\n}\n\nmpl.figure.prototype.send_draw_message = function() {\n    if (!this.waiting) {\n        this.waiting = true;\n        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n    }\n}\n\n\nmpl.figure.prototype.handle_save = function(fig, msg) {\n    var format_dropdown = fig.format_dropdown;\n    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n    fig.ondownload(fig, format);\n}\n\n\nmpl.figure.prototype.handle_resize = function(fig, msg) {\n    var size = msg['size'];\n    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n        fig._resize_canvas(size[0], size[1]);\n        fig.send_message(\"refresh\", {});\n    };\n}\n\nmpl.figure.prototype.handle_rubberband = function(fig, msg) {\n    var x0 = msg['x0'] / mpl.ratio;\n    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n    var x1 = msg['x1'] / mpl.ratio;\n    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n    x0 = Math.floor(x0) + 0.5;\n    y0 = Math.floor(y0) + 0.5;\n    x1 = Math.floor(x1) + 0.5;\n    y1 = Math.floor(y1) + 0.5;\n    var min_x = Math.min(x0, x1);\n    var min_y = Math.min(y0, y1);\n    var width = Math.abs(x1 - x0);\n    var height = Math.abs(y1 - y0);\n\n    fig.rubberband_context.clearRect(\n        0, 0, fig.canvas.width, fig.canvas.height);\n\n    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n}\n\nmpl.figure.prototype.handle_figure_label = function(fig, msg) {\n    // Updates the figure title.\n    fig.header.textContent = msg['label'];\n}\n\nmpl.figure.prototype.handle_cursor = function(fig, msg) {\n    var cursor = msg['cursor'];\n    switch(cursor)\n    {\n    case 0:\n        cursor = 'pointer';\n        break;\n    case 1:\n        cursor = 'default';\n        break;\n    case 2:\n        cursor = 'crosshair';\n        break;\n    case 3:\n        cursor = 'move';\n        break;\n    }\n    fig.rubberband_canvas.style.cursor = cursor;\n}\n\nmpl.figure.prototype.handle_message = function(fig, msg) {\n    fig.message.textContent = msg['message'];\n}\n\nmpl.figure.prototype.handle_draw = function(fig, msg) {\n    // Request the server to send over a new figure.\n    fig.send_draw_message();\n}\n\nmpl.figure.prototype.handle_image_mode = function(fig, msg) {\n    fig.image_mode = msg['mode'];\n}\n\nmpl.figure.prototype.updated_canvas_event = function() {\n    // Called whenever the canvas gets updated.\n    this.send_message(\"ack\", {});\n}\n\n// A function to construct a web socket function for onmessage handling.\n// Called in the figure constructor.\nmpl.figure.prototype._make_on_message_function = function(fig) {\n    return function socket_on_message(evt) {\n        if (evt.data instanceof Blob) {\n            /* FIXME: We get \"Resource interpreted as Image but\n             * transferred with MIME type text/plain:\" errors on\n             * Chrome.  But how to set the MIME type?  It doesn't seem\n             * to be part of the websocket stream */\n            evt.data.type = \"image/png\";\n\n            /* Free the memory for the previous frames */\n            if (fig.imageObj.src) {\n                (window.URL || window.webkitURL).revokeObjectURL(\n                    fig.imageObj.src);\n            }\n\n            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n                evt.data);\n            fig.updated_canvas_event();\n            fig.waiting = false;\n            return;\n        }\n        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n            fig.imageObj.src = evt.data;\n            fig.updated_canvas_event();\n            fig.waiting = false;\n            return;\n        }\n\n        var msg = JSON.parse(evt.data);\n        var msg_type = msg['type'];\n\n        // Call the  \"handle_{type}\" callback, which takes\n        // the figure and JSON message as its only arguments.\n        try {\n            var callback = fig[\"handle_\" + msg_type];\n        } catch (e) {\n            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n            return;\n        }\n\n        if (callback) {\n            try {\n                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n                callback(fig, msg);\n            } catch (e) {\n                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n            }\n        }\n    };\n}\n\n// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\nmpl.findpos = function(e) {\n    //this section is from http://www.quirksmode.org/js/events_properties.html\n    var targ;\n    if (!e)\n        e = window.event;\n    if (e.target)\n        targ = e.target;\n    else if (e.srcElement)\n        targ = e.srcElement;\n    if (targ.nodeType == 3) // defeat Safari bug\n        targ = targ.parentNode;\n\n    // jQuery normalizes the pageX and pageY\n    // pageX,Y are the mouse positions relative to the document\n    // offset() returns the position of the element relative to the document\n    var x = e.pageX - $(targ).offset().left;\n    var y = e.pageY - $(targ).offset().top;\n\n    return {\"x\": x, \"y\": y};\n};\n\n/*\n * return a copy of an object with only non-object keys\n * we need this to avoid circular references\n * http://stackoverflow.com/a/24161582/3208463\n */\nfunction simpleKeys (original) {\n  return Object.keys(original).reduce(function (obj, key) {\n    if (typeof original[key] !== 'object')\n        obj[key] = original[key]\n    return obj;\n  }, {});\n}\n\nmpl.figure.prototype.mouse_event = function(event, name) {\n    var canvas_pos = mpl.findpos(event)\n\n    if (name === 'button_press')\n    {\n        this.canvas.focus();\n        this.canvas_div.focus();\n    }\n\n    var x = canvas_pos.x * mpl.ratio;\n    var y = canvas_pos.y * mpl.ratio;\n\n    this.send_message(name, {x: x, y: y, button: event.button,\n                             step: event.step,\n                             guiEvent: simpleKeys(event)});\n\n    /* This prevents the web browser from automatically changing to\n     * the text insertion cursor when the button is pressed.  We want\n     * to control all of the cursor setting manually through the\n     * 'cursor' event from matplotlib */\n    event.preventDefault();\n    return false;\n}\n\nmpl.figure.prototype._key_event_extra = function(event, name) {\n    // Handle any extra behaviour associated with a key event\n}\n\nmpl.figure.prototype.key_event = function(event, name) {\n\n    // Prevent repeat events\n    if (name == 'key_press')\n    {\n        if (event.which === this._key)\n            return;\n        else\n            this._key = event.which;\n    }\n    if (name == 'key_release')\n        this._key = null;\n\n    var value = '';\n    if (event.ctrlKey && event.which != 17)\n        value += \"ctrl+\";\n    if (event.altKey && event.which != 18)\n        value += \"alt+\";\n    if (event.shiftKey && event.which != 16)\n        value += \"shift+\";\n\n    value += 'k';\n    value += event.which.toString();\n\n    this._key_event_extra(event, name);\n\n    this.send_message(name, {key: value,\n                             guiEvent: simpleKeys(event)});\n    return false;\n}\n\nmpl.figure.prototype.toolbar_button_onclick = function(name) {\n    if (name == 'download') {\n        this.handle_save(this, null);\n    } else {\n        this.send_message(\"toolbar_button\", {name: name});\n    }\n};\n\nmpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n    this.message.textContent = tooltip;\n};\nmpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n\nmpl.extensions = [\"eps\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\"];\n\nmpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n    // Create a \"websocket\"-like object which calls the given IPython comm\n    // object with the appropriate methods. Currently this is a non binary\n    // socket, so there is still some room for performance tuning.\n    var ws = {};\n\n    ws.close = function() {\n        comm.close()\n    };\n    ws.send = function(m) {\n        //console.log('sending', m);\n        comm.send(m);\n    };\n    // Register the callback with on_msg.\n    comm.on_msg(function(msg) {\n        //console.log('receiving', msg['content']['data'], msg);\n        // Pass the mpl event to the overriden (by mpl) onmessage function.\n        ws.onmessage(msg['content']['data'])\n    });\n    return ws;\n}\n\nmpl.mpl_figure_comm = function(comm, msg) {\n    // This is the function which gets called when the mpl process\n    // starts-up an IPython Comm through the \"matplotlib\" channel.\n\n    var id = msg.content.data.id;\n    // Get hold of the div created by the display call when the Comm\n    // socket was opened in Python.\n    var element = $(\"#\" + id);\n    var ws_proxy = comm_websocket_adapter(comm)\n\n    function ondownload(figure, format) {\n        window.open(figure.imageObj.src);\n    }\n\n    var fig = new mpl.figure(id, ws_proxy,\n                           ondownload,\n                           element.get(0));\n\n    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n    // web socket which is closed, not our websocket->open comm proxy.\n    ws_proxy.onopen();\n\n    fig.parent_element = element.get(0);\n    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n    if (!fig.cell_info) {\n        console.error(\"Failed to find cell for figure\", id, fig);\n        return;\n    }\n\n    var output_index = fig.cell_info[2]\n    var cell = fig.cell_info[0];\n\n};\n\nmpl.figure.prototype.handle_close = function(fig, msg) {\n    var width = fig.canvas.width/mpl.ratio\n    fig.root.unbind('remove')\n\n    // Update the output cell to use the data from the current canvas.\n    fig.push_to_output();\n    var dataURL = fig.canvas.toDataURL();\n    // Re-enable the keyboard manager in IPython - without this line, in FF,\n    // the notebook keyboard shortcuts fail.\n    IPython.keyboard_manager.enable()\n    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n    fig.close_ws(fig, msg);\n}\n\nmpl.figure.prototype.close_ws = function(fig, msg){\n    fig.send_message('closing', msg);\n    // fig.ws.close()\n}\n\nmpl.figure.prototype.push_to_output = function(remove_interactive) {\n    // Turn the data on the canvas into data in the output cell.\n    var width = this.canvas.width/mpl.ratio\n    var dataURL = this.canvas.toDataURL();\n    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n}\n\nmpl.figure.prototype.updated_canvas_event = function() {\n    // Tell IPython that the notebook contents must change.\n    IPython.notebook.set_dirty(true);\n    this.send_message(\"ack\", {});\n    var fig = this;\n    // Wait a second, then push the new image to the DOM so\n    // that it is saved nicely (might be nice to debounce this).\n    setTimeout(function () { fig.push_to_output() }, 1000);\n}\n\nmpl.figure.prototype._init_toolbar = function() {\n    var fig = this;\n\n    var nav_element = $('<div/>')\n    nav_element.attr('style', 'width: 100%');\n    this.root.append(nav_element);\n\n    // Define a callback function for later on.\n    function toolbar_event(event) {\n        return fig.toolbar_button_onclick(event['data']);\n    }\n    function toolbar_mouse_event(event) {\n        return fig.toolbar_button_onmouseover(event['data']);\n    }\n\n    for(var toolbar_ind in mpl.toolbar_items){\n        var name = mpl.toolbar_items[toolbar_ind][0];\n        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n        var image = mpl.toolbar_items[toolbar_ind][2];\n        var method_name = mpl.toolbar_items[toolbar_ind][3];\n\n        if (!name) { continue; };\n\n        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n        button.click(method_name, toolbar_event);\n        button.mouseover(tooltip, toolbar_mouse_event);\n        nav_element.append(button);\n    }\n\n    // Add the status bar.\n    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n    nav_element.append(status_bar);\n    this.message = status_bar[0];\n\n    // Add the close button to the window.\n    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n    button.click(function (evt) { fig.handle_close(fig, {}); } );\n    button.mouseover('Stop Interaction', toolbar_mouse_event);\n    buttongrp.append(button);\n    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n    titlebar.prepend(buttongrp);\n}\n\nmpl.figure.prototype._root_extra_style = function(el){\n    var fig = this\n    el.on(\"remove\", function(){\n\tfig.close_ws(fig, {});\n    });\n}\n\nmpl.figure.prototype._canvas_extra_style = function(el){\n    // this is important to make the div 'focusable\n    el.attr('tabindex', 0)\n    // reach out to IPython and tell the keyboard manager to turn it's self\n    // off when our div gets focus\n\n    // location in version 3\n    if (IPython.notebook.keyboard_manager) {\n        IPython.notebook.keyboard_manager.register_events(el);\n    }\n    else {\n        // location in version 2\n        IPython.keyboard_manager.register_events(el);\n    }\n\n}\n\nmpl.figure.prototype._key_event_extra = function(event, name) {\n    var manager = IPython.notebook.keyboard_manager;\n    if (!manager)\n        manager = IPython.keyboard_manager;\n\n    // Check for shift+enter\n    if (event.shiftKey && event.which == 13) {\n        this.canvas_div.blur();\n        // select the cell after this one\n        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n        IPython.notebook.select(index + 1);\n    }\n}\n\nmpl.figure.prototype.handle_save = function(fig, msg) {\n    fig.ondownload(fig, null);\n}\n\n\nmpl.find_output_cell = function(html_output) {\n    // Return the cell and output element which can be found *uniquely* in the notebook.\n    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n    // IPython event is triggered only after the cells have been serialised, which for\n    // our purposes (turning an active figure into a static one), is too late.\n    var cells = IPython.notebook.get_cells();\n    var ncells = cells.length;\n    for (var i=0; i<ncells; i++) {\n        var cell = cells[i];\n        if (cell.cell_type === 'code'){\n            for (var j=0; j<cell.output_area.outputs.length; j++) {\n                var data = cell.output_area.outputs[j];\n                if (data.data) {\n                    // IPython >= 3 moved mimebundle to data attribute of output\n                    data = data.data;\n                }\n                if (data['text/html'] == html_output) {\n                    return [cell, data, j];\n                }\n            }\n        }\n    }\n}\n\n// Register the function which deals with the matplotlib target/channel.\n// The kernel may be null if the page has been refreshed.\nif (IPython.notebook.kernel != null) {\n    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n}\n",
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAADpHSURBVHhe7d0JnBx1mf/xX+bMTCYXIZBAEsKREEBBJKAYrgABwyFJDILLiiCwK4KiZF2CwsJ6LCigIIcoh4Cgfy8EkYVdEVASEpFDCIoLiCSZEA4TMrkmk2Qy/+dbXQWdnu6e7qmuma6qz/v1+jpVPUmGdq6nn9/lAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADJMsB/i97R/387WNZ4dwAAIC4GW16zdHl3KUMBGM6OltbMJQAAiJkxlmWZy3ShAAxniKVt6dKlbsgQXQIAgGq3evVqN3bsWF0O1a0u0oYCMByvADQUgAAAxIQKwKFDVfultwCs8d8CAAAgJSgAAQAAUoYCEAAAIGUoAAEAAFKGAhAAACBlKAABAABShgIQAAAgZSgAAQAAUoYCEAAAIGUoAKvcc62r3Me/v9B7CwAAUAlxKQAPsdxnec3SZZlhyXabRY9nZ6GlmNMsuX9HGWipGnc/vcwteGWF+5c7nqIIBAAAFRGXAnCQ5VnLud5dfg9aRmflGEtPdP5f9t9RNlj6Vevb663wa3VHXvU798OFi73HXl+9wZ1+2x/d3U+1eu8HAADorQH+2zhRl26m5R7vLkMdwGGW3M5gMeoAXm3R3+utIZY244YM0WVljJ97v39V2LwLproxw5v9OwAAUKrVq1e7oUOH6lL/o2ZQ6iRpDuBhljctL1pusmxn6UmLRS22VsuvLfta+o06e7/9yxvuuL3ViCzuoG88QicQAAD0SlI6gCdZ1lpUzO1s+aqlzrKfpcOSzwctu1kWWdS+O8+iYeN9LC9Z8mn0Exhsaa1UB7CUzl8uOoEAAJSHDmByOoA/sah6et6ixSLTLRMtx1oK0SKROy2aW/iY5WMWdQ8/aynkQktbVtQ5rJirT3qfqymzJKcTCAAAypXUbWCWW9QNnODdlWaL5Y+WYn/nMoteLQQZY6mYGfvu6H517kH+XelUBAIAAJQqqQXgCMtYiwrBUqn39j5Lsb+j4WS1ioOssfS7Uw4YRxcQAACULC4FoBZrqDhTRPP8dD3OovddaTnQMt6ixSAaBv6H5ZeWwB0WdfACl1iOtuxi0b91i//2Rku/GdHS4LZpbnB1NQPcjsOa3PaDG/z3FHbXE0voAgIAgJLFZRGIirp8Fc7tlrMtWhCiFbza0kUdPP3Ziy1LLYFHLa9atP2LfNsyyzLKovl8z1gutSywlCqSbWA6Nne6rq4u11hX692/9OYaN+O6+W79Jo1S56cu4MmW4YPqWRQCAEARLAKJ5yrgahJJAZjP4hXr3AlWBK5q3+Q/Utirlxdb+wIAQLpRACZ3DmDi7DRikLv19Mn+XX6q5udMm8h8QAAAUBQFYIyMHtrkRrZkb0O4NW2QeNVvXmRrGAAAUBQFYIyoAJw3d6qbM63n3W1YFAIAAAqhAIwZLQyZPXmsa67PLBDJRyuItak0AABAPhSAMaRO4J1nHeDfdXflifu4XUe2MAwMAADyYhVwOH22CjjX8rZ2d+BlD/t3hbEiGACArbEKmA5gbKkLeNGxe/h33TEMDAAACqEAjLEzD97F3fiJ9/t3W/vcERMYBgYAAHkxBBxOvw0BB8bPvd+/KoxhYAAA3sUQMB3A2Lv0+D39q+5qrbxnY2gAAJCLAjDmTpuys3/VXWfXuxtDAwAABCgAE2DOURP9q+5YDAIAAHJRACbA7P3GuCGNdf7d1tgTEAAA5GIRSDj9vggkwGIQAABKwyIQOoCJUWwxCMPAAAAgGwVgQmgxSKE9ARkGBgAA2RgCDqdqhoCFYWAAAHrGEDAdwEQpNgysSp9hYAAAIBSACaJh4Ls/c6B/t7WamgGuq6vLLWptYygYAICUowBMmIbaWv9qa51butwXfvqsO/66eWwMDQBAylEAJsyIlgY3uMCegMLxcAAAgEUg4VTVIpBAx+ZOt/tFD/p3hbEgBACQRiwCoQOYSI11tRwPBwAACqIATCiOhwMAAIUwBBxOVQ4BB9gXEACA7hgCpgOYaBwPBwAA8qEATDCOhwMAAPkwBBxOVQ8BC8PAAABsjSFgOoCJxzAwAADIRQGYcAwDAwCAXAwBh1P1Q8DCMDAAAO9iCJgOYCowDAwAALJRAKYAw8AAACAbQ8DhxGIIWBgGBgAggyFgOoCpwTAwAAAIUACmRLFh4J23HeS2bOliGBgAgJRgCDic2AwBC8PAAAAwBCx0AFOk2DCwnHLAOLeotY1OIAAACUcHMJxYdQDl6SUr3awbFvh3hdEJBAAkFR1AOoCp01Bb61/lV2svCeZMm0gXEACABKMATJkRLQ1uZEujf9ddZ5dzV/3mRXfQNx7xHwEAAElDAZgyo4c2uXlzp7o50yb4j3THtjAAACQbBWAKNdbVutmTx7phTfX+I1u755wpbsa+O/p3AAAgaSgAU0qdwMtnv9e/29rLb65lNTAAAAnGKuBwYrcKOBv7AgIA0ohVwHQAU03z/Gpr8r8GYB4gAADJRQGYYprnd9Mn9vPvtnblifu4XUe2MAwMAEACMQQcTqyHgIVhYABA2jAETAcw9YodD6dNoRkGBgAgeSgAU+60KTu7Gz/xfv9ua1N225ZhYAAAEogh4HBiPwQspQwDz7tgqhszvNm/AwAgvhgCpgMIU2wYOMDRcAAAJAcFILxh4POO3M2/644tYQAASBYKQHim7THKv+qOLWEAAEiWuBSAh1jus7xm6bLMsGS7zaLHs7PQUqqTLfo793h3KTSipcG/6u7zP/mTO/66eQwDAwCQEHEpAAdZnrWc693l96BldFaOsZRiJ8uVlse8u5TS2cAXHbuHf9cdw8AAACRHXArABywXWe727vLrsLyelZWWntRa7rJcYnlFD6TZmQfvUnBLmJ23HeS2bOliGBgAgARI0hzAwyxvWl603GTZztKT/7C8ZbnFu4P79A+f9q+29tKba935P3uWYWAAABIgKQWgOoSnWA63zLHsb3nY0mgpZIrlDMtZ3l1p9O9pw78ggy2J0tOWMHOmTXSLWtvoBAIAEGNx3AhaizVmWoot2NAcwMUWLe7IN2yswu05y2csKh5FC0mGWXIXmGS71KLh4q3EfSPoXE8vWelm3bDAvyuMM4IBAHHERtDJGgLOttyiAnCCd9fdrpbxFq0s3uznVMtH/Gu9P5/LLPpiCTLGkjgNtZoaWRgLQgAAiLekFoAjLGMtKgTz+avlvRZVMUF+ZdEEN10vteSjhSZ6pRBkjSVxtCXMyJZGt9OI/Ee/sS8gAADxFpch4BZLcFTFM5bzLSrWtNJX0dDsLywq+NTZ+y/LOIv2NQmKtDssyywXenfdlTIEnEvjvrE/Czifjs2dbveLtLNOcQwDAwDihiHg+HQAJ1tU+CnyLYuuv2LptKibd69FK4Bv998eaMnu0Kkg1NxAlKCxrrboghCGgQEAiK84LgKpJontAAYe/PPyvFvDqPjTMPDwQfVuzPD8Q8UAAFQjOoAUgGElvgAcP/d+/6owhoEBAHFCAZjcRSCoEIaBAQBIHgpAFHXalJ0LHg/3uSMmsBoYAIAYYgg4nMQPAQvDwACAJGEImA4gSlBsGFhfQDoeji4gAADxQQGIHmkYuJAtlqt+86I76BvalhEAAMQBBSBKMueoif5Vd3QBAQCIFwpAlGT2foWPPaYLCABAvFAAoiSjhza5844MTuPrrnYAXUAAAOKCAhAlO3l/naaXX2cXXUAAAOKCAhAl66kLyMbQAADEAwUgyqIu4LCmev9ua/ecM8XN2HdH/w4AAFQrCkCURV3Ay2e/17/b2hd+8id391OtzAMEAKDKcRJIOKk4CSQXJ4MAAOKMk0DoAKIXNM+vtqbwa4d9xw5zi1rb6AQCAFCl6ACGk8oOoDy/rM0dd+08/66weRdMdWOGN/t3AAD0PzqAdAARMbaFAQCg+lAAoldGtDS4kS2Nbvshjf4j3bEtDAAA1YkCEL2i1cDz5k51N5862X+ku88dMcHtOrKFuYAAAFQZ5gCGk9o5gIHlbe3uwMse9u8KY1UwAKBaMAeQDiBCUifwomP38O+644xgAACqDwUgQjvz4F38q+44IxgAgOpDAYiKmHPURP+qOxaDAABQXSgAURGz9xvjhjTW+XdbYzEIAADVhUUg4aR+EUg2jogDAMQBi0DoAKKCLj1+T/+qO32hsRgEAIDqQAGIijltys7+VXdbLCwGAQCgOlAAoqKKLQaRUw4YRxcQAIB+RgGIitJikGLuemIJXUAAAPoZBSAqShtDn3fkbv5dfuoCLmptoxMIAEA/YRVwOKwCzkPHw02/+jG3qn2T/0hhrAoGAPQ1VgHTAUQE1AW89fTJ/l1+HBEHAED/oQBEJFQEjmxp9O+644g4AAD6DwUgIqECcN7cqW7OtAn+I91xRBwAAP2DAhCRaayrdbMnj+WIOAAAqgyLQMJhEUgJOCIOAFBNWARCBxB9gCPiAACoLhSAiBxHxAEAUF0oANEnih0RRxcQAIC+RQGIPlHsiDi6gAAA9C0KQPSJno6IowsIAEDfoQBEnzl5/3H+VXd0AQEA6DtRFIC7+2+L+bz/FilCFxAAgOoQRQH4v5bCE76cO8/yzcwl0oYuIAAA/S+KAvApy28sI7y7rX3WcqXlU94dUqenLqCccsA4uoAAAEQoigLwZMsyy4OWFj3gO8fyLcsZljv1ANKpWBdQ7npiCV1AAAAiFEUBuNFygqXL8itLg+VsyzWWf7HcYUGKldoFXNTaRicQAIAIRHkW8DaW31s2W/ayfNpyiyVJOAu4l5a3tbvpVz/mVrVv8h8pjHOCAQCVxFnA0XQAj/Iz2aIh30mWuy1LLcH7FKSYuoC3nq4vkcJq7eUJq4IBAKi8KDqAWszZEw0P12YuY40OYAjqAn7k2vnurbUd/iOF0QUEAFQKHcBoOoBNJaTZgpRTF3De3KluzrQJ/iPdqQt49Unv8+8AAEAlRFEAqp1TSgDXWFfrZk8e64Y01vmPbK1lYL3bsqWLYWAAACooigIQKIs6gas7tFaou7b2Te78nz3rbQtDEQgAQGVQAKIqXHr8nv5VYRSBAABUBgUgqsJpU3bucW9AYYNoAADCi0sBeIjlPstrFq0gnmHJdptFj2dnoaWYWZYnLass6yx/snzCgn4ybY9R/lV+dTUDWBACAEAFxKUAHGR51nKud5efjp4bnZVjLMWstHzdcqBlb8sP/BxtQT8Y0dLghjXV+3fdfe6ICW7XkS0MAwMAEFKUJ4EU84LlYct/Wt7UA2VQd2+m5R7vLkMdwGGW3M5guZ623G+52LvrGfsAVtjiFevcoVc86t8Vxr6AAIDeYh/A/usA3mvRUXFPeHeVcZhFxeSLlpss21lKpUL4CMvuFh1fV0ijRZVekMEWVNBOIwa5i47dw7/rTp8oTgcBACCc/uoABvTx1dErR74O4EmWtZbFlp0tX7VoY7n9LMX2HFTlv8yiwq7T8hnLrZZCLrVckrl8Fx3Ayhs/V43Y4uZdMNWNGc6e4gCA8tAB7JsOoD7GTpZ8R7+VW/wV8hOLKobnLVosMt0y0dLTOOEai1YV7G/5skVnF6uTWMhlFn2xBBljQQTmHKVPX3GsCAYAoHeiLAAHWq63tFv+ZlERKCqyzs9cRma5Rd3AwmeMZejc4pctWgF8leXnlgsthaibqFcKQVRAIgKz9xvjmusLHxet1jErggEA6J0oC8CvWaZYtBp3gx7waY7dKZnLyIywjLWoECyH6goNB6Of6XSQO886wL/rrqZmgOvq6nKLWtuYDwgAQJmiLABnW86x/NaSPdT7Z0vPO/5urcWidk/Q8tE8P12Ps+h9V1q0nct4i4ZwNQz8D8svLYE7LBrCDajTN82yi2WSRV3JUy13WlAFVAQW2hamc0uX+8JPn3XHXzePoWAAAMoUZQGoVbjauDlXk6XcxSeTLc/4EQ0j6/orFi3eeK9FK4u1Avh2/60KwuwhWhWL2h8woL0Fb7CoIH3cooL1ny03W1AFVADee+4UN7hR63ny0xcwq4IBAChPlKuA51nusnzXokJsH8srFhVvKtjUfYs7Lf1lH8CIdWzudLtfpH2+i2NvQABAKVgFHG0H8EuWb1q+bdFs/n+1aGj2bMtFFqAkjXW1Pa4KPuWAcXQBAQAoUZQFoBZ7aD7eDhYNBZ9o0SpaLQz5gwUomVYFF3PXE0uYCwgAQImiKgA1aUubMy/132rRhxZbaJ6djlsDyqL5gOcdWXzt0D9/YJx76IU36AQCANCDKOcAav8/ra7VfnxJxRzAPrS8rd0deJmOkO4Z8wEBAIUwBzDaIeA/WvbOXALhldIFZINoAAB6FmUBqMUf2p/vTMu+Fs3izw5QtpP3H1dwb0A5eMK2bteRLQwDAwBQRJRDwDpmLVuwGbQ+pq4Ln/MVHwwB94Onl6x0s25Y4N8VNu+CqW7M8Gb/DgCADIaAo+0A7pGTPf0E10CvaCh4UEPPrx9YFQwAQH5RdgDTgA5gP9Hm0Dc88rK75rcv+490p70BT7YMH1RPJxAA8A46gNEWgB/z3xbyU/9tnFEA9qPnl7W5467VgTM9Y1UwACBAARhtAahtYLJpzE77A+rs3k2WJLRkKAD7kbaFmX71Y25Vu76c8tMX+PnTJrqZ79+RLiAAwEMBGO0cwKacNFp0HvBjluMsQCiaC3jvuVPcyBZ9aeWn1UZX/eZFbz4gK4MBAMiIsgDMpd/FiyxzLdfqASCsnUYMcvPmTnVzpk3wHymMRSEAAGT0ZQEY2GAZl7kEwmusq3WzJ491zfWFVwbrC33OtIl0AQEAMFHOATzKfxvQxxptOc/yliX3/XHEHMAqUur+gCwIAYB0Yw5gtB3AB3PygOVqy98tZ1iAitKcwGJdQL0CoQsIAEC0HcDcmfk6GaTwcs14ogNYZRavWOcOveJR/64wTgkBgPSiAxhtB1BDvfr3O/wExd9Ay79nLoHK0qKQ847czb8rjAUhAIA0i7IAvMwyOHO5lUEWvQ+IxMn7j2MoGACAIqIsAPV7Vlu/5NJZwG9nLoHK01zAO886wL/rjr0BAQBpF0UBuNzymiXY90/XQd6waOztlxYgMj0tCAkwFAwASKMoFoH8q0X/7g2WCyzZkys3Wl61JOW3LotAqpgWhOiouPWbdPpgfqccMM6dbBk+qJ5FIQCQEiwCiXYV8NGWhy1JW/mbjQKwypW6N6CwPyAApAMFYLRzAP/HEhR/+jgNOQEip6HgYU31/l1+LAoBAKRNlAWgtnu50rLEoqHf9pwAkVMBeO+5U9zIltxtKd/FohAAQNpEWQBebjnB8iWLCsBzLNr+RQtBPmUB+oT2Bpw3d6qbM22C/0hhLAoBAKRBlAXgTMunLXdaNAv/IctFli9bPmoB+kxjXa2bPXlsjyuDtSiELiAAIOmiLAC3tbyUufQmWA7PXDqd0zU1cwn0nZ72B5S7nljCUDAAIPGiLAD/bhmXuXR/tczKXHqrg1O54gb9r5z9ASkCAQBJFWUB+EPL/plL9w3LFywq/K63XGMB+pwKwAc+fzCbRAMAUi3KfQBz7Wr5gOVlyxN6IAHYBzCmetofUK+MvjBtopv5/h3ZIBoAEoZ9AKPrAGrjtQcs2csu/2b5kSUpxR9irKf9AbdYgq1hAABImqgKQG0AreFfbbEGVJ1gf8CehoJbGmvdQ3953b8DACAZopwDqG7fqZlLoPpof0DNByxmbUenO/OOp9yi1jYWhQAAEiPKOYBXWc60PG950rLOkk0bRMcdcwAT4NsP/Z+75iFNTe3ZvAumMicQAGKOOYDRdgA/ZPmLRR9Dm69p778gh1mAqnDy/uN6PC84wJxAAEAS9OUq4CSiA5gQPa0KFr2SYWUwAMQfHcBoO4CBMZZDLQO9O6AKaVHIyJZGt9t2Lf4j3bEyGACQFFEWgMMs91uWWB627GCRWyzaGBqoGioA582d6u741P5FVwarZT5n2kQWhAAAYi3KAlCLQJosEy3Zvy1/Zjk2cwlUj8a6WrfDsOaiK4O1r1HQBaQIBADEVZQF4HTLHEvu8soXLeMzl0D10fYw5x25m39XGEPBAIC4irIA1KqINZnLrQy3bMxcAtVJK4N72iR6+yGN7ubHXqETCACInSgLwHmWj2cuPcGpIOdbfpe5BKqT5gTeeZZ2LyrsjdUd7mv3v0AnEAAQO1EWgF+0fN7yS0uD5auWP1mOtsy1AFVNRWBPXUDRcXF3P9VKJxAAEBtR7wOoLWA+a9nPomLzacs1lqWWJGAfwIRbvGKdO+G6+W5Vu4637hknhQBA9WMfQDaCDosCMAVK2SQ6G0UgAFQ3CsBoh4BlsOVcy/WW6yznWKiUECsaCt6mucHV1Qxwgxvr/EcLY04gAKDaRdkB1FnAv7JstmjoV/a1aD7g8ZbH9UDM0QFMiY7Nna6rq8u99MZad/x18/1Hu9MrKo6LA4DqRgcw2g7gdy2/toyzHONH1yoK9T4gNrRJ9MD6Orft4EY3rKnef7Q7josDAMRBlAWgTgC53JK9559m0usYOL0PiB0NB9977hSOiwMAxFqUBeCzlnzHKUywLMpcAvGjk0I4Lg4AEGdRFoBXWr5j0SKQyX50/W3LNy3qAgYBYoXj4gAAcRblIhBNhyokOBVEH1/XPe+2W51YBJJiy9va3RFX/s6t39TpP9Ld2OFN7pLj93RH7jnKfwQA0N9YBBJtAbi7/7YU/+e/jRsKwJQrdY/Am0/djyIQAKoEBWC0Q8Aq6kpNTw6x3Gd5zaKO4QxLttssejw7Cy3FnGV5zPK2n4csxQ9/BXKUelzcmXc85Ra1tjEnEABQFaLsAMpIy4GW7Sy5xeb3/belmG6ZYtF+gr+wzLTcYwmoANzecrp3l6HVxyszl3ndZdGGbtqPcIPl3y2zLHtZlllKQQcQ3nFx069+rOhQcDZOCgGA/kUHMNoC8J8st1hU+KnDpq5cQNc7ZC7Lpr+brwAcZsntDJZDbRz9d2qhyh16oAQUgPCUe1zcq5cf618BAPoaBWC0Q8CXWa6wDLJo8tPorPS2+CvmMMublhctN1nUdSyHWjLa4bdY17DRokoviI66A7yh4GIbRGc75YBxDAUDAPpVlAWgCqTbLToKLmoPWE6xHG6ZY9nf8rBFBVuptGm1hn41F7CQCy1tWWm1AO9sED2ypdHttl2L/2h+dz2xhD0CAQD9KsohYO33t9TyLe+ucvINAedSl3Gx5WTL3XqgB5r/N9eiLuJzeqAAFZTZRaU6gK0MASOgM4NXrO1wR171+5LmBDIfEAD6HkPA0RaAdRat3NXH0MkfOgYu25f8t+UqpQCUlyw3W3T0XDH/ZrnIcqTlST1QBlV9zAFEN+UsDKEIBIC+RQEY7RDw+ZajLdoPUOdmTc2KOm1RGmEZa1nu3RX2RcvFlg9byi3+gIJ0UsidZ5W2qxCnhQAA+lqUBeAFlk9bdrZ80KLtYIJ8yFIOTap6nx/Rv6nrcRa9T8fO6d8db1Fxqc7jPyy/tAS0slcLUwIa9v2a5VOWVy1aqKIUn8AFlKjUhSE6LeShv7zu3wEAEL0oC0CNfVWqtaFzhJ/xI5pXqOuvWPRx3mu516IVwFp4orcqCNdYAioWNTcw8BlLg+XnFnUKg2hIGAgtWBjSUxG49O12b6NoFoUAAPpKlHMAL7FokUSSCyrmAKJH5cwH/NaJ+7gDdtmGOYEAECHmAEZbAP7YohM81FV73pK7CEQbRccdBSBKUu5pIWwUDQDRoQCMdghYtD/fnyzaC1DFZnaA1ChnUcj0vbZ3D73wBkPCAIDIUIiFQwcQJVve1u51AVe15zbDC6MTCACVRwcw+g6g/v2DLJ+0BKtrt7U0ZS6B9Mg+LWT7IT0fUjNn2kS6gACASETZARxjud8yyVJrmWh5xXKtRYXhOZa4owOIsgWnhRxzzbySuoFsFA0AlUUHMNoO4DWWFyzDLO16wKej2aZlLoH0aayrdTsMa3a3nq7djXrGRtEAgEqLsgA8xPKfluziT/5uUXcQSDUNCW/T3NBjG36bQfXu0CseYbNoAEDFRFkAFtr9dgfL2swlkF4qABd86XD387O1Z3lhK9dtcotXrPc2iwYAoBKiLAB/Y8me59dl0eIPbRD9oB4A0i4zHNzkmus1Tba4mgHOXfPQS2wRAwAILcpFIDp67VHLKst7LI9bdressxxs0QbRccciEFSENoo+4br5bBEDAH2ARSDRdgCXWPa23GS5w6IVwF+37GtJQvEHVIw2ii51UYiccsA4uoAAgF6LogN4q+U8yxrvLtnoAKJitFH0sdfMc6s3bHLDmxvcW2s7/PcUxhYxAFA+OoDRFIA67HS05U3vLtkoAFFR2iOwq6vLvfTGWnf8dfP9R4ujCASA8lAARjMEHOW8QiDRtChkYH2d23ZwY0lbxAj7BAIAyhXVHECt+AXQS6VuESMf3nN7N+uG+e65Vq23AgCgZ1EVgC9aVvYQAEUEW8QMayq0pWbGg395wz29ZJW7Y8Fi/xEAAIqLYrh2i+XzljbvrrDb/bdxxhxARE5bxEy/+jG3fpOm1xZ38bF7uAN2HuGGD6pnXiAAFMAcwOgKwFEWFoEAFVJOERhgn0AAyI8CMJohYOb/ARWmfQLvPOsA/65nLY21nB0MACiIVcBATGhhSE/zAQNrOzq9s4NvfuwVNowGAHRDsRYOQ8DoU705Mk7YKxAA3sUQcHSrgAFEQEPB9547xY1saXSTRg129bWlvYbTXoF0AgEAATqA4dABRL/QiSENtTXumaVvu1k3LPAf7RkLQwCADqDQAQRiSHsEDhgwoKx5gdsMqneHXvEIi0MAABSAQJypANSQsI6N62k0eOW6TW7xivXe4hAAQLpRAAIxp3mBOjbut3MOdc31tf6jhdVYoXjNQy+5h154g3mBAJBSzAEMhzmAqCq9WSXMvEAAacMcQDqAQKKoG3jr6ZP9u56xYTQApBMFIJAwmheoOYF1NQO87WKKCTaMZjgYANKFIeBwGAJGVdI2MV1dXe6lN9a646+b7z/aMzaMBpAGDAHTAQQSSdvEDKyvc9sObvS6gaW+0tOG0QCA5KMABBJMw8FaIfzzsw/0HymOvQIBIB0oAIGEUzdwh2GlbRidvVfgc62r/EcBAElDAQikQPaG0aUOB3/ilifczY+9wuIQAEggFoGEwyIQxIoWhyxf1e5mXP94WXsF3nzqfu7IPUf5dwAQbywCoQMIpIqGg8dv21LWXoGiIWE6gQCQHBSAQAqVs1dgQCuEP3z171kgAgAJwBBwOAwBI7Z6u1egcHwcgDhjCJgOIJBauXsF1pb4crDG/twtj73iFrW2MSwMADFFBzAcOoBIhN4uDhG6gQDihg4gHUAAJlgcoq1iStkvMNDSWMucQACIIQpAAO/YacSgd/YLLGVIeG1Hp7dCmP0CASBeGAIOhyFgJJKGhJ9f1uY++t0F/iOlmXfBVDdmeLN/BwDViSFgOoAA8giOj9MWMZNGDXb1Ja4Q0VYxdAIBoPrRAQyHDiASTZ3Ahtoa98zSt92sG0rvBk4a1eK+OXsft/eYYf4jAFA96ADSAQRQhDqBAwYM8DaOLmdxyF9fX+vOuO1J91zrKv8RAEA1oQAE0CMVgOUsDpG31na4T9zyBAtEAKAKMQQcDkPASJUw+wWyQARAtWAImA4ggDL0dr9AYYEIAFQPCkAAZQv2C+xNEcjG0QDQ/ygAAfRKUARqq5jdtmvxH+2ZNo6edcN8FogAQD9iDmA4zAFE6mle4Iq1He6478x3be0bXWeX/44eDG6sdTsOb2K7GAB9jjmAdAABhJTZNLrZLfjS4e63cw4teVh4TUent13MTb9nlTAA9DUKQAAV0dsFIvc9t9ybG8h2MQDQd+JSAB5iuc/ymkUDTDMs2W6z6PHsLLQUs5flF5ZXLfrzn7cACKm3C0S+dv8LXiEIAIheXArAQZZnLed6d/k9aBmdlWMsxWhDslcscy0sSwQqqLdFYFN9rfvKfX9mkQgARCyOi0DUrZtpuce7y1AHULPIczuDpVIX8Go/5WARCFDE4hXr3MzrH3erN2xyw5sbvNNBSqXVxbecNpkFIgAqjkUgyZoDeJjlTcuLlpss21kA9CN1ArU45Pn/PMrde+6HyuoIqlg85aaF7sirfsfegQBQYUkpAB+wnGI53DLHsr/lYUujpZL076nVF2SwBUARWhwysL7OWymsYeFyzhPWSuGX31rr7R1IEQgAlZOUAvAnlvstz1u0WGS6ZaLlWEslXWhpy0qrBUCJgo7gTz99oP9I6VQEslIYACojSUPA2ZZbFlsmeHeVc5lF8wWCjLEAKENm38CmsjqBgWClMEUgAIST1AJwhGWsRYVgJWkGuyaLBlljAVCm0UObyt44OpuKQFYKA0DvxaUA1EGj7/MjO1t0Pc6i911p0ZjSeIsWg2gY+B+WX1oCd1jUwQs0WIJ/U9c7+te7WQBELHvjaK34nTRqsLcNTKmeXrLKnXHbkxSBANALcdkGRkVdvh1ib7ecbdGWMPtatF+Eun76sxdblloCj1q03ctp3l2mWPx75nIrv7Po45WCbWCACtB5wg21NW7JyvXuhOvmuVXtm/339ExnCm8/pMnNnb67O3LPUf6jAFAY28DEcx/AakIBCFRYsHdgW/tG16ldP8swaVSL++bsfdg7EEBRFIAUgGFRAAIRUEdw+ap2N8MKwVXtm/xHSzO8ud5dMXtvuoEACqIApAAMiwIQiJC6gSdcN7/sIlA+NWW8+9PSVe7Sj+xFRxDAVigAk7sKGEAC9PZMYbl1/qveQpE7FmhHKABANjqA4dABBPqAOoGzv7vAjWhpsOv1rn1Tp/+e0owa2uhaGupZKALAQweQAjAsCkCgj4RZKZztV+dOYUgYSDkKQIaAAcSE9g0cMGCAPyx8UK9OEpGTv7/QHXrFI5wtDCDV6ACGQwcQ6CfBSuFZNyzo1ZYx8v5xw1gkAqQQHUAKwLAoAIF+FmbLGNFG0jsOb2L/QCBFKAAZAgYQc2GPlFvT0en++vpad/oP/siwMIDUoAMYDh1AoIoEC0WeWfq2NzTcG1oxrG1n6AgCyUUHkA4ggAQJFoqMHtr0ziKRcl/lvt7W4XUE2T8QQJLRAQyHDiBQpdQN7OrqsoJug5txvU4TKX/bmI/sM9o929rmmupr6AgCCUIHkAIwLApAIAa0kfTM6x93qzdscsObG9xbazv895Ru1JCB7vun7kcRCCQABSAFYFgUgEBMBB3Bles2umOu0UbS5a8YHtRY67ZtaXQXH7sHJ4oAMUYBSAEYFgUgEENBR7C3+wfKpFEtDAsDMUUBSAEYFgUgEFNh9w+U4c317pzDdnW//etb7sJjJlEMAjFBAUgBGBYFIBBz2d3ALV3O9bIh6G0dM3abZvf1me+hEASqHAUgBWBYFIBAAgTzA//82mr30e/2bv/AwIhBDe4bH30vcwSBKkYBSAEYFgUgkCDL29rdR66d70a0NLglK9e79Rs7/feUjw2lgepFAUgBGBYFIJAwwWkiGzu3uNfeXu9m3qA5guXvIRg4fNJ2rt0KSeYIAtWDApACMCwKQCDhKrFiWAY31rodhzfREQSqAAUgBWBYFIBACgQrhnW+cNhCcOjAOjd62EAKQaAfUQBSAIZFAQikSFAInnjjQm+eoLqD7Zu2+O8tjzqC29i/MWRgAyuHgT5GAUgBGBYFIJBCwTxBLRSZef1819a+KVRXkHmCQN+iAKQADIsCEEi5oCsYZkPpAF1BoG9QADpX478FAPRCY12tG79ti7v33Clum+YGV2svq2ssdb346bqmo9MtXtHuFi1rc2fc9qR7rnWV/x4AqCw6gOHQAQTwDnUDtaG0rFjb4Y79zrxQW8g0N9S6IU117CcIVBgdQArAsCgAARRUqS1kRPME31rT4V0zPAyEQwFIARgWBSCAoiq5hUyAc4eBcCgAKQDDogAEUJLcQlB1oOYKbu7dLjIeFo0AvUMBSAEYFgUggLJUep5ggCFioHQUgBSAYVEAAggle57gFqsLKzBC7A0Ra6Pq5oY6ikEgDwpACsCwKAABhBZ0BV9v2+BmXD+/Ih3BwMiWRnfLaZMpAoEsFIAUgGFRAAKoqOyOYCXmCYq2kxk5uNF94gPj3L3PLvceozOINKMApAAMiwIQQMXlzhM87tr5bnXI4+ZyjRoy0H3/1P0oApFKFIAUgGFRAAKIXBRbycigxlo3eGCda6yrYSUxUoUCkAIwLApAAH0miq1kstEVRFpQAFIAhkUBCKDPRTlEvOfowW7dxk676nJ1NZkDjVlNjKShAKQADIsCEEC/y+0MVmo7mWx0B5EkFIAUgGFRAAKoGkFnMIrtZETzBU/cb4yb//I/XENdLV1BxBYFIAVgWBSAAKpSsJ3M6g2b3LYtje6N1Rsq3hUcMajB7TCsiUIQsUMBSAEYFgUggKoVdAQb62q96yVWFH7sewsr3xnkTGLEDAWgc5kZvgCAxFHhN7C+zg0YMMB7O3HUUHfvuQe5bZobXK29/K9UB2BNR6dbvKLdLVrW5s647Un3i6eWuo9/f6F7rnWV/ycAVBs6gOHQAQQQO9mriLV45MTvLXRDm+rd395a5z0W1vaDG90bazreOZNYWEmMakIHkAIwLApAALGngtDbTuY78yNbRSyaMzisud4rBk/70E7u508tcxceM4miEH2OApACMCwKQACJkd0Z1NulK9dHMmdQdhvZ4l5+ay0LSdAvKAApAMOiAASQaFpNPPu7C7yh3CVWEK73NomuvKED69ywQfXe5tMMFyNqFIAUgGFRAAJIPHUGG2pr3MbOLe61t9e7WVYQVurkkUK08fQXj57IMDEiQQFIARgWBSCA1FFBGPXJI8JiEkSFApACMCwKQACpFcwZ1MkjJ964wNsLcMmK9a590xb/T0QjmDeohSS3Pb7Ye4yiEOWgAKQADIsCEABMMEyseYIzr9dq4k2RdQYDY4Y3uda3273r7BXGFIPoCQUgBWBYFIAAkCN7NbE3VOzPGdQjNfZbZ3O0DcJ3howpBlEIBSAFYFgUgADQg+yCUF1Cbf9y0vcWRr6QRA6ftJ17a02HW7/x3a1sKAxBAUgBGBYFIAD0Ql8tJCmE/QfTjQIwPgXgIZYvWvazjLbMtNxjCdxm+WTm8h1/sHwwc1nQRy1ftexq+Zvly5ZfWkpFAQgAIeQOF+tYOm/PwT5YTCLDm+vdOYft6u59drnXJaQ7mA4UgPEpAKdbplietvzCkq8A3N5yuneXsdGyMnOZ14GWxywXW1T06d/8iuUgi4rHUlAAAkAF9cdiklzBghKhIEwmCsB4DgHr50C+AlDfnTO8u9L8xKKqTcVl4EHL25aPe3c9owAEgIjkW0yyZsNmr0P45uqOPisKj9t7tFu8QqegZOYRUhTGHwVgsgpAFX/q+q2y/M6i4dw3LYUssXzbT+ALls9bdvLuekYBCAB9JCgIG+tqveslK9a5k77/B28xSZRdwrqaAW6zPkCWfF1Cuey//8rJJTFAAZicAvAky1qLdgTd2aJ5fXUWzRnssOSjYvE0y4+8u4x/svzA0ujddafHs9832NJKAQgA/SMoCrUZ9awbHveGjPVLoi+2m8mmLuGTr77tXl+9gT0JY4ACMDkFYC4tFFExeLLlbj2QhwpALRz5sXeXcYrlFstA7667Sy2XZC7fRQEIAP0ve8h4xdoOd9y18yPvDgZq7bdpvi1tsvckDE4uYSi5/1EAJrcAlJcsN1u+4d1115shYDqAABAT2QWh3i5duf6dIeOo9x/MZ7eRLd4eiNmyu4Ucbdd3KACTWwCOsCyz/IvlDj2QhxaBqIA7xrvLeMCiOYQsAgGABFJR6G03c6O/3czKdW79xj4cKy4iu0DMnmMYCIrEnz+1jHmGIVEAxqcAbLHslrl0z1jOtzxi0TYvioZmtT3Mcst4y39Zxln2sKyxiApBFYUXenfOfcjye4sWi9xrOcHyNQvbwABAwqkQ1HYzGzu3uC1btrilb7e7k7+3sF+2nSnHHqMGuxdeX9NtEQqFYXkoAONTAB5mUcGX63bL2RZ1A/e16KteRaD+rPb3W2oJPGp51aKFH4HZFhV9u1iCjaALzRnMhwIQABIid8jYW1jin2NczUWhNDfUuvUbO92oIQPdF4+euNVcw2z55iLmqsSfqeTHiWI4nAIwnkPA1YQCEAASLLsoDPYirPaCcMzwJtf6drt/113QRSymEn+mUh/ntA+Nd5d+ZC//rjIoACkAw6IABIAUydcl1PF1wwfVu2VWdFXLfMIkaaqvdf929ES304hBbpIVjGOGN/vv6T0KQArAsCgAASDlis0nVKnY13sSJt2rlx/rX/UeBSAFYFgUgACAbrI7hSoO13Zs9u4HDBjg7VEYDCVTIJZO/z9962PvczP23dF/pPcoACkAw6IABACULV+B+Hpbe58cbRdXv/7sQe49O3pFW2gUgBSAYVEAAgAqJneO4Ua/NdhQV+MNLZ/0vYWpLRApACuLAjAcCkAAQJ/JLhC9Da39BSjLV21wmzu3uA2bk1UWDrAqpcb+Z2hTvbv/cwe50UOb/PeEQwFIARgWBSAAoN9kL0Cprxng1m3sdHU1zm3q7HL1tQPe6SCKuoh6PJiLmPt+qcSfqeTH0ftEbxvrar3rSqAApAAMiwIQAICYoQB0zl4nAAAAIE0oAAEAAFKGAhAAACBlKAABAABShgIQAAAgZSgAAQAAUoYCEAAAIGUoAAEAAFKGAhAAACBlKAABAABShqPgwvGOglu6dClHwQEAEBM6Cm7s2LG65Cxg9MqOltbMJQAAiJkxlmWZy3ShAAxH///tYFnj3VXWYIuKS31xRvHv97ekPz/hcxh/fA7jj89h/EX1HPXvvmbp8u6AKqExZX1RJnVsOenPT/gcxh+fw/jjcxh/aXiOfY5FIAAAAClDAQgAAJAytf5bVKdOy6P+2yRK+vMTPofxx+cw/vgcxl8aniMAAAAAAAAAAAAAAAAAAAAAAACAOPiM5e+WDZanLAdb4uoQy32WYLf1GZZsOk3lUove327RKq+9LHFxoeWPFu1O/6blHsvulmyNlmst/7Css/zKoh3t4+Bsy3MWnZWpLLBMtwTi/NwK0edUX6tXe3cZcX6e+v7S88nO65ZA3L8HAzqa807LCst6y58s+1kCcX6er1pyP4fK9RZJwvdhneVrFv3u0+fnFct/WLK3q0vK1yqQ10mWjZYzLXtY9EtorWWcJY5ULOibepZFP7ByC8ALLCos9P73WP6fRd/cOqInDh60nGbRD6F9LL+2LLYMsgS+a9ExRkda9rU8bNEvpzhsw3S85RjLRD9ft+jrM/ihG+fnls/+Fv0CetaSXQDG+XnqF+bzllFZGWkJxP17UIZbVCT9wHKAZbzlCMuulkCcn6c+X9mfP30d6ufpYRZJwvfhly0qYI+16PM326IX1udZAkn4WgUK+oNF38zZXrBclrmMtdwCUK/mllv0TR3QK9lVln/17uJHP6j1PNX5lKEWFUwq7AM6P1p7WR3t3cXPSssZlqQ9txbLixb9ElVnISgA4/48VQCqGMgnKd+Dl1sey1zmlbSfNfrafNmi55WU70O9eL4lc/mOX1h+mLlM5O+LfsVJINWlwaIhi//17t6l+w9lLhNlZ4tezWY/3w7L7yxxfb76YSwqkkSfz3pL9nPUK1Z1ZOL2HNVNONmi7qaGgpP03ETDafdbHvLu3pWE5znBov9mdTfVNdnFIkn5HvyI5UnLzyyaivGM5SxLIEk/a/R74p8tt1r0YjMp34fzLOraaqRBNKJykOW/vbtk/r7oVxSA1WVbi37JvuHdvUv3+sJPmuA5JeX56hXqtyz6QaYfvqLnoVfnb3t374rTc3yvRdMQ9MP2RstMy18sSXhuARW2+kWq+X+54v48NapwqkXdIBVF+m9+3DLCv5a4fw+qoNV81Zcsep76Ov2ORc9bkvI8RaMowyy3eXeZ//4kfB9+w/Jjy18tmywq4tXp1GOSpM9hVaAArE56VZdNhUXuY0mSlOd7nWVvy8e9u+Li9Bz/z/I+ywctmp5wu2VPSyFx+/yNtVxjOcWihVelisvzfMCiobRFFnU3NcdKPum/lbh/D+p32dOWL1lUOHzPcpNFRWG2uD9P0fQLfU7V5Ssmbs9NQ9jqbP6T5f0WfX3+m/82WxI+h1WBArC6aAKs5m3kvprZzpL7qicJgpWISXi+WoGnYaipFk3GDug5ashGk9Szxek5qrug+UYaYlOHTAskNDE7Cc9N1PnTf7NW3G/2c6jlc/61nksSnmdAq0RVDGpYOCnfg5obpq50Ns2dDhbPJeV57mTRHNWbvbuMpHwfXmHRXE5NUdDXp+b+fdsSdOWT8jmsGhSA1UW/aPVLaJp39y7da8gmaTQfSd/U2c9XP8j0yzcuz1evPtX506q0wy16Ttn0+dRwRvZzHG3RCra4fk71nDX5OinP7bcWDXOryxlExe5dWddJ+hzqc6cdBlQ0JeF7UOZbcrdf0lwyrciXpDzP0y2a46i5qoGkfB82W7ZkLt+hhkhQpyTlcwgUpDa4CsFPWfRDWq+ANP9Kr/ziSCsrg1+qatN/wb8OXplrRZdWcWlemX5g/ciioY24LOu/waL/fv0Q0ivTIE2WgIZNl1o0wVlbNKjgiMsWDf9l0T6U2pZBRZK2gdEP5eCHcJyfWzHZq4Alzs/zSou+PjWJ/gMW7cuprTSCnylx/x4Ubd+jIkhDwLtZNIyoTqeG9QNxf54qhFTQqkuWKwnfh5rTqNGTYBsYfZ7esmhuYCAJX6tAUdoIWntaadK9Xt0FW4rEkfapUuGXm2ACs7pJ2qZC3QjNv9KKLn1jx0W+56Zob8DAQIuGiIMNavULWPPO4kDbMgRfi+o8aA5Z9ivwOD+3YnILwDg/z2CvNL2wXGbRfMDsOZxx/x4MHGfR0KGeg4Z/s1cBS9yf51EW/WwJVslmS8L3oYo4fc+pyNUmz3+zaA9ZdfkCSflaBQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIDE0IkC2jxXp88AAAAg5rJPXskXnTajo7B0NF+dBQAAADGXffbyeZa2nMeGWgAAAJBQOntZh8Tnyh0CDs6mPtryjEVnjj5s2c4y3aIzZFdbfmxptgR0Bum/W16x6O88a5ltKcc+lkcsayz6GDrne7IFAAAAvVBuAbjAMsWyr+Uly6OW//HvD7b8w3KBJfB1i4pDFY67WPTxdBj9oZZSPW/5oWWSZYLlRIuKQgAAAPRCuQXgEd5dxlyLHlNhF7jR8mDm0g2yqOt3oHf3rpstP8pclkRdv09mLgGgf9X4bwEgTZ7z38oblvUWDe8G9JiGhWVPy0DLbyxrs3KqZVdLqb5lUdH4kEVFZzl/FwAqigIQQBpt8t+Kun/Z96LHgp+PwdtjLeokBlFhWM48wEste1nutxxu+YtlpgUA+hwFIAAUp0KtwzLO8nJOllrK8aLl25ajLHdbTrcAQJ+jAASA4rRq90qLCjfN4dPQrRaLnOPfl6LJcp1FcxB3smgByv4WLSwBgD5HAQgAPbvY8hXLhRYVbVoxfLzl75bAqxYN8+bTaRlhucOiLuBPLQ9YLrEAAAAghtTh00rhqd4dAAAAEu/DFi3uAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAQBI59/8BxVyQO/RZzX0AAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "lake.heater_1.control_mode(\"Manual PID\")\n",
    "lake.heater_1.output_range(\"Low (2.5W)\")\n",
    "lake.heater_1.input_channel(\"A\")\n",
    "# The following seem to be good settings for our setup\n",
    "lake.heater_1.P(400)\n",
    "lake.heater_1.I(40)\n",
    "lake.heater_1.D(10)\n",
    "\n",
    "\n",
    "lake.heater_1.setpoint(15.0)  # <- temperature\n",
    "live_plot_temperature_reading(lake.sensor_a, n_reads=400)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Querying the resistance and heater output "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "25"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# to get the resistance of the system (25 or 50 Ohm)\n",
    "lake.heater_1.resistance()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "50"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# to set the resistance of the system (25 or 50 Ohm)\n",
    "lake.heater_1.resistance(50)\n",
    "lake.heater_1.resistance()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.0"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# output in percent (%) of current or power, depending on setting, which can be queried by lake.heater_1.output_metric()\n",
    "lake.heater_1.heater_output() # in %, 50 means 50%"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.8"
  },
  "nbsphinx": {
   "execute": "never"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}